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Introduction 


The scientific goal of the experiment is to test the equality of gravitational and inertial 
mass (i.e., to test the Principle of Equivalence) by measuring the independence of the rate 
of fall of bodies from the composition of the falling body. The measurement is 
accomplished by measuring the relative displacement (or equivalently acceleration) of 
two falling bodies of different materials which are the proof masses of a differential 
accelerometer. The goal of the experiment is to measure the Eotvos ratio 6g/g 
(differential acceleration/common acceleration) with an accuracy goal of several parts in 
10 15 . The estimated accuracy is about two orders of magnitude better than the present 
state of the art. 

The main goal of the study to be carried out under this grant is part of the flight definition 
of the experiment and laboratory testing of key components. The project involves an 
international cooperation in which the responsibility of the US side is the flight definition 
of the experimental facility while the responsibility of the non-US partners is the flight 
definition and laboratory prototyping of the differential acceleration detector. 

In summary, the experiment to be designed is for taking differential acceleration 
measurements with a high-scnsitivity detector (the sensor) during free fall conditions 
lasting up to 30 s in a disturbance-free acceleration environment. The experiment 
strategy consists in letting the sensor free fall inside a few meters long (in the vertical 
direction) evacuated capsule that is falling simultaneously in the rarefied atmosphere 
after release from a helium balloon flying at a stratospheric altitude. 

Description of Experiment Concept 
Need for Picogravity Environment 

The accuracy of the Weak Equivalence Principle (WEP) tests with laboratory proof 
masses on the ground is limited by the Earth's seismic noise and the weakness of suitable 
signal sources. Previous experiments include the famous torsion balance experiments of 
Eotvos (1 890- 1922) 1 " as well as the classical tests of the Equivalence Principle by Roll- 
Krotkov-Dicke (1964) Hi utilizing a torsion balance which had an accuracy of 3 parts in 
10" and the 1. 1. Shapiro et al. (1976) iv and Williams et al. (1976) v lunar laser ranging 
experiment with an accuracy of 3 parts in 10 12 . The present state of the art is at a few 
parts in 10" both for lunar laser ranging tests (Williams et al., 1996) VI and for torsion 
balance tests (Adelberger et al., 1999)™. 

By conducting the experiment in free fall, the signal strength increases by about three 
orders of magnitude because the full strength of the Earth’s gravity is sensed in free fall. 
Seismic noise is also absent in free fall. Seismic noise is replaced in orbit by the noise 
sources of the space environment which require drag free stages or drag compensation in 
order to achieve the promised improvements in the test accuracy. An orbital free fall, 
with a "drag-free" satellite, is one approach that has been under development for many 
years. A small number of orbital tests of the WEP has been proposed with estimated 
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accuracies of 10 1,11 , 10 lx , and 10 1 \ An alternative to the free fall in space is the 
vertical free fall inside a drag-shielding capsule released from a balloon at a stratospheric 
altitude as proposed in this experiment* . 

An orbital free-fall has advantages and disadvantages with respect to vertical free fall. On 
the one hand, orbital free-fall tests can achieve an even higher accuracy than vertical free- 
fall tests thanks to the longer integration time and lower resonance frequency of the 
detector. On the other hand, vertical free fall tests have some key advantages over orbital 
tests. First of all, in a vertical free fall (from a balloon) the experiment can be repeated at 
relatively short intervals of time (a few weeks) and at a more affordable cost. The ability 
to repeat the experiment is important for the success because these high-accuracy 
differential detectors can not be tested on the ground at the accuracy that they can achieve 
in free fall conditions. Therefore, modifications and improvements have to be expected 
before the detector/experiment performs at the estimated free-fall accuracy. 

Both orbital and vertical free fall are Galilean experiments in which the differential 
displacement or rate of fall or acceleration is measured between two bodies of different 
materials falling in a gravitational field. However, classic Galilean experiments, in which 
the relative displacement of two bodies falling side by side is measured (with drops 
ranging from lm to 140m) have yielded an accuracy* 1 in testing the WEP of order 10‘ 10 . 
The limitation mostly stems from relative errors in initial conditions at release which 
propagate over time due to gravity gradients. This problem can be overcome in orbital 
and long vertical free falls (i.e., from stratospheric heights) thanks to two provisions: (1) 
the initial relative motion of the two sensing masses inside the detector is abated during a 
damping phase preceding the measurement phase and (2) the detector is rotated with 
respect to the gravity field in order to modulate the signal (at a frequency f s ) and move 
the frequency of key gravity gradient components to 2f s . 

The test of the Equivalence Principle requires a differential measurement of acceleration. 
This fact has a positive consequence in terms of the rejection of accelerations that affect 
the two proof masses equally (common-mode type) and their effects on the differential 
acceleration. Typical values of the common-mode rejection factor of differential 
accelerometers are of order 10' 4 . Consequently, for an experiment that aims at measuring 
differential acceleration of order 10‘ 15 g, the common-mode acceleration perturbations 
external to the detector must be of order 10' 12 g or less. 

Drop Facility 

The following is a preliminary description of the drop facility the design of which will 
evolve as a result of the analyses carried out during the flight definition phase. The free 
fall facility (see Fig. 1) consists of: (1) the gondola that stays attached to the balloon; (2) 
a leveling mechanism that keeps the capsule vertical before release; (3) the capsule, 
which houses a large vacuum chamber/cryostat; (4) the instrument package which free 
falls inside the cryostat and contains a small, high-vacuum chamber which in turn houses 
the detector; and (5) the parachute system to decelerate the capsule at the end of the free 
fall run. 
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The capsule is kept vertically leveled and stabilized in azimuth by the gondola before 
release. Upon reaching an altitude higher than 40 km, the capsule is released from the 
gondola and immediately afterwards (< 1 s) the instrument package is released from the 
top of the capsule. The analysis indicates (see later on) that with a 1-3 m long vertical 
space available inside the capsule, the instrument package will span that space in 25-30 s 
while the capsule, that is slightly decelerated by the rarefied atmosphere, falls by a few 
km over the same time. The capsule shields the instrument package from external 
perturbations and allows it to free fall under acceleration conditions which are close to 
ideal. The differential acceleration between the two falling test masses is measured 
during the free fall time. At the end of the free-fall run the capsule is decelerated by a 
parachute system for recovery in water or, alternatively, over land. 


Small 

Cryostat 


Spin axis 

Sensitive axis 

Large 

vacuum 

chamber 



Detector 

release 

mechanism 


Instrument 
package with 
detector 


Figure 1 Pre-definition-study configuration of capsule in free fall after detector release 

Detector Concept 

The following is a preliminary description of a differential acceleration detector concept, 
the design of which (carried out in cooperation with our non-US partners) will evolve as 
a result of the analyses conducted during the flight definition phase. 

The detector that we plan to use for the experiment is a differential accelerometer 
that will be developed at the Institute of Space Physics (IFSI) in Rome (Italy), under the 
sponsorship of the Italian Space Agency in the framework of the participation in this 
project of non-US investigators (with V. Iafolla, PI). This detector technology xl " has 
been pioneered by V. Iafolla and the late F. Fuligni and applied to the construction of a 
number of high-sensitivity, low-frequency accelerometers over several years. In the 
following we give a brief description of the detector conceptual design at this stage of the 
project. 
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Figure 2 Longitudinal section of conceptual instrument and sensing masses. 

The two sensing masses are constrained by torsion springs to rotate independently 
about the twist axis (which is parallel to the spin axis of the instrument) and their 
resonant frequencies are electrostatically controlled for frequency matching. The 
displacements generated by the rotations are sensed by the capacitive pick-ups of the 
instrument as explained later on. Sensing mass 1 (in dark color) is a hollow cylinder 
mostly made of a given material while sensing mass 2 (in light color) is a dumbbell- 
shaped cylinder made of a different material. Each sensing mass constitutes the moving 
part of a capacitor w ith symmetric fixed plates on cither side of the sensing mass (see 
Fig. 3). Capacitor 1 is formed by sensing mass 1 and the fixed plates marked A and B 
while capacitor 2 is formed by sensing mass 2 and the fixed plates marked C and D. The 
fixed plates A and C are used for signal pick-up and the fixed plates B and D for feed- 
back control. The displacement of sensing mass 1, for example, is detected by the series 
capacitances As (one fixed plate on each side of the sensing mass). These plates form one 
branch of a capacitive bridge in which two additional reference capacitors form the other 
branch. The bridge is pumped by a quartz oscillator at a stable frequency of 10-20 kHz. 
reducing the relevant noise temperature of the preamplifier. The difference between the 
output signals from capacitors 1 and 2 is amplified by a low-noise preamplifier, sent to a 
lock-in amplifier for phase-detection, and then to a low-pass filter. 

The cross sections of the ellipsoids of inertia about the spin axis of the instrument are 
circular so as to minimize, within the construction tolerance, the mass-moment torques x,v . 
In the detector shown in Figs. 2 and 3. the inner cylindrical mass is made mostly of a 
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The differential-acceleration detector (see Fig. 2) measures the relative displacement, 
along the sensitive axis, between two sensing masses of different materials. The centers 
of mass of the sensing masses are made to coincide within the attainable values in order 
to minimize the effect of gravity gradients, rotational motions and linear accelerations 
upon the differential output signal. 


high-density material (e.g., Platinum-Iridium) while the outer dumbbell-mass is made of 
a low-density material (e.g.. Aluminum). 



Figure 3 Interior of conceptual differential acceleration detector. 

In order to achieve an experimental accuracy of several parts in 10 15 in 25-30 s 
integration time, the detector must have an intrinsic noise (expressed in terms of 
acceleration) of about 10 14 g/Hz l/2 . Earlier analyses indicate that this level of noise is 
attainable with an instrument refrigerated to a temperature close to that of liquid Helium 
and with state-of-the-art low-noise preamplifiers. 

Experiment Sequence and Communication Links 

Figure 4 shows the preliminary timetable of the experiment sequence. The 
experiment starts with the loading of the sensor into the vacuum chamber/cryostat about 
2 weeks before the planned launch. This operation is then followed by the pumping 
down of the chamber and the refrigeration of the sensor. After connecting the capsule to 
the gondola and the balloon, the balloon is launched. The estimated time to reach altitude 
is of order 3 hours. Upon reaching altitude, the attitude of the capsule is stabilized by the 
leveling mechanism attached to the gondola, the sensor is spun up, and the dynamics of 
the system is analyzed. When the dynamics is within the acceptable bounds, the capsule 
is released from the gondola and the sensor is released from the top of the 
chamber/cryostat immediately afterwards. The science data is taken during the free-fall 
phase in which the sensor spans the length of the chamber. Shortly after the sensor has 
reached the bottom of the capsule, the blut (first stage of the deceleration system) is 
released and. when the speed has decreased below the required value, the parachute is 
deployed. The chamber is vented before the capsule hits the surface/water and the 
locator beacon is turned on. 
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Figure 4 Preliminary timeline of experiment 

Figure 5 is a schematic of the communication system between the ground, the 
gondola and the capsule through radio links and the communication between the sensor 
(during free fall) and the capsule through an infra-red link. 
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Analysis of Free Fall Phase 


The free fall time of the instrument package inside the capsule can be computed after 
considering that the capsule is slowly decelerated by the air drag while the instrument 
package (after release) moves inside the vacuum chamber at low relative speed and 
consequently it is unaffected by air drag (it is indeed in free-fall conditions). 

Free fall time 

The free fall time and vertical size of the vacuum chamber/cryostat can be computed 
from the equations of motion of the instrument package in free fall and the capsule in 
decelerated fall. The equations of motion are as follows: 


2f = g 

» = s--!-c d spk (1) 

2 fti^ 

where z is the vertical distance from the time of release (the subscript 1 stands for 
instrument package and 2 for capsule), S is the frontal cross section of the capsule, C D the 
air drag coefficient of the capsule and p = f(h) is the air density with h the altitude above 
the Earth’s surface. Equations (1) can be solved analytically only if C D is assumed 
constant and the atmospheric density exponential. We will not spend time on the analytic 
solution because it is valid only for relatively-short drops. 



Figure 6 Atmospheric density in the stratosphere per US Standard Atmosphere 1976 
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After adopting a density profile from the US Standard 1976 Atmospheric Model (as 
shown in Fig. 6) the exponential fit of the density for the altitudes of interest is: 


P = P o* 


-h/H 


( 2 ) 


Where H is the scale height, p 0 is the reference density which is taken at the Earth’s 
surface and h the altitude above the Earth’s surface. The relative distance can be 
obtained as the double integral of the acceleration difference between the capsule and the 
instrument package d$<= which is equal to the deceleration of the capsule due to 

aerodynamic drag: 



(3) 


where (5 = m/(C D S) is the frontal ballistic coefficient of the capsule, h D is the drop 
altitude, z * (l/2)gt 2 the distance traveled by the capsule and g the Earth’s gravitational 
acceleration. The drop velocity of the capsule is assumed equal to the free fall velocity gt 
only for the purpose of computing the air drag deceleration (which is a valid 
approximation at high altitudes). 

The air drag coefficient C D is fairly constant in the non-compressible regime but then 
it grows substantially with the Mach number for speeds approaching the transonic regime 
as shown in Fig. 7 for an aerodynamically-shaped cylinder with fineness ratio D/L = 6 * 
0.25. The Mach number M is the ratio between the actual speed of the capsule and the 
speed of sound at the local altitude: 


M 


V 

4r 


(4) 


where R = 287 J/(kg-K) is the gas constant of air, T the local air temperature and y * 14. 
The Mach number vs. the drop time is shown in Fig. 8 together with the drop distance vs. 
time. 

As a result of the functional dependence C D = f(Mach), we can separate the ballistic 
coefficient into two components as follows: 

ft = Ad + AjSS (5) 


where p (l = m/(C D0 A) is the low-speed ballistic coefficient (i.e., its minimum value) and 
Aft is the fractional variation of the ballistic coefficient due to the increase of the Mach 
number. 
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Figure 7 C D vs. Mach number for an aerodynamically-shaped cylinder with 6 = 0.25 (see 
text) 



Figure 8 Mach number and capsule drop distance vs. time 

After taking into account the functional dependence C D = f(Mach) and integrating 
eqn. (3) twice for different values of the ballistic coefficient, we obtain the numerical 
results shown in Fig. 9, for a fineness ratio 6 ■» 0.25, a drop altitude of 40 km and a (low- 
speed) ballistic coefficient |} 0 ranging from 2000 kg/m 2 to 10000 kg/m 2 . 
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The numerical results clearly indicate that it is possible to attain a free fall time between 
25 s and 30 s with very reasonable lengths of the chamber and capsule. Appropriate 
values of the (low-speed) ballistic coefficient in the range of greater interest of 6000- 
10000 kg/m 2 can be readily obtained with capsule masses si 500 kg and external 
diameters smaller than 1 .8 m. Designs options will be investigated later on in this report. 



Figure 9 Relative distance for various initial ballistic coefficients 
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Figure 10 Capsule deceleration due to air drag (p„ = 6000 kg/m 2 ) 
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Altitude (km) 

Figure 11 Limit velocity and actual velocity of capsule vs. altitude ((5 0 = 6000 kg/m 2 ) 

For completeness the deceleration of the capsule caused by air drag is shown in Fig. 10 
for p o = 6000 kg/m 2 . Figure 1 1 depicts the limit velocity and the capsule actual velocity 
vs. altitude for the same case. The limit velocity is the velocity at which the gravitational 
and air drag force are equal and, consequently, for a limit velocity much greater than the 
actual velocity the air-drag deceleration is very small. 

Effect of Wind Shear 

The horizontal velocities of the capsule and the instrument package (attached to the 
capsule) are the same at the start of the fall. The inertial horizontal velocity is determined 
by the rotational velocity of the Earth at the latitude of capsule release and by the local 
wind. The former (which is much bigger than the latter) simply makes the falling bodies 
follow a parabolic trajectory rather than a fall along the local vertical. The maximum 
lateral displacement is of order a couple of hundred meters over a fall distance of 4.4 km 
which is consistent with a 30-s fall time. It is also worth pointing out that this lateral 
displacement does not generate any acceleration on board because the displacement is 
due to an initial non-null velocity and not to external acceleration acting on the falling 
body. 

The diameter of the capsule is important for tolerating vertical gradients (wind shear) 
of the lateral wind without the need for a propulsion system to compensate for their 
effect. The balloon will move at the speed of the local wind once the floating altitude has 
been reached, i.e., the capsule will be at zero relative speed with respect to the local wind. 
If the wind vertical profile were constant, the capsule and the instrument package would 
move laterally during the fall with the same initial lateral velocity and hence maintain the 
same lateral distance with respect to one another. But, if the wind vertical profile 
changes, the capsule will experience a lateral force that will change its lateral speed while 
the instrument package will not experience such force. 
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The lateral displacement y of the capsule due to the wind shear V z = dV/dz over the 
vertical drop distance (Z = l/2gt 2 ) is as follows: 


30/3, g : 


( 8 ) 


where p is the atmospheric density, P, = m/(C D |S,) is the lateral ballistic coefficient of the 
capsule, S, the lateral area of the capsule, and g the Earth’s surface gravity. Equation (8) 
is simplified because the atmospheric density has been assumed constant over the drop. 
The equation, however, provides a good estimate of the lateral displacement of the 
capsule due to wind shear after adopting the average value p/ of the density along the 
drop. After calling S the frontal area of the capsule and assuming that C D1 « 10C D , S, « 
3S/(jt6), which are valid in approximation for a cylinder with aerodynamical ly-shaped 
nose and tail, we can relate the lateral ballistic coefficient to the frontal (low-speed) 
ballistic coefficient as follows: 



Consequently, eqn. (8) yields: 


. ±£? v ? 

8 ' 


(9) 


( 10 ) 



Figure 12 Lateral displacement due to wind shear vs. p o 
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For/V = 6xl0' 3 kg/m 1 (i.e., average atmospheric density between 40 km and 35 km of 
altitude), V z = 0.005 s' 1 , and 6 = 0.25, we obtain the results shown in Fig. 12. The 
accuracy of eqns. (8) and (10) was also checked through numerical integration of the 
equation of lateral motion after assuming an exponential air density profile. The 
displacement error is less than 5% by using the average density value. The value adopted 
for the wind shear of 0.005 s' 1 is equivalent to a vertical gradient of 10 knots per km. 
This value is twice as high as the maximum wind shear reported 1 " for the Air Force 
balloon base at Holloman, New Mexico. 

The lateral displacements due to wind shear are relatively small for rather 
conservative values of wind shear, free-fall times and ballistic coefficients greater than 
6000 kg/m 2 , which are easy to obtain. These results show that there is no need for a 
thruster system to compensate for the effect of the lateral wind acting on the capsule. 
The geometry and mass of the capsule can be chosen in a way to accommodate the 
presence of wind shear. Furthermore, if the balloon is launched during the periodically- 
occurring wind reversal times (in April-May and September-October) the vertical wind 
gradient is much smaller than the value adopted for the computations shown here. 

The capsule displacement due to wind shear has to be taken into account when computing 
the internal diameter of the capsule (where the instrument package falls). However, it 
will be shown later on that other factors (e.g., gravity gradients) are more important in 
determining the capsule internal diameter. 
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Science Concept Analysis 


The error analysis has been extended with respect to what is reported in Ref. xvi to 
include a thorough analysis of the gravity gradient generated by the distributed mass of 
the capsule, the concentrated masses on board the capsule, and the Earth’s mass for 
generic positions and orientations of the sensor. The acceleration noise inside the capsule 
has also been revisited after considering that in the new reference design (see later on) the 
chamber is fully cryogenic. The intrinsic noise components of the detector (Brownian 
and preamplifier noise) were also recomputed based on new information from the IFSI 
laboratory. As a results of these new analyses, requirements have been derived for: (a) 
the mass distribution of the capsule; (b) the tolerable mass and location of equipment on 
board the capsule; (c) the orientation of the sensor during free fall and the centering of the 
two sensing masses; and (d) the characteristics parameters of the detector that affects its 
intrinsic noise level. 


Acceleration Noise inside Capsule 

The experimental package moves at very low speed inside the capsule. Consequently, 
the residual gas inside the vacuum chamber produces a minute force on the free-falling 
package with a frequency content centered at f = 1/tf s' 1 where t f is the free-fall time. This 
gas thus affects the acceleration of the instrument package in a frequency range well 
removed from the signal frequency. The acceleration, as a function of pressure inside the 
chamber, is as follows; 


a 


D 


V 2 p 
2m RT 


(ID 


where A and m are the frontal area and mass of the instrument package, respectively, V 
is the maximum velocity of the instrument package with respect to the falling capsule, R 
is the gas constant, T the temperature of the residual gas, and p the pressure inside the 
chamber. Because of the new reference design (see later on) with a fully cryogenic 
vacuum chamber, the residual gas in the chamber is refrigerated Helium. After assuming 
C D = 2.2 (for a free-molecular regime), A = 0.1 m 2 , m = 30 kg, V = 0.5 m/s (obtained by 
integrating eqn. 3 once up to 30 s), R = 2078 J/(kg-K), T«5K, and p = 10' 6 mBar, eqn. 
(11) yields a D = 10 12 g. The spectrum of this acceleration is centered at a frequency 
0.033-0.05 Hz for free-fall times t f in the range 20-30 s. Consequently, the magnitude of 
the acceleration at the signal frequency f s , which is in the range 0.2-0.5 Hz, is well 
smaller than 10 12 g. This acceleration is a common-mode acceleration which is further 
reduced by the common-mode rejection factor (CMRF) of order 1 0' 4 . 

Furthermore, the vacuum strongly attenuates the propagation of perturbations from the 
walls of the capsule to the free-falling instrument package. The estimate of the 
acceleration at the falling instrument package produced by the vibrating walls of the 
capsule are based on the experimental data measured on board the system Mikroba*™. 
This system shares the fall from a stratospheric altitude; it is not, however, a free-falling 
experimental package inside the shielding capsule. In Mikroba, the measurement 
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package or experiment is solidly attached to the walls of the falling capsule. Moreover, 
Mikroba is not propelled during the first 30 s (like our experiment) although it is 
propelled downwards during the next 30 s. Once the magnitude of the acceleration at the 
walls la wall l is known, the magnitude of the acceleration at the falling package lal can be 
readily computed as explained in the following. The motion of the vibrating walls 
increase the kinetic energy of the gas molecules above the thermal velocity. The kinetic 
energy variation is then expressed as a pressure variation Ap of the gas after equating the 
increase in kinetic energy to the work done on the gas molecules by the vibrating walls. 
We then assume, conservatively, that the pressure perturbation Ap acts on one side only 
of the instrument package in order to compute an upper bound of the acceleration 
disturbance imparted to the package. The upper bound of the acceleration at the 
instrument package is as follows: 


lal - -^l„ 
mv 


wall 


Ap 

nuav 


I a 


( 12 ) 
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Figure 13 Acceleration measured on board the Mikroba capsule (with accelerometer 
solidly attached) during the fall [Kretzchmar, 1999] 

In equation (12), A and m are the cross section and mass of the instrument package, v 
is the thermal velocity of the residual gas, a and a n . a „ are the accelerations of the package 
and the wall, respectively,/? is the pressure inside the capsule, u wall = a wall /oj is the velocity 
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and to the angular frequency of the vibrating walls. If, for example, p = 10" mBar and the 
temperature of the gas inside the capsule is slightly higher than LHe, we obtain an 
attenuation factor a/a wall = 6xlO' 10 at the signal frequency of our experiment (~0.5 Hz). In 
other words, the high vacuum provides an excellent attenuation of the wall vibrations. 
Moreover, the capsule walls do not necessarily vibrate at the signal frequency of the 
instrument. The cryostat will be designed with structural frequencies much higher than 
the experiment signal frequency. However, in this early stage of the design and for 
conservativeness, we assume that the wall acceleration has a component at the signal 
frequency with an intensity equal to the largest magnitude of the acceleration recorded on 
Mikroba during the first 30-s of fall, that is, a wal , < 10 4 g (although not very visible in Fig. 
13). We then obtain an acceleration at the instrument package of order 10' 13 g under 
rather conservative assumptions. Consequently, the free-falling capsule reduces the 
acceleration noise to values unmatched by any other Earth-based drop facility and 
comparable to values achieved on board the Triad drag-free satellite xv "'. 

The acceleration components above are common-mode-type (i.e., they affect equally both 
sensing masses) thus they can be further reduced by the common-mode rejection fact of 
the differential accelerometer. With a typical value of 10‘ 4 for the CMRF, the influence 
of these accelerations on the differential measurement is made negligible. 

The acceleration noise components produced by the residual gas in the capsule are 
proportional to the pressure inside the capsule. The pressure can be reduced in 
successive flights if, for any unanticipated reasons, its influence on the measurement 
proves to be greater than expected. It is, in fact, well within the state of the art to obtain 
pressures at room temperature as low as 10' 8 mBar in large volumes. 

Internal noise of detector 


The most important internal noise sources for a high accuracy mechanical detector like 
the one proposed for this WEP test are: (1) preamplifier noise; and (2) thermal noise 
(Braginsky, 1974 x,x ; Giffard. 1976 xx ). The combined effect of these two noise 
components upon the acceleration spectral density S a of the detector's output is given by 
the following equation for an instrument with the measurement frequency smaller than 
the resonant frequency coo, a measurement bridge pumped at the frequency co P (of 
typically tens of kHz) and a preamplifier that matches the transducer impedance xxl : 


= 


4Wl + r 2a) 0 ' 

m eff [Q n Xoj P/ 


-. 1/2 




(13) 


In equation (13), the two terms in round parentheses correspond to the Brownian noise 
and the preamplifier noise, respectively; coo is the detector resonant frequency; k B the 
Boltzmann's constant; T the ambient temperature; T n the preamplifier noise temperature; 
Q the quality factor; m e ff the effective mass of the sensing mass; and X the 
electromechanical transducer factor. The effective mass is used to convert a rotation of 
the sensing mass into a translation of equal energy. Its relationship to the mass m is: m e ff 
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= m(r\ { /l) 2 where r| t and l are, respectively, the radius of inertia and the lever arm with 
respect to the rotation axis of the sensing mass. With the geometry of the sensing masses 
in our detector m e ff * 1 .8 m. 

Clearly, from eqn (13), we see that the sensitivity of the detector increases by 
decreasing the resonant frequency and the temperature, and by increasing the mass of the 
sensing masses and the Q-factor. Liquid He (4.2 K) refrigeration will be used to provide 
low Brownian noise and a high Q-factor. These are necessary conditions to achieve the 
desired measurement accuracy. In order to derive requirements for the detector, we 
assume that the contribution of the Brownian noise is about equal to the contribution of 
the preamplifier noise. In this case, if we set our experiment accuracy goal to several 
parts in 10 15 (with 95% confidence level and a 20-s integration time) each one of the two 
noise components should be about 10" 14 g/^fHz or smaller. These noise requirements 
imply the following (see also the section on Requirements Development): 

T < 10 K; T„ < 60 mK; m > 2 kg; to^Q < 6n/ 10 5 rad/s (14) 

These requirements do not exceed the state-of the-art but they do require a very 
careful construction of the detector with low dissipation and the use of very-low-noise 
preamplifier. Key quantities like the Q-factor at low temperature and the preamplifier 
noise will be measured experimentally by our partners at IFSI once a prototype 
laboratory detector is built and operated at LHe temperature. Noise contributions other 
than the intrinsic noise components of the detector should be kept at a lower level in 
order to make them smaller than the intrinsic noise. 

Gravity Gradients 

Capsule Gravity Gradients 

The gravity gradients generated by the distributed mass of the chamber/cryostat and 
their effects on the differential measurement are analyzed in the following for a generic 
position of the detector inside the capsule and a generic orientation of its spin axis with 
respect to the gradient field. 

Gravity gradient for a mass distribution with cylindrical symmetry 

For a mass distributed with cylindrical symmetry, the resultant gravitational 
acceleration has two components: 

a =■ acceleration component along the cylinder axis 

a,= acceleration component along the cylinder radius 
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Figure 14 Reference frames for gravity gradient analysis 

After erecting a Cartesian reference system xyz as depicted in Figure 14, the 
components of the gravity gradient tensor are computed according to the following 
transformation formulas: 

dx = cos(0) ■ dr - r ■ sin(d)- dd 

dy - sin(d) • dr + r ■ cos(6) • dd (15) 

da x = cos(0) • da r - a r • sin(d) • dd 
da v = sin{d)-da r +a r - cos (d)-dd 

Setting 0 = 0 and indicating the spatial derivative with a second subscript: 


a y = 0 

and 


a yy ~ a r I r ( 16 ) 

a xz = a zx ~ a rz 

a r- =«,v =° 

Hence the gravity gradient tensor has the form: 

0 a rz ' 

a,./r 0 (17) 

0 a,. 
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As a result of the Laplace equation, the trace of the gravity gradient tensor is equal to 
zero, that is 


a rr + a r lr + a zz =0 

In the singular case of r = 0 the limit calculation yields: 


(18) 


^(r=0> ~ 


a rr 0 0 

0 a„ 0 

0 0-2 a. 


(19) 


Gravity gradient matrix of a rotating body 

In general a gravity gradient matrix has the form: 


r = 


The rotated matrix F after a 0 = cot rotation is: 
T©= R^ritf 

where R a is the rotation matrix and R H T its transpose. 


r *v 


r v. 

r vz 


r a 


( 20 ) 


( 21 ) 


After a rotation about an axis (i.e., the x axis), the rotated matrix has four components 
modulated at co, four components modulated at 2co and one component that is not 
modulated. 


The to-modulated components of the transformed matrix are: 

r; 2 =r;,= r*. sm(cy • t) + cos(a> • t) 

r ,' 3 = r 31 = r t . cos(cu • t) - r sin(a> • t) 


(22.1) 

( 22 . 2 ) 


where the rotated axes are labelled x’ = 1, y’ = 2, z’ = 3. In summary, the off-diagonal 
components F )2 = T’ 2 i and F l3 = T’,, of the gravity gradient matrix produce components 
that are modulated at the rotation frequency. 
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Gravity gradient matrix projected onto body axes 

In a general case the body reference frame placed at the CM of a sensing mass of the 
detector can be identified with respect to the cryostat frame by means of 3 successive 
rotations as follows: 

1 - Rotation +a around z axis (azimuth rotation) 

2 - Rotation +(3 around y’ axis (elevation rotation) 

3 - Rotation +oat around x” axis (spin rotation) 

In the computation of gravity gradients, these rotations can either be rotations of the 
sensing mass with respect to the cryostat or, equivalently, rotations of the cryostat with 
respect to the sensing mass. In the former case, and solely for the reason of pointing out 
a typical geometrical situation, the first and second rotations could, for example, be 
caused by the detector dynamics during free-fall (e.g., precession of its body axes) while 
the third rotation is the cot rotation of the detector about its longitudinal axis aimed at 
modulating the signal. 

Clearly, we are mostly concerned about the components of the gravity gradient matrix 
that contain a frequency to equal to the modulation frequency of the signal. We can 
choose the body axis y’ = 2 to coincide with the sensitive axis of the accelerometer and, 
consequently, we are only concerned with the component r' 2l of eqn. (22.1). In general, 
the moduli of the two components r’ 12 and r’, 3 are the same and they can be written as 


follows: 

< 23 ) 

After rotating the original matrix by two rotations a and P (where a is the azimuth of the 
spin axis with respect to the radial and (3 is the elevation with respect to the capsule 
equatorial plane) the expressions of T xy and T xz in eqn. (23) are as follows 

r vv = ^!sin(/5)sin(a)-A:2 cos (/3)sin(2a) (24.1) 

r A . = kj sin(2/I) + k t cos(a)cos(2)S) + k 2 sin(2/3)cos 2 (a) (24.2) 

k l = a x . (24.3) 

k 2 = —(a xx - a yy ) (24.4) 

= ^(«vv - «-) (24,5) 
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where the a u are the matrix components before the rotations are carried out. In the case of 
a body with cylindrical symmetry and for cylindrical coordinates, eqn. (16) yields a xx = 
a„, ay V = a/r, a^ = a u , a X2 = a^ and the other components are null. 

After carrying out a numerical analysis of eqns. (24) and taking into account that 
inside a cylinder k 2 is always at least one order of magnitude less than k, and k 3 , we find 
that the maximum value for x occurs for a = 0 (where a is the azimuth). This result 
implies that the maximum disturbance of the capsule gravity field on the differential 
accelerometer is produced when the capsule moves radially with respect to the sensor 
(see Figure 15) in such a way that the spin axis is oriented along the radius of the cylinder 
through the sensor and the capsule has been displaced radially with respect to the sensor 
(e.g., by wind shear). 



Figure 15 Geometry of sensor and capsule (viewed from the top) for strongest gravity 
gradient affecting the measurement 

On the opposite end, if the motion of the capsule is such as to keep a close to 90°, 
that is the spin axis is orthogonal to the radial, the disturbance is minimum. In any case 
since the translational motion of the capsule is not predictable nor controlled, the worst 
condition is analyzed setting a equal to zero and varying the angle p. After doing so 
eqns. (23) and (24) yield: 

X = (k 3 + & 2 )sin(2/3) + k 1 cos(2/J) (25) 

Equation (25) summarizes the disturbances induced by the cryostat mass modulated 
at the measurement frequency. This equation is important for the cryostat/capsule design. 
The variations of the quantities k,, k 2 , k 3 inside the cryostat for different shapes and sizes 
are analyzed numerically in the following subsection. 

Variation of k,, k 2 , k 3 in the cryostat/capsule 

A numerical code has been developed in Matlab to compute the gravity gradient 
matrix inside a distributed, massive cylinder. The program, which uses a very large (of 
order 10 4 ) number of mass points, can map the desired components of the gravity 
gradient matrix inside the enclosed surface. The code has been applied to a number of 
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cylinders with different H/D ratios (where D is the diameter and H the height) and 
different shapes of caps. Each run takes a couple of hours on a Pentium III PC. 

Cylindrical cryostat without caps 

We first analyzed the gradient field for various cases of cylinders with various H/D 
(height over diameter) ratios to conclude that the gradient field is strongly reduced (for 
the components of interest) for H/D > 1 

The following results are for a cylinder of uniform mass distribution with the following 
characteristics: 


Cylinder mass = 500 kg; Dimensions: 1 m (dia) x 1 m (height) 



Figure 16 Schematic of cylinder and reference frame 
Vertical profiles at r = 0 (along centerline) 



(continued) 
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Vertical profiles at 
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Figure 18 All components of gravity gradient along a vertical profile at r = 10 cm (units 
are f„ = ma r2 in kg/s 2 = s‘ 2 for 1 -kg test mass) 


Vertical profiles at r = 20 cm 



(continued) 
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The following is a contour plot of the k, component inside the cylinder and the radial 
profiles of this component along radii at various distances from the cylinder’s equatorial 
plane. 
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Figure 22 Contour plot of capsule gravity gradient component k, = a T/ (s' 2 ) for a cylinder 
with H/D = l-m/l-m 
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z = 10 cm (above cylinder’s equatorial plane) 
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(continued) 
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Figure 23 Gravity gradient (a r , = k, component) for latitudinal sections at different 
distances above the cylinder’s equator (units are f r/ = ma„ in kg/s 2 = s' 2 for 1-kg test mass) 





k, (1/s 2 ) 



Figure 24 Contour plot of capsule gravity gradient component a re (s' 2 ) for a cylinder with 
H/D = 1.5-m/l-m 


Figure 24 shows the contour plot of the k, = a„ component for a cylinder with H/D = 
1 .5-m/l-m. Note that the strength of the gravity gradient in the area of interest (near the 
cylinder's centerline) is strongly decreased thanks to the lengthening of the cylinder in 
the vertical direction. 
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As the cylinder is stretched along the vertical, the effect of the cylinder edges (which 
produce stronger gravity gradients) is smaller the closer the detector is to the cylinder’s 
centerline, where the detector free falls. In the following plots we will show all the 
gravity gradient terms k, , k 2 and k 3 (modulated at oa) inside a cylinder with dimensions 
close to those that we are considering at the present stage of development of the design. 
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(continued) 
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Variation of lq. k 2 . and k 2 along free-fall trajectories 


Based on the evaluation of the dynamics of the instrument package relative to the 
capsule (see also later on), we can map trajectories of the instrument package inside the 
capsule and evaluate the maximum values of the gravity gradients that the detector will 
experience during the fall. 

It is important to evaluate the strength of the gravity gradient field inside the cryostat 
along a worst-case trajectory of the instrument package (sensor) that moves with respect 
to the cryostat/capsule during free fall. Based on worst-case wind shear conditions the 
trajectory (in z-r coordinates) of the sensor with respect to the cryostat can be expressed 
as follows: 


z = z 0 


-f t h exp{ct 2 ) 
Ho 


r ~ r o ~ 


dg 2 f + 

8*A> 


(z-Zo)siny 


(26) 


with a = 0.149636, b = 0.001692, c = 3.084, and d = 6xl0" 7 (see section on Optimization). 
In equation (26), p o is the low-speed ballistic coefficient, z 0 and r 0 are the coordinates of 
the point of release in the cryostat coordinate frame and y is the angle of the capsule’s 
longitudinal axis with respect to the local vertical which (at this stage of the project) 
represents a reasonable upper limit for the verticality error of the capsule during the fall. 
This trajectory will beused in the next subsections. 
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Cylindrical cryostat with flat caps 


D 

The Matlab routine has been modified to include 
cryostat caps of different shapes. The routine 
creates a two-dimensional mesh of point masses 
uniformly distributed on the average surface of 
the cylinder and its caps. In the case of flat caps 
the mass distribution results in a closed 
cylindrical surface of height H and diameter D. 
The gravity gradient field has been mapped on 
the z-x plane where x coincides with the 
cylinder’s radial and z with the longitudinal axis, 
as far as s = 10 cm from the top and bottom and c 
= 20 cm from the side walls. 

The point of release P lies on the symmetry axis of the cylinder and at d = 40 cm from the 
top. The sensor trajectory obeys eqns. (26) with p o = 10000 kg/m 2 and the capsule 
verticality error y = 5° has been conservatively assumed to produce a lateral displacement 
of the sensor in the same direction of the wind shear. The other parameters in eqns. (26) 
also represent a worst-case scenario for lateral displacements. The key quantities k,, k 2 
and k 3 are plotted as contour plots on the x-z plane in Figs. 26. 
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The quantities k,, k 2 , k, are then computed along the trajectory of the instrument 
package and the power spectral densities are computed over the free-fall time of 26 s. 
Results are shown in Fig. 27 where the peaks with frequency l/t f = 1/26 Hz due to the 
free-fall duration, are clearly visible in the spectra. 



Msi 


Figure 27 k, Jt 2 Jc_, along the trajectory and their spectra. 
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Figure 26 shows that the quantity k 3 is relatively larger than k, and k 2 . Consequently, 
eqn. (25) poses a limit for the angle (3 that defines the capsule attitude with respect to the 
spin axis of the sensor. After neglecting the much smaller k 2 (and k, which depends on 
cos(2|3)) we find that the maximum allowable p is: 

<27) 

With a Xmax of 10* 9 s" 2 and the results shown in Fig. (27), we obtain a limit of 4.8 deg for 
P, which is a relatively large value. 

Cylindrical cryostat with hemispherical caps 



A cylindrical cryostat with hemispherical caps has also 
been analyzed. The gravity gradient distribution has been 
mapped on the x-z plane, where x is the cylinder’s radial 
and z the longitudinal axis, between the base of the upper 
and lower hemispheres and as close as c = 20 cm from 
the side wall. 

The release point P (and starting point for the simulation) 
lies on the symmetry axis of the cylinder and at the base 
of the upper hemisphere. The sensor trajectory is the 
same as in the previous case. 
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k3 (1/s 1 ) 



Figure 29 k,, k 2 , k 3 along the trajectory and their spectra. 

The analysis carried out here and additional results not shown in this report lead to 
the definition of a stand-off distance of about 40 cm between the sensor CM and the 
heavy part of the cryostat walls to provide w-modulated components of the gravity 
gradient that are sufficiently low for the sizes and masses relevant to this project. 
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We compute here the Earth's gravity gradient tensor and we analyze the effects of 
Earth’s gravity gradient components on a rotating detector with a generic orientation of 
its spin axis with respect to the gravity gradient field. Let us consider the gravitational 
potential per unit mass at a point (x, y, z) with respect to the detector’s center of mass: 



where R x , R Y , R z are the components of the radius vector R from the Earth’s center to the 
detector’s CM (in which Z is the local vertical) and p is the Earth’s gravitational 
constant. After projecting about the detector’s body axes in which x is the spin axis and 
calling 0 = tot the rotation about the spin axis and <]> the elevation angle of the spin axis 
with respect to the horizontal plane: 


R, = R(t)s'm(tp) 

R y « /?(/)cos(0)sin(atf) (29) 

R. = R(t)cos(<p)cos(oM) 


The gravitational acceleration in body axes is obtained by substituting eqns (29) into 
eqn. (28) and computing the gradient: 



(30) 


The components of the gravity gradient matrix in body axes are finally computed by 
taking a further derivative with respect to the spatial coordinates, to yield 


8xx ■ -- &[-2 + 3cos : (0)] 

K 

(31.1) 

8„ = 3-^ r sin(ctf)cos(0)sin(0) 

A 

(31.2) 

g xz = 3-^cos(cur)cos(0)sin(0) 

R 

(31.3) 

gyy = — £j[l - 3cos 2 (0) + 3cos 2 (fttf)cos : (0)] 

A 

(31.4) 

gy. = 3-^-cos(a>f)sin(ft)t)cos 2 (^) 
R 

(31.5) 
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(31.6) 


= -^-[-1 + 3cos 2 (cot)cos 2 (0)] 

Here again, the components modulated at to are g xy and g X2 (in which x is the spin 
axis). In other words, if the spin axis lies on the horizontal plane, the detector only sees 
components modulated at 2o) but if it is not, components modulated at the frequency 
o) appear. The strengths of these components are proportional to the tilt angle with 
respect to the local horizon. Note also that the effect of the Earth’s gravity gradient on a 
rotating body can be readily applied to the space-based tests of the Equivalence Principle 
in which the only difference from the balloon-based experiment is the slightly larger 
value of the radial distance from the space-based sensor to the Earth’s center. 

An alternative way of portraying the origin of the co-modulated components of the 
Earth’s gravity gradient field is by considering the following. If the z-axis of the body 
reference frame is directed along the local vertical (that is the spin axis x of the sensor is 
on the horizontal plane) then the gravity gradient tensor is 


r £ 

1 


GM < 


R 3 


-1 

0 

0 


0 0 

-1 0 

0 2 


(32) 


The tensor does not change under a rotation about the z-body axis due to its structure 
which reflects a symmetry about the radial line. Consequently, we can choose the 
azimuth orientation of the spin axis at will (let us call it the y-body axis). A rotation § 
about an axis perpendicular to the spin axis produces terms g’ X2 = g’ 2X in the transformed 
tensor: 


r/ = R,v l R', 


GMe 

R 3 


g<& 

0 

gQ 

0 

-1 

0 


0 



where 


g<§. = -l+ 3sin 2 <p 
g© = -l+ 3cos 2 0 
= g zx = -3sin0cos0 


(33) 


(34) 


The g’ X2 and g’ 2X terms are subsequently modulated at the frequency <o by the rotation 
0 = cot about the spin axis as shown previously. 

Since there are terms modulated at the signal frequency to, we have to make sure that 
they are kept lower than the accuracy with which we want to measure the signal. From 
the detector point of view, there will be requirements imposed on the centering of the 
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sensing masses and their orientation with respect to the Earth’s gravity field as shown 
later on. 


Disturbances induced bv concentrated masses on board the capsule 

Let us consider the reference system (x, y, z) fixed to the sensor with origin at the 
center of mass and with x oriented along the spin axis and let us indicate the position of a 
point mass nip in proximity of the sensor in spherical coordinates (6,y,p) (see Fig. 30). 



Figure 30 Geometry of the sensing mass and reference frame 
The position of nip is then expressed as: 
x p = pc os<5 cosy 

y,, = p sin 6 cosy (35) 

z p = Psiny 


The gravity gradient matrix at the detector due to the gravity field induced by a mass 
point m p is: 


G-m p 

— 


r-. 2 2 

3jc - p 

3yx 

3zx 


3xy 3xz 

3x 2 - p 2 3 yz 

3yz 3x 2 - p 2 


(36) 


Considering a sensor that rotates with respect to a fixed point mass in its proximity 
we obtain the two w-modulatcd components already shown in the previous paragraphs: 

r i 2 = H, = r i3 sin(ft> • t) + r i2 cos(cu t) (37. 1 ) 

r,3 = r;, = r l3 sin(cu t) - r 12 cos(a> • t) (37.2) 
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The moduli of the two oo-modulated gravity gradient components T’ 12 and r ’ 13 are the 
same and can be expressed as follows: 


X = 


P 5 


MV 2 ' +J r 


(38) 


which shows that the masses located on the plane y-z (i.e., x = 0) do not generate 
disturbances with the same frequency as the measured signal. Substituting eqn. (35) into 
(38) and extracting p yields the minimum distance for a point mass to produce a 
disturbing gradient equal to or less than the critical gradient a gg . max : 


Pmin = 


G' m, 


l a 


-|cos<5cos y\^ sin 2 <5cos 2 y + sin 2 y 


gg- max 


1/3 


(39) 


Setting a limit of 10' 9 s' 2 for a max we plot the locus f(p,y,6) = 0 of the points in space 
with a gg = a gg _ max in Figure 31 for a disturbing point mass of 1 kg. Next, meridian sections 
(rotated about the z-axis by the meridian angle 6) of the same locus are plotted for 
different values of the angle 6 in Fig. 32 where r is the radial direction. 


z(m) 



Figure 31 Locus of o-modulated gravity gradient component with strength = 10‘ 9 s' 2 
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Figure 32 Meridian sections of locus in Fig. 31 for different angles 6 [6 = 0 (i.e., y-z 
plane) - solid black line; 6=15°- blue dots; 6 = 30° - red dash; 6 = 45° - gray dash dot] 

The worst case meridian section at 6 = 45° is plotted for different values of the 
perturbing mass in Fig. 33. 



Figure 33 Meridian sections for 6 = 45° and different values of perturbing mass nip 
(m p = 1 kg - solid black line; m,, = 10 kg - blue dots; mp= 100 kg - red dash) 

The previous analysis defines exclusion zones for concentrated masses on board the 
capsule. In general, masses can be placed rather freely on the y-z plane (perpendicular to 
the sensor spin axis). Masses lying on this plane generate only 2co-modulated components 
whose strength only needs to be reasonably smaller than the upper bound of the dynamic 
range of the sensor. The Earth itself produces such 2o)-modulated components with a 
strength equal to 3X10 -6 s 2 that is well stronger than the sensor sensitivity (for realistic 
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values of the sensing mass CMs errors), but is about 3-orders of magnitude lower than the 
upper bound of the sensor dynamic range. For the 2to-modulated term, the equivalent 
Earth is a mass of 22,500 kg at 1 -m or a 22.5 kg at 10 cm from the sensor. 

Requirements related to the to-modulated components are more stringent and, 
consequently, we will concentrate on these components which have been dealt with in 
this analysis. In summary, concentrated masses should be placed as close as possible to 
the y-z plane (perpendicular to the sensor spin axis). For masses away from the y-z 
plane, Fig. 33 defines the exclusion zones from the sensor for different mass values under 
the worst possible condition of masses placed on the 45° meridian plane. 


Effect of Gravity Gradients on Differential Acceleration Measurement 

The differential accelerometer consists (from the mechanical point of view) of two 
sensing masses with ideally coincident centers of mass (CM). The equivalence violation 
signal is measured as a differential displacement along the y-body axis of the sensor 
which is orthogonal to the spin axis along the x-body axis. In reality the two centers of 
mass (or more appropriately centers of gravity) do not coincide and CM 2 (i.e., the CM of 
mass 2) is displaced by a position error vector 6 with respect to CM, as follows 


6 = 


< 5 , 

<5. 


(40) 


We can place the body reference frame at CM, and compute the differential 
acceleration due to gravity gradients by simply multiplying the gravity gradient matrix in 
body axis , that is 


da f 



da y 

= (r £ + r c + r M ) 


daz 


d z 


(41) 


where T E , T c , r M are the gravity gradient matrices of the Earth, the distributed capsule 
mass and concentrated masses on board the capsule, respectively. Since the differential 
accelerometer measures only the component along the y-body axis, we obtain finally: 

&v = at + 1$ . + r" )<5, + (rj + r£ + r" )<5, + (r* + r£ + r" )S : (42) 

in which r E , T^ , and T^ are the components modulated at the signal frequency to while 
the other terms in eqn. (42) are modulated at 2to. In conclusion, the disturbing 
differential acceleration along y produced by gravity gradients can be expressed as 
follows: 
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da y < 3-^jCOs(<p)sin(<p)sin((ot)6 x + [(& 2 + & 3 )sin(2/3) + fc, cos(2/S)]sin(a>?)(5 t 

R ' (43) 

+ a M {(o) + f y (2(o,8 y ) + f.(2wA) 

where k,, k,, and k 3 are the gravity gradient terms (see previous subsections) generated by 
the capsule in the capsule-body reference frame, <(> is the elevation of the spin axis with 
respect to the local horizon, 0 is the elevation of the spin axis with respect to the 
equatorial plane of the capsule, a M (u)) is the oo-modulated disturbing acceleration (in 
functional form) produced by concentrated mass on board the capsule (see previous 
subsection) and f y (2to, 6 y ) and f z (2co, 6 Z ) represent all the other 2to-modulated components 
which have been separated in eqn. (43), according to the centering error components. 
Note that the 2w-components depend on the centering errors 6 y and 6 Z while the o>- 
components depend only on the centering error 6 X . The less-than sign in eqn. (43) is due 
to the fact that, on the right hand side of the equation, we have adopted the strongest 
value of the to-modulated gravity gradient component of the capsule, that is, for a = 0 
(see Fig. 15). Moreover, from the analysis of the capsule gravity gradients, we have 
concluded that if we keep the sensor (at the CM) about 40 cm away from the heavy part 
of the chamber/cryostat walls, the to-modulated gravity gradients are well below the 
critical value of about 10' 9 s' 2 . Based on similar reasoning, we assume that the 
concentrated masses on board the capsule are placed outside of the exclusion zones 
(defined in the previous subsection) in order to keep them below the critical value. In 
other words, an appropriate design and a careful mechanical construction of the sensing 
masses (6 X of order microns) will make the gravity gradient contribution of the capsule 
and the concentrated masses on board the capsule negligible. 

To attenuate the effect of the gravity gradient of the Earth we have to make sure that 
the product sin(^)cos(^»)^is sufficiently small. In other words we can trade the position 
error between the CMs of the sensing masses along the spin axis 6 X for the tolerable angle 
<|) of the spin axis with respect to the local horizontal. For small values of <J>, we readily 
compute that for the first term on the right hand side of eqn. (34) to be smaller than, let 
us say, 10' 15 g, the product <j)6 x must be smaller than 0.1 deg-pm. This requirement must 
be considered in the design of the detector, the release mechanism, and the capsule 
leveling system of the capsule. The complexity of some subsystems can be traded for the 
simplicity of other subsystems among those three devices. 
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Theoretical Impacts 


Contribution of Parity Non-conserving Weak Interactions 

Introduction 


A theoretical question that we sought to address can be simply stated: What is the 
contribution to the mass-energy of each material due to the parity non conserving part of 
the weak interaction? The reason that we sought an answer to this question was to find 
out if our experiment would be sensitive enough to determine whether or not this 
contribution obeys the weak principle of equivalence. 

Unfortunately, the present level of development of physics does not allow us to 
address our question reliably. Knowledge of the physics of nuclear matter is too 
primitive. We have therefore had to make a number of "reasonable," but nonetheless 
somewhat arbitrary, assumptions to carry out the calculations. The discussion below 
mentions each of these assumptions. The results of the calculations indicate that our 
experiment will not be sensitive to the contribution to mass-energy of the parity non 
conserving part of the weak interaction. 

Evaluation of contribution 


The materials to be compared in the experiment should have binding energies stored 
in forms which are as different as possible. For example, if gravity couples differently to 
protons and neutrons, we should compare elements with different proton to neutron 
ratios. A new long range force could also be detected by comparing such elements. A 
force coupling to baryon number would cause an acceleration proportional to the total 
number of protons and neutrons divided by the mass, or for a single nucleus, to 

(m„ - m n ) 

(Z+N)/(m Z + m„ N)~l/m N (l x). Here x is the ratio Z/A and Z, N and A are 

m N 

the number of protons, neutrons, and combined nucleons. The mass m N is used for the 
common mass of protons and neutrons. In order to observe this effect, it is best to 
compare heavy elements with x~l/3 to light elements with x~l/2. 

As we will see below, the energy of a nucleus can depend very sensitively on the 
wave functions of the protons and neutrons. We will discuss a force coupling to the 
product of proton and neutron densities. Such a force is more significant for nuclei in 
which proton and neutron wave functions have greater overlaps. This suggests choosing 
a nucleus with a magic number to compare to a less stable nucleus. The filled shell 
structure of a nucleus with a magic number may imply greater overlap of the wave 
functions of the protons and neutrons. 
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Suppose gravity couples differently to the energy produced by the parity violating 
part of the weak force. We must calculate what fraction of a nucleus’s energy is stored in 
this form, although we will see that it is not very well defined. First of all, right-handed 
neutrinos do not appear to exist, so it is not just the weak force which violates parity. The 
model used in Ref x ’ u ' to calculate the total weak force energy of a nucleus is a current- 
current interaction. To discuss this, let us consider the Lagrangian for the interaction of 
the weak force with quarks and leptons. Schematically, it is A ft (J 5 " + J") , where J and 

J 5 are the weak vector and axial currents, and A is the field of the weak force. (This 
expression must be summed over the three types of W particles; this formula ignores the 
fact that the observed W and Z particles are actually linear combinations of these three 
fields with a fourth field.) By including both terms, we ensure that the coupling to right- 
handed particles cancels. The first term here is the parity violating term. However, if we 
use a contact interaction, taking into account the large mass of the W particles, we find an 
interaction of the form L e(r = (J u + J* ){J + J Sfl ) . Multiplying this out, we get four 

terms — only the terms coupling the axial current to the vector current seem to violate 
parity in this description of the weak force. Another difficulty is that the separate vector 
and axial currents are not actually well defined if right-handed neutrinos do not exist. 

Since the parity conserving parts of L cff contribute one part in 10 8 to the mass of a 
nucleus, current limits imply that gravity’s coupling to these terms must differ by less 
than 10 4 from its coupling to other matter. The parity violating terms of L eff contribute 
much less to the mass of a nucleus. If li> is the state of the nucleus with weak forces 
neglect, then the first order perturbation of its energy due to these terms is <il •/"./* Ii>, 

which vanishes by symmetry. (Ii> is a parity eigenstate, so the operator changes its parity, 
and gives a state orthogonal to <il.) This ensures that the contribution of the parity 
violating terms is second order in perturbation theory, and therefore small enough to be 
unconstrained by previous experiments. However, L t . r , is only useful for calculations. It 
would be strange for gravity to couple differently to the parity violating part of the weak 
force unless it coupled differently to in the original Lagrangian. In this case 

^ would also couple differently to gravity, and this would be a much larger effect. 

While we are deciding which terms of the effective Lagrangian will be considered 
to be parity violating, we should notice that there are other parity violating interactions 
involving pion exchange™" 1 ™ 4 . Ref. xxiv probably overestimated the contribution of 
J iJ* 1 to the weak energy of a nucleus, because of the assumption that the nucleons are 

distributed independently of one another in a nucleus. In fact, protons and neutrons are 
generally not closer to one another than their radius of 10 ,s m, or (160 MeV) 1 as 
compared to the range 1 GeV 1 of the weak force. Thus, the weak interaction of nucleons 
occurs only very rarely. The pion exchange force has a range larger than the size of a 
nucleon, so the assumption that nucleons are independently distributed is more accurate 
when considering this force. 

The interaction between nucleons and pions is described by the Lagrangian 
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An. = -‘SnNN Ny 5 f ■ jfN - N(f xH) 2 N 


(Ref xxv) where we have used isospin notation, so that N is an isospinor of proton and 


neutron fields: N= 


r is a vector of three Pauli matrices, and jt is the vector of pion 


fields, which are usually written as (- 


n +n 


IT' 1 V2 


. Jt JT 




The first term in the 


Lagrangian represents pion exchange due to the strong force (£*^=13.45), the second to 
the weak force {f R =4.54xl(T 7 ). This term violates parity. We will calculate 
contributions to the energy of a nucleus arising from this term. 


Any interaction in which a light particle like the pion is exchanged can be 
approximated by interaction potentials. For example, we begin with a simplified model 
of spin zero protons and neutrons which can exchange a neutral spin 0 pion (of mass mj. 
Let the interaction Hamiltonian be g(n + n+p + p)Jt. Then the amplitude for pion emission is 
-ig. Let us calculate the potential acting between protons and neutrons. The Feynman 
diagram is: 



The amplitude associated with this is A = y(-ig) 2 — — — — 7^ — • We have 

r-m„~ q-+m / 

omitted the time-like component of q, since this is the energy transferred between the 
nucleons, and is small compared to the momentum in the non-relativistic limit. 

Compare this to the calculation of scattering under the influence of a potential V. 

i ( p\ -r„ +p 2 -r p ) H h 'K + P 4 r ‘ P ) 

Let the initial and final states be It//, >= — and |t// f >= — - — — — . In time 

(2*r) 3 f (2 Jt) 3 

independent perturbation theory, we solve for the coefficient 
c f =< / I T(-cc,t) |/>* j< f\V\i> e' (E/ - E ' )l dt' > V fl 2nd(E f - E,) . 

-oc 

It is not hard to check that y (p x + p 2 - p 3 - p 4 )fd 3 xV(x)e'^ ■ 
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The Feynman amplitude is closely related to this coefficient; however, it is 
calculated with the wave functions normalized so that 

< PiP 4 I Pi Pi >= 4£, E 2 (2tt ) 3 6 (p ] - p } )6 (p 2 -p 4 )~ 4m N \2ji ) 6 d (p, -p 3 )d(p 2 - p 4 ). 

Thus we must multiply our value for c f by 4m N 2 {2nf . The relationship between c f and 
A is thus 

iA(2ji)*d(E 3 +E 4 -E ] -E.Jdip^ +p 2 -p 2 - p 4 ) = 4m N 2 (Inf c f . 

Therefore 

A = -4m N ijV (x)e u<x d i x , and the inverse Fourier transform of this gives V. We find 

that V (x) = ; — , where the transform is evaluated using spherical coordinates. 

{4m N -)4n\x\ 

This is the Yukawa potential, but it is not actually the true interaction potential of the 
protons and neutrons, since protons and neutrons really have spin 1/2, and have axial 
couplings to the pions. 

Expanding the Hamiltonian given above, we find the parity conserving and parity 
violating pion-nucleon interactions 

h pc =i 8^mf d3 x(py 5 p-ny 5 n)^ 0 + ^fl(pY 5 nji* +ny 5 pjt ') 
and 

H py = -jrfd } xpn7i' > — TTpji~ . 

The jfs are pseudoscalar particles and py^p and similar terms are scalars, and so it is the 

first of these interactions which is parity conserving. This is actually a convention, since 
if the second interaction had been discovered first, the pions would have been called 
scalar particles. A consequence is that an interaction between nucleons by exchange of 
pions is equivalent to a parity conserving potential as long as both pion-nucleon vertices 
are governed by the same interaction. If the pion is emitted according to the parity 
conserving interaction and absorbed according to the parity violating interaction, then the 
equivalent potential distinguishes between left and right. We call the three potentials 

F pcpc ’ F PV py , V PVPC , where the subscripts indicate the nature of the pion exchanges (parity 
violating or parity conserving). 
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Since the PVPV interaction is mediated only by charged pions, the only possible 
interaction is the above. The proton turns into a neutron when it emits the jt* in order to 
conserve charge. The calculation of the potential associated with this interaction is very 
similar to the calculation above for the simplified model. There are two differences. 
First, the extra factors proportional to the mass which arise from normalization 
conventions do not occur for spinors. Second, since the proton and neutron have two 
spin states, we must determine how the interactions affect their spins. Since the pions are 
spinless, a nucleon’s spin does not flip when it emits a pion. Thus if the proton initially 
has spin up and the neutron is spin down, then the scattered proton has spin down and the 
scattered neutron has spin up. The effective potential for this interaction is 

f 2 

V(r. - F, ) = — X , where X is an operator which turns protons into neutrons 

2 4jr | F, - r 2 \ 

and neutrons into protons. 

Since particles can change types, we will use second quantization. The operator 
n *(A)P, ('i) + w I(^)/ 7 i (' ; i) turns a P roton at point r, into a neutron at the same point 
with the same spin. So the operator X is equal to 
[ P | (Ti )« t ) + P\ iT\ )«* ) ][ n* (F, )p f (F, ) + n * ( r 2 )p l (F, ) ]. 

The interaction Hamiltonian is therefore 


-m a \r x -r 2 \ 


H PVPV = 


2 8jt . — 


■ 1 p\ (F, )« t (F, ) + p\ (F, ) Wi (r, ) ][n| (F 2 )p f (F 2 ) + ( r 2 )p j (F 2 ) 


This interaction does not seem to mediate a force between two separated 
nucleons, since the nucleons can’t retain their identity. In particular, it cannot mediate a 
force between two separate nuclei, even though it does contribute to the binding energy 
of a single nucleus. Suppose the nucleus can be modeled as consisting of two Fermi seas, 
one of protons and one of neutrons. (Of course H PCPC , which is much larger than any 

parity violating interactions, produces correlations in nuclei.) There are protons of both 
spins and of every momentum less than k p , and there are neutrons up to k n . We will 
assume the nucleus is very large so that we can approximate the waves by plane waves, 
and we will calculate the energy density due to H pvp Then since the proton and 
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neutron which interact are in the same nucleus, the interaction does not change the 
nucleus type. 

The energy shift is <il H pypv li>. Consider the action of H pvpv on the Fermi sea. 

It acts on two particles — a proton of momentum k , and a neutron of momentum k -, , 
changing their momentum and interchanging their particle types. Since the resulting state 
is multiplied by <il, it must have the same occupied states as li> in order to give a nonzero 
contribution to the energy. Thus the final proton must have the same momentum as the 
initial proton, as in the following diagram: 


P fC * OT1 proton 



2 2 


Also, the particles must have the same spin, since the interaction does not flip spin. 

In order to see at least one failure of our approximations, we consider two general wave 
functions for a proton and neutron and calculate their interaction in these states. We 
ignore the spin wave function, and set 

I ^ >= )/a ( p 2 > I P n> -fii p \ )f\ ( p 2 )\np>). 

The wave function here is antisymmetric in space and isospin, according to the Pauli 
exclusion principle for multiple types of fermions. Thus neither rj nor r 2 is the 

coordinate of the proton. But f, must be the wave function of the proton. Assuming both 
spins are up, 

Hpypy \y >= -x 2 )d 3 x t d 3 x 2 p^(x 1 )n 1 {x ) )n*(x 2 )p 1 (x 2 )[f ] (F l )/ 2 (F 2 )lpn > -f 2 (F t )f t (F 2 ) \ np >] 

)f\<j\)f 2 {F 2 )\ np> - V(F, -F 2 )f 2 (F i )f i (F 2 )\ pn >]. 

In evaluating the integrals, we must choose *, and x 2 to coincide with F } and f 2 in one 
order, so that x, represents the location of the initial neutron. Now, noticing that 
<pnlpn>=<nplnp>=l while <pnlnp>=<nplpn>=0, we obtain 

< V i H PVPV |V’ ^4jpV 3? 2^-^)[/i(n)*/2(f 2 )*/2('i)/i(f2) + /2('i)*/i(' : 2)Vi(n)/2(' : 2)] 

= 4-JpV^ v ( p } - h )./)('■] )/] (F 2 )f 2 (F l )*, 
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where the two terms can be combined by exchanging r, and ? 2 in the second term. 


This has the form of an exchange force— notice that if f, and f 2 are never nonzero at the 
same point, then the integrand vanishes. This is why we claimed above that the shell 
structure might affect the parity violating energy; a simple form with plane wave states 
ignores all the structure of the nucleus. 

In order to find the total energy, we must sum over all occupied proton and 
neutron states f, and f 2 respectively. There is also a self-interaction due to H PVPV in which 
a nucleon changes type and then changes back. This gives an infinite contribution to the 
energy, and is also more significant because a particle is certain to be at the same location 
as itself, so the comment about overlap at the end of the last paragraph doesn’t apply. 
We will discuss self-energy contributions at the end. 

Substituting V(r } - r 2 ) = j and plane waves which are normalized so that their 

* I n - r 2 | 

integral over the volume V of the nucleus is 1 , we obtain an energy proportional to the 
Fourier transform of the Yukawa potential. It is equal to 


< H pypy >— 


fn 


2V ni„ 2 + \k i - A 2 


- — — , and the sum over all pairs of nucleons is 


6Epy = 2 


L 


'jl- y \ l 

2V m » 2+ \ k \- k i 


where the factor of 2 arises because the proton and neutron can be both spin up or both 
spin down. Converting the sum to an integral in the usual way, we obtain 


5E Pl , = 


L 


dk.dk-, 


,T , ' \ 2 r r “ "1^ 

y w 3 


First integrate over k.-, , taking the z-axis along k ] . 

r , 3 r 1 „ k " r , 2 „ * c sin 6dd 

f d k, — , = — = — = 2jt f A:, dk , f — , , 5 

J " m/+\ A, -k 2 |- { " { m/ + A,' + A,' -2A,A, cos0 

*’ 2 \n(m„ 2 +A, 2 +A 2 2 -2A,A,cos0) 

- 2jt I k 2 dk 2 lo 

J 2A,A, 

Jrt, Ah w ^ 2 +(^i +k 2 ) 2 
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We now integrate over the direction of k , , giving a factor of 4 ji, and then give the 
expression for the integral: 
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The integral is evaluated by repeated integrations by parts. 

We would now like to evaluate this expression as a function of x = Z/A. The 
volume of a nucleus is proportional to A; that is, each nucleon occupies a fixed volume. 
The volume is given by 4^'n A where r„=1.2 fm xxv . As usual, the Fermi momenta are 
given by k p =(3jr~n p ) ]} where n p is the proton number density, so 

k p - T-tfW* = 320\[x MeV and similarly k n = 320^1 - x MeV . Furthermore, the mass 
of a nucleus is roughly Am N , so 
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(Notice that our 


* N r "‘N 

estimates, for x = 1/2, give a Fermi momentum of 250 MeV, which is somewhat small 
compared to the nucleon mass, so our nonrelativistic approximations (in particular, 
neglecting the time-like component of q in the derivation of the effective interaction) 
seem reasonable.) 


6E 

The following figure shows — — as a function of x; it is on the order of 10 16 and is 

M 

5E 

almost exactly equal to — — = 9x(l - x)xl0~ 16 , which is much easier to interpret than 

M 

the formula given above; for a contact force, the number of interactions occurring in the 
nucleus is proportional to the product of the densities, which is proportional to x(l-x). 
(The suggested derivation of this formula— assume the mass of the pion is very large and 
replace the interaction by a contact interaction — gives an incorrect answer, 
30x(l - x) xl(T' 6 . The mass of the pion is not large enough to regard V pv as a contact 
interaction.) 
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The fraction we have just calculated is very small, and (since we can only 
compare different nuclei) its variation from nucleus to nucleus is even less— 2x10 17 , if x 
varies from 1/2 to 1/3. There are similar smooth empirical formulae for the total mass of 
a nucleus, but the mass differences are somewhat larger from nucleus to nucleus because 
the mass depends in a more jagged way on the atomic number. Thus, the plane wave 
description of the protons and neutrons cannot suggest which nuclei have unusually large 
or unusually small contributions from the weak force. There are many other 
contributions to the parity violating energy as well— for example, the interaction V PC p V 
produces an energy shift at second order in perturbation theory. This can occur only at 
second order because <il V pcpv li> vanishes— after all, li> is a parity eigenstate since the 
strong forces which determine the structure of the unperturbed nucleus are parity 
conserving. However, by introducing an intermediate state of opposite parity and very 
similar energy, one can hope to obtain the largest possible weak parity violating 

5E 

contribution to the energy. For l9 F, Ref xxv ‘ estimates — — =7xl0‘ 17 , which is very 

M 

similar to what we obtained above. Ref. xxvi did not calculate the matrix elements from 
the theoretical formula for V PCPV but from an experiment mentioned in Ref xxv. The 
latter reference also calculated the matrix element from a theory based on exchange of 
pions and other particles, giving a value which was off by only 20% .) 

Ref xxvi’s calculation did not take the direct exchange of single pions in which 
both the emitted and absorbed pion have parity violating interactions, leading to an 
impression that only nuclei with narrowly separated partners of opposite parity have large 
6E 

values of — — . However it leads to a very clear case in which one expects a jump in the 
M 

parity violating energy— choose a nucleus li> which is close in energy to another nucleus 
li’> with the opposite parity. Fluorine is not the only choice. If one were to calculate the 
matrix element of V PCPV between these states theoretically (the calculation would be 
similar to the above), one would find that < 1 I Vpcpv I ‘ >oc vo ^ ume a M . Thus, 
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&E pi ] !< 1 I Vpcpv 1 1 > l „ M . . 

— = — — — « — — , where the denominator here is just the energy 

M Mj-Mj. 

denominator (with energy replaced by mass). (M is either of the two masses.) So we 
expect an especially large parity violating energy if the mass difference of the nucleus 
and its excited state is small compared to the mass of the nucleus. 

Finally, Ref. ” v " calculates the weak-force self-energy of a proton or neutron and 
finds that its fractional contribution to the mass of a nucleus is 10 times larger than the 
energy from the interactions of different protons and neutrons in the nucleus calculated in 
Ref. xxiv. Most likely, the parity violating self-energy of the protons and neutrons is 
more significant than anything we have already calculated. This returns to one of the 
most obvious models for equivalence principle violations — assume that gravity couples 
differently to protons and neutrons. This gives gravitational mass differences which are 
linear in x (just as we found in the model of a long range force coupled to the total 
number of baryons in a nucleus). So a simple approach is to choose nuclei with the 
largest and smallest possible values of x. (Hydrogen and Uranium are suggested in Ref. 

xxviiij^ 

There are three reasons why it is unlikely that gravity couples differently to the 
“parity violating part of the weak force.” First, this depends on whether one uses the 
fundamental Lagrangian or the current-current Lagrangian as discussed above. In 
practice, the definition of the parity violating and parity conserving energies is very 
technical— a summary of the definition we used is “it is the part of the weak force which 
is second order in perturbation theory.” Second, in the electroweak theory, photons are 
linear combinations of a W particle and another gauge particle, so if gravity couples 
differently to the W particles, it probably couples differently to the photon as well. Third, 
all self-energies are infinite in any case. This is why the self-energy of the protons and 
neutrons calculated in Ref. xxix is much larger than the energies of the nuclei calculated 
in Ref. xxiv. The interactions of the quarks inside the protons and neutrons are larger 
than the interaction energies of different protons and neutrons because the quarks are 
closer together. If the quarks are point particles, then they contribute infinite self- 
energies due to the weak force besides these other energies. Fourth, energy cannot be 
separated into different forms of energies. For example, the energy we calculated above 
includes the mass of the virtual pions, which is due to the strong force. 


Addendum to Theoretical Impact Analysis 

In two recent papers"^ “\ Damour et al., provide a theoretical justification, based on the 
string theory for expecting violations of the Equivalence Principle at an accuracy above 
one part in 10' 12 which is slightly below the present state of the art in EP tests. Their 
argument centers on the inflationary rate of expansion of the early universe and on the 
decoupling of scalar fields from matter during the expansion of the universe. When the 
recent estimates of the early expansion rates are factored into the theory, they conclude 
that violations of the EP may well occur at accuracy much lower than previously thought. 
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The conclusion of the analysis by Damour et al. gives strength to EP experiments with 
improved accuracy even if the improvement is not as dramatic as those promised by the 
space-based tests (of several orders of magnitude) of the EP. Another implication is that 
our goal should not be to push the accuracy as high as possible (and consequently 
increase the difficulty and cost of the experiment) but rather achieve a significant 
improvement in accuracy with respect to the present state of the art. 


4 
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Capsule/Cryostat Design Approach 


Introductory Remarks 

During the early conceptual development of this experiment, we were focusing our 
attention mostly on free-fall spans of a few meters to achieve free-fall time as long as 
possible. This situation, however, precluded the use of a fully cryogenic capsule because 
of mass considerations and pushed us towards the adoption of a small cryostat at the top 
of the vacuum chamber to refrigerate the instrument before release. However, the small 
cryostat creates additional problems which are highlighted by the previous analyses as 
follows: (1) the small cryostat is a concentrated mass placed very close to the detector 
and, cosequently, produces relatively strong gravity gradients; and (2) the instrument 
package spends most of the free-fall time in the vicinity of the small cryostat and moves 
away from the small cryostat only in the later portion of the fall. 



Figure 34 Vertical motion of package relative to capsule for p o = 7000 kg//m 3 (shown 
for distances < 1 m) 

Figure 34 shows clearly that the relative distance between the instrument package and the 
capsule (for a typical p o = 7000 kg/m 2 ) is less than 10 cm during the first 14 s, increases 
to about 1 m at the 24 s mark and becomes greater than 3 m at the 30 s mark (not shown 
in the figure). In other words, the small cryostat can not be small because it needs to be 
moved farther away for reducing the self-generated gravity gradients and also for 
allow ing for lateral motions of the capsule with respect to the falling package. The free 
fall time (once several seconds are spent in attenuating the initial transient motion of the 
sensing masses after release) affects the experiment accuracy only through its square root 
(i.e., a weak function). Consequently, it is sensible to sacrifice a few seconds of free-fall 
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time in order to make the experimental area sufficiently short to open up the option of a 
fully cryogenic vacuum chamber for the detector to fall in. 

These considerations are formalized in the next section by developing an optimization 
process aimed at identifying the size and mass of the experimental chamber and the 
capsule. 

Optimization of Capsule Size and Mass 

A design optimization would appear premature at this point of the flight definition study. 
This optimization process, however, is aimed exclusively at evaluating in a consistent 
way the size of the capsule and its overall mass. The goal of the optimization process is 
to devise a mathematical formulation for evaluating capsule sizes and masses that satisfy: 
(1) geometrical/dynamical constraints and (2) engineering constraints. The first category 
includes the constraints determined by the strength of the gravity gradient inside the 
vacuum chamber (generated by the capsule mass itself) and related to the motion of the 
instrument package with respect to the capsule during free fall. The second category 
includes size constraint due to transportability and handling of the capsule and also 
minimum and maximum allowable mass. The minimum mass is determined by the 
ability of building a vacuum chamber that does not buckle under the atmospheric 
pressure while the maximum mass is limited by the carrying capacity of the balloon. 

The independent variables to work with (as it will become clearer in the following) are 
the (frontal) low-speed ballistic coefficient of the capsule and the free-fall time. The 
output variables are the capsule diameter and related capsule length, free-fall span, 
cryostat mass, and capsule mass. 

In more detail, we first need an expression of the free-fall span of the package inside the 
capsule as a function of the ballistic coefficient |3 0 and the free-fall time T f . The free-fall 
distances were derived numerically because the relevant functions are not integrable. 
Consequently, we have to fit the numerically-derived function over the two parameter 
space ((3 0 , T f ). The fitting was done successfully over the range of interest of the 
parameters by starting from the analytical formulation that approximates the solution for 
a constant ballistic coefficient and then evaluating the fitting coefficients for the general 
case with variable C D . The two-parameter fitting process yields: 

z f = aT f b exp(cTj 2 )/ P 0 (45) 

where a = 0.149636, b = 3.084, and c = 0.001692 are the fitting parameters. The results 
are accurate within a few percent for ballistic coefficients varying from 2000 kg/m 2 to 
20000 kg/m 2 and free-fall times from 0 to 30 s. 

The free fall vertical distance readily determines the minimum length of the 
experimental chamber which must at least equal the free-fall span plus the vertical size 
(outer diameter) of the instrument package, that is: 

Lm.n = TIZf + D (46) 


69 


where D is the outer diameter of the instrument package and q is a safety coefficient 
(typically equal to 1 .5) that takes into account the uncertainties in estimating the drag 
deceleration due to the rarefied atmosphere. 

The next step is to evaluate the minimum internal diameter of the chamber/cryostat 
that can accommodate the lateral and attitude dynamic of the capsule and also provide 
low values of the gravity gradients along the free-fall trajectory. In order to keep the 
equations manageable we adopt for this analysis the constant stand-off distances which 
were derived previously and boundary values for the capsule attitude dynamics, rather 
than using an (on-line) numerical computations of the gravity-gradient field as a function 
of the capsule position (as done in the previous section). The goal is to compute a 
minimum value of the internal diameter and external diameter (where the two are related 
through the cryostat thickness) that keep the free-falling package close to the capsule 
center-line (far away from the cryostat walls) under the action of a maximum wind shear 
and for reasonable assumptions about the capsule attitude (pitch or roll) dynamics during 
the fall. The point above can be translated into the following equation for the external 
capsule diameter: 

^min — pilch y shear 4 gg A'ctvo ) (47) 


y P «ch = z fSin(0 max ) where z f is given by eqn. (45) and 0 max can be interpreted as either a 
limit value of the maximum pitch motion of an uncontrolled capsule or a requirement 
imposed on the control system of the capsule attitude during the fall. The lateral 
displacement due to wind shear y shear =f(0 o , T f ) is given by eqn. (10) while y cryo is the 
thickness of the cryostat wall of 10-15 cm (i.e., typical values for large-size helium 
cryostats). The stand-off distance y gg is the upper bound of the minimum distance 
between the detector CM and the heavy walls of the cryostat that defines an area 
(cylindrical and centered at the cryostat center-line) where the gravity gradients generated 
by the cryostat are sufficiently low. An upper bound of 0.4 m can be adopted across a 
variety of cylindrical cryostats for sizes and masses of interest to this project. 

One obvious consideration is that the internal diameter must be large enough to contain 
the instrument package. This constraint, however, is already accounted for by the fact 
that 2y gg > w P where w P is the width of the instrument package that is presently estimated 
at about 40 cm. The overall length of the capsule is related to the external diameter 
through the fineness ratio that, as explained previously, needs to be about 1 :4 in order to 
keep low the value of the drag coefficient, especially in the compressible flow regime. 

The next step is to estimate the minimum mass required for the cryostat. Most of the 
cryostat mass is for the metal of the vacuum chamber that has to withstand the external 
atmospheric pressure (at ground level). From engineering formulas for the buckling of 
cylindrical vacuum chambers we obtain that the critical thickness of the vacuum chamber 
walls (at buckling) is given by: 


s cri, JpALI^I ] 25 

R 0.92 E 


( 48 ) 
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where R and L are the radius and length of the vacuum chamber, s cril is the (minimum) 
thickness of the load-bearing shell, E is the Young’s modulus of the material and p e the 
external pressure. The cryostat length is readily computed as: 

L = rfZff + 2y gg (49) 

The minimum cryostat mass is then computed as follows: 

M'ayo = vs crit npD(D I k + L) (50) 

where p is the density of the material, D the cryostat average diameter, v is a load safety 
factor and k is a shape factor that is equal 2 for flat cryostat caps and 1 for hemispherical 
caps. After adopting v = 3 and using eqn. (48) to compute the thickness, we find that 
eqn. (50) provides masses which are in line with actual values of large cryostats (as 
verified with Janis Research). Since the cryostat is the heaviest component of the 
capsule, we can estimate the overall (minimum) mass M crit of the capsule as a proportion 
of the cryostat mass M = yM cryo . 

The actual mass of the capsule M, however, is defined by the ballistic coefficient |3 0 
once the external diameter of the capsule is determined as follows: 


M - PqCpqA 


(51) 


where A is the frontal area of the capsule computed through eqn. (47). 

Consequently, we must verify that M > M crit . If the inequality is true we can simply 
add ballast to the capsule in order to preserve the values of the geometrical variables and 
meet the actual mass requirement. In this case, the vacuum chamber of the ballasted 
capsule will satisfy the strength requirements captured by eqn. (48). If the inequality is 
violated, the total mass determined by the ballistic coefficient is too small to build a 
cryostat that is strong enough to withstand the atmospheric pressure. The violation of 
this inequality will be dealt with in the optimization process as a (strong) penalty on the 
cost function as explained later on. 

Figure 35 shows an example of a contour plot of two of the key output variables, that 
is, the capsule mass and the capsule external diameter vs. the low-speed ballistic 
coefficient and the free-fall time for realistic values of parameters. 
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Figure 35 Contour plot with grid lines showing capsule mass and external diameter vs. 
low-speed ballistic coefficient and free-fall time. 


The goal of the optimization process is to define the size and mass of the capsule that 
minimize the value of a cost function. The cost function consists of the overall capsule 
mass, the capsule diameter (which is related to the overall capsule length), and the free- 
fall time. The cost function for this optimization problem can be expressed (and the 
choice is not uniquely defined) as follows: 


CF = 


a,M + a-,D, 




— + p + p 

+ r M + r D 


With '/ * 'min 


( 52 ) 


where M and D ext are the overall mass and external diameter of the capsule, z f is the free- 
fall span inside the capsule (which defines the cryostat length), t^,, the minimum desired 
free-fall time, t f the free-fall time, a, and a 2 are weighting coefficients and P M and P D are 
penalty functions. P M is the penalty function for the total mass and P D for the external 
diameter. 


The structure of the cost function is readily understood after considering that the 
experiment accuracy depends on the square root of the integration time. The minimum 
free-fall time has to include the time required to damp initial transients (a few seconds) 
and also a minimum number of cycles needed for a reliable extraction of the signal from 
noise. The weight coefficients adjust the relevant weights of the two quantities at the 


72 




numerator by making them comparable and/or making one component more important 
than the other one. By adopting meters for the units of length and metric tons for the 
units of mass, the two weighting coefficients have the same order of magnitude. 

Going back to the penalty functions, obviously the overall mass of the capsule must be 
greater than the minimum (or critical) mass M crit determined previously and smaller than 
the maximum mass M max than a reasonable size balloon can carry to an altitude of 40 km. 
Note that M max is a fixed value while the critical mass depends upon the optimization 
parameters. 

In a similar manner, the external diameter (which also determines the overall length of 
the capsule according to L c = D ext /6) must be large enough to accommodate the 
instrument package plus the cryostat walls (i.e., larger than D min ) and smaller than a 
maximum value D max . In this case both values are fixed and moreover the first inequality 
(larger than) is superseded by eqn. (47) if 2y gg > w P . 

The penalty functions are bowl-shaped functions which are equal to zero for values 
within the two boundaries of the inequality and whose value rapidly increases as the input 
variable approaches the boundary values. Figure 36 shows the mass penalty function for 
two values of M crit = 0.3 ton and 1.0 ton, respectively, and M max = 2.2 tons. The 
maximum value of the penalty function is quite arbitrary. The value simply needs to be 
substantially greater than the expected range of values of the first term that appears on the 
right hand side of the cost function in eqn. (52). 



Figure 36 Example of mass penalty function with M cri , = 0.3 ton and 1.0 ton, and 
M max = 2.2 ton. 

A similar penalty function with fixed-valued left and right boundaries was derived to 
penalize the external diameter of the capsule. Without going into excessive details, the 
penalty functions are obtained from a pair of exponential functions with a set of 
parameters and power exponents that shapes the function appropriately. 
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Optimization Results 

The optimization routine has been exercised on a large number of cases in order to 
inspect the influence of any given input parameter on the output/design parameters and 
the cost function. In the following we will show only the most important results relevant 
to: (a) a cryostat made of Aluminum and (b) a cryostat made of Steel. The numerical 
parameters adopted for computing the results shown in Figure 37 and Figure 38 are as 
follows: instrument package width w p = 40 cm; low-speed drag coefficient of the capsule 
Cuo = 0.1; fineness ratio of the capsule 6 = 0.25; wind shear = 0.005 m-s '/m; wall 
thickness of cryostat = 15 cm; stand-off distance for reducing gravity gradients y gg = 40 
cm; uncertainty factor on drag deceleration r| = 1.5; strength safety factor for vacuum 
chamber v = 3; mass ratio capsule-mass/cryostat-mass y = 2; density of material p = 2800 
kg/m 3 (Aluminum) and 7800 kg/m 3 (Steel); Young’s modulus E = 80 GPa (Aluminum) 
and 200 GPa (Steel). Moreover, the maximum capsule mass for a mass penalty in the cost 
function is equal to M max = 2200 kg. The maximum external diameter of the capsule for a 
penalty is = 2.2 m (which implies a maximum capsule length of 8.8 m). The weight 
coefficients have been as follows: a, = 1 (mass) and a, = 1 .5 (external diameter) in which 
masses are expressed in metric tons and diameters in meters. More weight has been 
given to the diameter rather than to the total mass because, as it will be seen later on, 
there is a fairly large mass margin (large balloons can carry a few tons) while a large 
external diameter leads to long capsules that are difficult to handle and transport. 

The value of the cost function is shown as black contour lines in Figure 37 and Figure 
38 and the contour regions are colored with the deep blue color corresponding to the 
minimum of the cost function. The grid lines of the free-fall span are also added to the 
figures in green color. The figures provide a host of useful information as the grid lines 
specify the values of key variables such as capsule mass, external diameter, and free-fall 
span (of the instrument package inside the capsule) for any point identified by the 
coordinates j5 0 and t f . The equations shown previously enable us to compute readily other 
key design parameters: the cryostat length L from eqn. (49) and its mass by dividing the 
capsule mass by y. 

Figure 37 shows the results for a cryostat made of Aluminum and Figure 38 for a Steel 
cryostat. The cost function identifies the area on which to focus the selection of the key 
design parameters. In the (desirable) deep blue area (inside the smaller circular contour 
lines) the cost function is rather shallow and, consequently, the region for selecting the 
design parameters is reasonably large. In the dark color region, one of the mass or 
geometrical constraints has been violated and consequently the penalty functions bound 
the design area. Note that the maximum capsule mass limits the design area on the 
right-hand side of the plots while the (critical) cryostat mass and the maximum external 
diameter of the capsule limit the design area on the left side of the plots. In the Steel- 
cryostat case the (critical) cryostat mass is the leading boundary condition on the left side 
of the design area. In the lower side of the plot, the design area is limited by the 
minimum free fall time that has been set equal to 1 5 s in order to allow for a reasonable 
numbers of signal cycles. 
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*Blue - Tot. mass (tons); ¥Red - Ext. dia. (m); UGreen - FreeFall span (m); §Black - Cost function 
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Figure 37 Cost function in the capsule parameter space (Aluminum cryostat) 

The figures show that an external diameter in the neighborhood of 1 .4 m is a valid 
selection for the capsule based on the assumptions of this analysis. The 1.4-m diameter 
contour line cuts across the region of the minimum value of the cost function. 
Furthermore, the minimum-value region is larger for an Aluminum cryostat than for a 
Steel cryostat because in the latter case the boundary of the minimum (and critical) 
cryostat mass encroaches into the low-valued area of the cost function. Consequently, 
Aluminum is preferable to Steel, from this analysis point of view, as it allows wider 
margins to work with in the capsule design. However, Steel is also a valid choice and it 
might be preferable to Aluminum for the vacuum chamber because of the low outgassing. 
After focusing on Figure 37, an “optimal’’ design choice could be an external capsule 
diameter of 1.4 m (with an overall capsule length of 5.6 m). The overall capsule mass 
could start at a minimum value of about a 1000 kg which provides a free fall time of 23.5 
s. Note from the figures that the free-fall span contour line is about parallel to the 1 .4-m 
external diameter line. This situation implies that we could move the design point along 
those two lines by increasing (with ballast) the capsule mass and increasing the free-fall 
time. Consequently, a cryostat of given length can provide longer free fall times by 
simply making the capsule heavier with ballast. The free-fall time can be increased by 1 
s for every 250-kg ballast (or other equipment) added to the capsule up to the point where 
the upper mass limit is reached. The low-speed ballistic coefficient (3„ of the capsule will 
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vary from about 6500 kg/nr for a capsule mass = 1000 kg to 13.000 kg/nr for a capsule 
mass of 2000 kg. 


*Blue - Tot. mass (tons); ¥Red - Ext. dia. (m); IJGreen - FreeFall span (m); §Black - Cost function 
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Figure 38 Cost function in the capsule parameter space (Steel cryostat) 

The previous choice of parameters is also valid for a Steel cryostat with the only 
difference that the design point for the lower value of ballistic coefficient (i.e., 6500 
kg/m : which is associated with a mass of 1000 kg) is closer to the left boundary, defined 
by buckling. The design point moves away from that boundary for higher values of the 
low-speed ballistic coefficient. 

This ‘‘design concept’* optimization does not intend to exhaust the design options for 
the capsule which will have to be revisited many times before the design is frozen based 
on additional system analyses. However, this process is a valid starting point to provide 
preliminary numbers for the analysis to be conducted by our cryostat subcontractor. The 
optimization process also provides a framework that lends itself to future refinements for 
evaluating system design parameters vis a vis performance input requirements and 
engineering constraints. 
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Updated Reference Configuration and Drop Time 
Capsule mass and Size 

One of the important results of the analysis carried out is that the small vacuum 
chamber at the top of the capsule has been eliminated in favor of a fully cryogenic 
vacuum chamber (see Figure 39). The new solution eliminates the problem of the gravity 
gradient produced by the small cryostat in the proximity of the detector and also provides 
more clearance to the instrument package during the early stage of free fall. 



Spin and release 
mechanism 


Instrument package 
before release 


Cryostat 


Figure 39 Schematic of capsule with instrument package attached before release 
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The new configuration also has the advantage of using a more conventional design of 
the cryostat than the old reference design with the small LHe cryostat. 

The preferred size of the vacuum chamber/cryostat was identified as a cylinder with 
internal diameter of 1 m and an internal height (at the cylinder edges) of 2.1 m which 
results in an external diameter for the capsule of about 1 .4 m and a free-to-chose capsule 
length of roughly 5.6 m. Preliminary data on large cryostats (from Janis Research), 
indicate that the mass of a cryostat of the size considered above will fit well within the 
mass limit of the system. Such chamber/cryostat will allow free fall times in the range 
24-28 s depending on the amount of ballast added to the capsule. 

The overall capsule mass could start at a minimum value of about a 1000 kg which 
provides a free fall time of about 24 s. The free-fall time can be increased by 1 s for 
every 250-kg ballast (or other equipment) added to the capsule up to the point where the a 
(limit) mass of 2000 kg is reached. The low-speed ballistic coefficient p n of the capsule 
will vary from about 6500 kg/m 2 for a capsule mass = 1000 kg to 13,000 kg/m 2 for a 
capsule mass of 2000 kg. Correspondingly, the free fall time will vary from about 24 s to 
28 s. 

Capsule Drop Time 

The free fall time of the instrument package inside the capsule was computed for several 
values of the low-speed ballistic coefficient in the previous sections. Figure 40 depicts 
(again) the capsule drop distance and Mach number vs. drop time. 



Figure 40 Capsule drop and Mach number vs. drop time 
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Following the computation of the preferable range of values for the capsule mass and 
size. Figure 41 shows the relative distance covered by the instrument package during free 
fall with respect to the capsule for a low-speed ballistic coefficient ranging from 6500 
kg/m 2 to 13500 kg/m 2 . 

These results clearly indicate that a capsule of relatively compact size can provide a 
free fall time between 24 s and 28 s with a capsule mass ranging from 1000 kg to 2000 
kg, respectively. The higher value of the mass is for the same capsule with ballast. 



Time (s) 


Figure 41 Relative vertical distance vs. free fall time for limit values of interest of 
ballistic coefficient. Dots mark the expected operating points for the light capsule (M = 
1000 kg) and ballasted capsule (M = 2000 kg). 


79 



Figure 42 Relative lateral displacement due to wind shear vs. free fall time. Dots mark 
the expected operating points for the light capsule (M = 1000 kg) and ballasted capsule 
(M = 2000 kg). 

An internal diameter of the capsule of 1 m was also computed from the optimization 
process. This diameter determines the tolerance to wind shear acting sideways on the 
capsule. Figure 42 shows the lateral displacement of the capsule with respect to the free- 
falling instrument package vs. free fall time for the limit values of interest for the low- 
speed ballistic coefficient. The value adopted for the wind shear of 0.005 s' 1 is equivalent 
to a vertical gradient of 10 knots per km. This value is twice as high as the maximum 
wind shear reported’™ for the Air Force balloon base at Holloman, New Mexico. The 
lateral displacements due to wind shear are relatively small in the parameter range of 
interest. Furthermore, if the balloon is launched during the periodically-occurring wind 
reversal times (in April-May and September-October) the vertical wind gradient is much 
smaller than the value adopted for the computations shown here. 
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Modal Analysis of Detector 


Introduction 

The following is the modal analysis of the differential accelerometer for the EP test. 
The sensor is composed by two masses of different materials constrained to an external 
case by means of elastic springs. The measurement is performed by monitoring the 
displacement of the two masses with respect to the external case through the 
measurement of capacitance variations caused by the displacement. 

This analysis evaluates the eigenfrequencies and eigenmodes of the detector and the 
sensitivity of the modes to variations in the detector parameters. 

The model adopts six degrees of freedom, that is, one rotation and one translation for 
each body. The system parameters considered for the sensitivity analysis are masses and 
moments of inertia of the bodies, stiffness of the springs and their ratio. 

The influence of these parameters is evaluated on both the separation of the natural 
frequencies and the vibrational modes, particularly for what concerns the differential and 
common modes of the detector. 


Detector Dynamics Model 

The detector consists of two sensing masses having their Centers of Mass (CM) 
coincident with the CM of the external case. In the current detector model, one mass has 
a dumbbell shape, while the other is a hollow cylinder. The shape of the flight detector 
sensing masses will likely be different from this early choice. However, the shape of the 
sensing masses does not affect the generality of this study. The two sensing masses are 
constrained to the case by means of elastic springs, having a pivot axis parallel to the 
longitudinal axis of symmetry of the three elements. Figure 43 shows a schematic of the 
sensor model. 



spin axis 


pivot axis 


Figure 43 Schematic of detector model. 
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For this analysis, we assume that the capsule containing the sensor is fixed to the 
Earth surface, so that only the sensor is falling into it. This assumption does not affect the 
generality of the modal analysis. The Z axis is the axis of fall and of the Earth’s gravity 
acceleration. 

The reference system Cf is a frame parallel to the geographic X, Y, Z axes and 
displaced by the vector r z (see Figure 44). Another coordinate system D is attached to the 
capsule and having axes parallel to the local normal frame. At this stage of the analysis, 
reference frame D is assumed to coincide with the local normal frame and the capsule is 
assumed stationary with respect to the Earth’s surface. This assumption doesn’t cause any 
loss of generality for the modal analysis. The coordinate that describes the fall of the case 
along the z axis is r 2 (t), where r is the position vector that identifies the body C (external 
case) with respect to the capsule and, finally, z indicates the axis of fall. 



Figure 44 Local normal (geographic) and Cf frames 

The body coordinate systems x, y, z are centered at the body CM, and fixed with the 
body itself. Each system has the x axis along the longitudinal axis of symmetry, 
corresponding also to the spin axis while the y and z axes are radial axes rotating with the 
body. 

A rotational and translational degrees of freedom have been considered for each 
mass. More precisely, the x rotation about the CM of each mass, and the y translation 
referred to the xyz frame. Therefore each body can rotate about the x axis through its CM 
and can translate along the y axis of our reference frame. Figure 45 shows the possible 
movements of each body. 
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Figure 45 Degrees of freedom of each body 


External forces 


The motion of the sensing mass consists of a rotation 0(t) about the pivot axis (related 
to the torsional stiffnesses of the springs) and a translation y(t) that depends on the 
flexural stiffness of the springs. The flexural stiffnesses are typically much greater than 
the torsional ones. The distance d in the figure identifies the distance between the CM of 
the sensing mass and the pivot axis. 

The external forces are written in the Cf reference frame, and consist of elastic forces. 
The elastic forces are produced by the springs connecting the sensing masses to the 
external case, and consist of a torsional torque and a force. In this analysis we assume 
that the rotational displacements are small, so that the displacement of the point S along 
the y axis due to a rotation around the CM is simply equal to 6d [i.e., we assume that 
cos(0) = 0]. 

The differential displacement of the point S between body A and body C and the 
elastic force along the y axis are therefore: 

Ay A = y A - 0A-d A - (y c - 0 c -d A ) 


FyA = -k yA Ay A = - k yA (y A - 0 A -d A - y c + 0 c d A ) 

Where y A and y c , are the displacements along the y axis with respect to the Cf frame, 
0 A and 0 C are the rotations about the x axis of masses A and C, respectively. k yA is the 
flexural stiffness along y of the spring connecting body A to the case C. 

The elastic torque about the CM of A is: 
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N a = -k tA (0 A - 0 C ) - ky A ’d A Ay A = -k w (0 A - 0 C ) - k yA d A (y A - 0 A *d A - y c + 0 c d A ) 
Where k tA is the torsional stiffness of the spring connecting A and C. 

Similarly, for the other two bodies we have: 

FyB = -k yB *Ay B = -k yB (y B - 0 B -d B - y c + 0 c -d B ) 

N B = — k^ (0 B — 0 C ) - k yB *d B *Ay B = -kt B (0 B — 0 C ) - ky B *d B (y B - 0 B *d B — y c + 0 c *d B ) 
fyc = -k yA (y c - 0 c *d A — y A + 0 A *d A ) — k yB (y c — 0<fd B — y B + 0 B *d B ) 

N c — — k^ (0 C — 0 A ) — k yA *d A (y c — 0c*d c — y A + 0 A *d A ) — k lB (0 B — 0 C ) — ky B *d B (y B — 0 B *d B 
-yc + 0c*d B ) 


These equations will be used in the next section to compute the stiffness matrix. 

Modal Analysis 

The mass and stiffness matrices are computed in this section. It has to be noted that 
out of 18 degrees of freedom (DoF) associated with the translational and rotational 
motion of the three masses only six of them are describing the motion of the degrees of 
freedom of interest. Specifically, only the translations along the y axis and the rotations 
about the x axis have to be considered for each body because the other motions do not 
affect the measurement. The state vector involved is, therefore, as follows: 

x = [>'a y B yc 6 A e B e C f (53) 

which contains the translations along y axis and the rotations about the principal x axis of 
each body. 

The mass matrix is a diagonal matrix in a system of rigid bodies as follows: 

m A 0 0 0 0 O' 

0 m B 0 0 0 0 

,, 0 0 mr 0 0 0 

M = ( 54 ) 

000/^0 0 K ’ 

0 0 0 0 I Bxx 0 

0 0 0 0 0 I Cxx 

The stiffness matrix, obtained from the force and torque equations by using the 
displacement method, is: 
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kyA 

0 

kyA 

-d A kyA 

0 

dAkyA 

0 

kyB 

~ k yB 

0 

-d B k yB 

d B k yB 

~kyA 

~k yB 

kyA + kyB 

d a kyA 

d B k yB 

~ d AkyA - d B kyB 

~^AkyA 

0 

d A ky A 

k lA + d A k yA 

0 

~k/A ~ d\kyA 

0 

-d B k yB 

d B k yB 

0 

2 

k tB + d B k yB 

-ktB ~ d B k yB 

d A kyA 

d B k yB 

-dAkyA ~ d B k yB 

2 

-k(A — dAkyA 

~k tB - dftk y ft 

k t A + k tB + d A k yA + d B k yB 


(55) 

We can observe that for an external body that is not constrained in its free fall, two 
rigid motions are possible, that is, the rigid body rotation around the x axis of the Cf 
frame, and the rigid body translation along the y axis. In order to carry out an analysis of 
the natural frequencies and natural modes of the system we can eliminate these two rigid 
motions so as to preserve only the elastic modes. This can be done by expressing one 
rotational and one translational coordinate of a body as a function of the corresponding 
coordinates of the other two bodies. 

A vector having the three translations equal to a constant and the three rotations equal 
to zero is a solution of the eigenproblem, so that this vector is an eigenvector. 
Consequently, all the others eigenvectors must be orthogonal to this one to yield: 

{<& 0 } T [M]{<D} = 0 

where <f>„ is a rigid eigenvector, and O is another eigenvector. After evaluating the 
matrix products we obtain: 

O 0 (m A <I>| + m B d> 2 + m c <J> 3 ) = 0 

where now cjx is the i-th component of a general eigenvector. Because O; = fx h where Xj 
is the i-th component of the vector x and f is a constant equal for all the components, we 
can also write 


m A y A + m B y B + m c y c = 0 

which allows to express one translational coordinate of a body as function of the 
translations of the other two. It can also be noted that in deriving the former equation we 
have obtained the conservation of the linear momentum for any elastic motion. 

We can obtain similar expressions for the rotational case as follows (from now on we 
drop the subscript “xx” in the moments of inertia): 

I A I B ® B "*■ Ic = 0 

We can now express the coordinates of a body as function of the analogous 
coordinates of the other two bodies. We can write therefore 

{x} = |C]{x c } 
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where x is the coordinate vector, C is a matrix, and x e is the new coordinate vector having 
the two coordinates of a body expressed as function of the others. In our case, we have 
eliminated the coordinates of body C, obtaining the following matrix [CJ: 


C = 


1 

0 

m A 

m c 

0 

0 

0 


0 

1 

m B 

m c 

0 

0 

0 


0 

0 

0 

1 

0 



0 

0 

0 

0 

1 



( 56 ) 


and the vector X e being 


>'a 



d B 


( 57 ) 


The stiffness and mass matrices can be modified using the new coordinates by 
[KJ = [C] t [K] [C] 

[MJ = [Cf[M] [C] 

where M e and K c are the transformed mass and stiffness matrices which are not shown 
here for the sake of brevity. The dynamical matrix [DJ = [MJ '[KJ that is needed to 
compute the eigenfrequencies of the elastic modes is as follows: 
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D = 


, i 1 1 

Ka\ + — , 

' \ m A "CJ 

Kb»’A 


m B tn C 

dA k yA( m A +m c) 

Ia™c 

d B k yB m A 

lB ni C 


Ka™b 


Kb 


m A m B 

— + — 

\m B me) 

d A k xA m B 


I A mc 

d B k x[} (m B + me) 

I B l11 C 


_d A ky A Uj\+J^ 

l C m A 

^bKb 1 A 

Ic m B 

{ l A + 1 c)[Ka + ^aKa) 

kJ~c 

1 a(Kb + ^bKb) 

l B l C 


<KKa]_b 

Ic m A 

Ic m B 

1 b[Ka +cI aKa) 

TJc 

(Jb + + ^bKb) | 

l B l C 


( 58 ) 


Once solved the eigenvalue problem for the constrained system (i.e. with the previous 
dynamical matrix), the total eigenvector including also the two variables eliminated can 
be obtained simply premultiplying the eigenvector obtained by the matrix [C]. 

Eigenvalues and Eigenvectors. 

Using the dynamical matrix evaluated before an analysis of the eigenvalues and 
eigenvectors have been carried out in order to evaluate the influence of parameters on the 
natural frequencies and modal vectors. 

The parameters are the following: 

M a = mass of body A 

M b = mass of body B 

M c = mass of body C 

I A = moment of inertia of body A about the spin axis (x axis) 

I B = moment of inertia of body B about the spin axis (x axis) 

I c = moment of inertia of body C about the spin axis (x axis) 
d A = distance between the CM and the constrained point of body A 
d B = distance between the CM and the constrained point of body B 
k tA = torsional stiffness of spring connecting bodies A and C 
k lB = torsional stiffness of spring connecting bodies B and C 
k yA = bending stiffness of spring connecting bodies A and C 
k yB = bending stiffness of spring connecting bodies B and C 
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Numerical values have been assigned to these parameters to evaluate natural 
frequencies and eigenvectors. Subsequently, the values of key parameters have been 
changed while keeping the others constant. 

The values used for the first evaluation are as follows: 

m A = m B = 5 kg 

I A = I B = 0.009 kg-m 2 (which corresponds to an external radius of body A of 0.06 m 
and a mass m A ) 


I c = 1 kg-m 2 


d A = d B = 0.06 m 

k tA = k,B = 9.593 Nm/rad (corresponding to a first torsional frequency for a single- 
degree-of-freedom system of 3 Hz) 

kyA - 18000 jt N/m (corresponding to a first lateral frequency for a single- 
degree-of-freedom system of 30 Hz) 

Using these values we obtain the following values: 

f, = 2.99 Hz; f 2 = 3.29Hz 

f 3 = 52.13 Hz; f 4 = 55.22 Hz 

and the associated eigenvectors: 

v, = {0.0427736, -0.0427736, 0.705812, -0.705812} 

v 2 = {0.0326524, 0.0326524, 0.706352, 0.706352} 

v, = {-0.0209932, 0.0209932, 0.706795, -0.706795} 

v 4 = {-0.0210127, -0.0210127, 0.706795, 0.706795} 

Remembering that the eigenvector components represent in order y A , y B , 0 A , 0 B , we 
can make a few considerations about the natural modes. 

The first vector is composed of a translation of the masses and a rotation about the 
CM. The sign is opposite for the two masses, so that one moves in opposition to the 
other while the signs of y A and 0 A are the same so that the pivot point tends to remain in 
its original position. In the first mode body C is at rest, the motion being auto- 
compensating. In summary, the first mode corresponds to a differential mode of the two 
test masses. 


88 




Figure 46 First natural mode (differential) 

The second vibrational mode has a similar type of motion, the only difference being 
that the two proof masses move in the same direction. In this mode body C moves in the 
opposite direction of the proof masses in order to equilibrate the motion. 



Figure 47 Second natural mode (common) 

The third and fourth modes involve the lateral stiffness of the springs as the sign of 
the rotation and translation of each mass is now opposite, so that the displacement of the 
constrained point is now the sum of the y translation and of 0d A term (i.e., y A - 0 A d A ). The 
y translation depends on the lateral stiffness of the springs that is much higher than the 
torsional stiffness. The first two modes are analyzed in detail in the following. 
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Parametric Analysis 
K y influence 

The following plots show the variation of the first and second natural frequencies as 
k y varies. In each plot the components of the eigenvector is also plotted. 



ky 



Figure 48 Dependence of first natural frequency on k y 
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Figure 49 Dependence of second natural frequency on k y 
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In the previous plots we can notice that the natural frequencies increase with k y , 
having as asymptotes 2.99 and 3.29 Hz respectively. Increasing the value of k y therefore 
contributes to separating the two frequencies up to the value of 0.3 Hz. The following 
plot shows the influence of k v on the displacement of the constrained point. 


(y-0d) /y 



Figure 50 Dependence on k y of translation of constrained point in first mode 
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Figure 51 Dependence on k y of translation of constrained point in second mode 

As k y increases, the displacement of the point S (constrained point) decreases, so that 
the differential displacement in mode 1 is reduced and likewise for mode 2 (see Figs 50 
and 51). 
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I c influence 




Figure 52 Dependence of first natural frequency on l c 
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Figure 53 Dependence of second natural frequency on Ic 

The value of I c affects the frequency of the second mode; decreasing the value of I c 
the natural frequency increases, yielding also a larger separation between the two first 
frequencies. I c has no influence on the displacement of the pivot point. 
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m c influence 

The effect of m c on the frequencies considered is to increase the separation of the 
frequencies as m c decreases, similarly to the effect of l c . The plot below shows this effect 
versus the m c /m A ratio. 


nat freq [Hz] 



Figure 54 Dependence of natural frequency separation on m c /m A 

Combining the effect of I c and m c , the separation of the frequencies increases from a 
value of 0.3 Hz when I c and m c are much greater than I A and m A , to a value of 2.17 Hz 
when the ratios I C /I A approaches unity . 

Relative variation off, m A , k yA , k lA 

We perturb the values of the parameters k yA , I A , k tA , d A , m A of the two sensing masses 
in order to test the influence of possible imperfections in the mechanical and physical 
characteristics of the sensing masses and springs. The ratios p A /p B , where “p” indicates 
one of the above mentioned parameters have been varied between ±5% of the nominal 
value. 

The main influence of these parameters is on the difference of displacement of the 
two sensing masses in the common mode (mode 2), resulting in a differential 
displacement. The plot below shows the differential displacement (y A -y B )/y A versus the 
p A /p B ratio. 
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Figure 55 Influence of bodies A and B parameter’s ratio on the differential displacement 

As it can be seen, the most influent parameter is the distance of the constrained point, 
and subsequently the mass, and inertia ratio. This kind of dependence is explained by the 
fact that d A changes the moment of inertia around the pivot axis with a squared 
dependence; m A influences both the translational and the rotational motion, while I A 
affects only the rotational degrees of freedom (DoF). 

Concluding Remarks 

The analysis performed allowed to describe the differential accelerometer’s modes 
and frequencies and their dependence on the inertial and mechanical parameters. The 
system has been modeled using 6 DoF, the rotation around the x axis and the translation 
along a radial axis for each body. 

The first two frequencies and eigenvectors have been analyzed in detail because these 
modes of rotation around the pivot axis affect the acceleration measurement the most. 
The separation of these two frequencies depends mostly on the mass and inertia moment 
of the C body, and slightly on the lateral stiffness of the springs. The lower C mass and 
inertia are, the greater the separation between the first two frequencies (if I*-.^ = 1 = 
mc/mA the difference is 2.17 Hz for a value of k y = 18000ji 2 , while it is only 0.3 Hz if 
I C /I A =110 and rric/mA = 6). Increasing k y contributes to reducing the displacement of the 
pivot axis both for mode 1 and mode 2. 

A variation of the parameters ratio between body A and B has also been checked, to 
evaluate the influence of possible imperfections of the inertial and mechanical 
characteristics of the proof masses. 
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The differential displacement of the CM of the sensing masses has been analyzed. 
The k lA /k lB ratio has the greatest influence, causing a Ay difference of 5% as the previous 
ratio varies by 0.7%. The k yA /k tB ratio has a similar influence. The ratio m A /m B has 
roughly half the influence of the spring stiffness ratio, causing a 5% Ay difference as the 
ratio varies by 1 .3%. The inertia ratio I A /I B produces a 5% variation of the Ay differential 
displacement as it varies by 2% from the unit value. Finally the lateral stiffness ratio 
k yA /k yB has a negligible influence compared to the others, i.e„ approximately two orders 
of magnitude less, reaching a 0.025% variation of Ay as the parameter ratio varies by 1%. 
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Dynamics Analysis of Instrument Package/Detector 


Introduction 

The following analysis is a study of the accelerometer and instrument package 
dynamics. Two different designs of sensor have been taken into consideration as follows: 
(1) a sensor with two sensing masses and a case (three bodies); and (2) a sensor with one 
sensing mass (made of two different materials) and a case. In both types the measurement 
is carried out by monitoring the differential displacements between masses, through the 
measurement of capacitance variations caused by those displacements. 

In this study we derive the differential equations of motion for each sensor and carry 
out numerical integrations of cases of interest, starting from different initial conditions of 
the positions of the centers of mass (CM) and initial rotational velocities of the bodies in 
order to evaluate the influence of those parameters on the differential output of the 
accelerometer. 

Each body is defined using all its six degrees of freedom (DoF) and elastic forces and 
torques are applied to each of them at the attachment points of the constraining springs. 

Description of sensor configurations 

Three-bodv sensor. 

The sensor is composed by two sensing masses having the CM coincident with the 
CM of the external case. One mass has a dumbbell shape, while the other is a hollow 
cylinder. These two masses are attached to the case by means of elastic springs, having a 
pivot axis that is parallel to the longitudinal axis of symmetry of the sensor. All the 
bodies spin about the symmetry axis x (during the fall). Figure 56 shows a schematic of 
the 3-body sensor. 



spin axis 


pivot axis 


Figure 56 Schematic of three-body sensor 
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The motion of the sensor is essentially a rotation of the sensing masses about the 
pivot axis. The translational motion of the sensing masses, which depends on the high 
lateral stiffness of the springs is much smaller than the rotational motion. A rotation of 
the mass implies a variation of the distance between the sensing mass and the external 
case that changes the capacitance of the motion pickup system. 

Two-bodv sensor. 

The two bodies sensor consists of an external case (body A) and an internal sensing 
mass. The internal mass is actually formed by two equal masses of different materials 
which are rigidly attached so as to behave like a single rigid body. The internal mass is 
connected to the external mass by means of elastic springs which have a lateral stiffness 
much higher than the torsional one. All the bodies spin about the symmetry axis x 
(during the fall). 

I 



A violation of the EP during the fall would cause a rotation about the pivot axis of the 
sensing mass with respect to the sensor case. The measurement is then carried out by 
detecting the differential displacement (and therefore the capacitance variation) of the 
central disk of body B with respect to the four fixed capacitor plates attached to body A 
(see Figure 57). The capacitor plates are arranged in such a way that a rotation of B 
about the pivot axis z unbalances the capacitive measurement bridge and produce an 
output signal. 

Reference frames 

A body reference frame has been attached to each body. The attitude of each body is 
identified by Euler’s angles, following the sequence 1-2-3 (i.e., by building the rotation 
matrix using sequential rotations about x, y, and z axes). 
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The relationships between the different reference frames have been written by means 
of rotation matrices and translation vectors. All the matrices are described in details in the 
following. The matrix transforms the inertial frame (X*,, Y 0 , Zq) into the body A 
reference frame (X A , Y A , Z A ); 


Z A 



Figure 58 Inertial and body reference frames 

The body frame is attached to the geometrical center of the body which, in ideal 
conditions, coincides with the CM of the body. The position of the CM (of body A for 
example) in the A reference frame is defined by {x CMA , y CMA , Zc MA }. The coordinates of 
the origin of the A frame (R A ) with respect to the inertial 0 frame are: {x A (t), y A (t), z A (t)}. 

The matrix yielding the transformation from system X„, Y,„ 7 ^ to X A , Y A , Z A is 
composed by the following sequence of transformations: 

- rotation 0 xA (t) around X axis; 

- rotation 0 >A (t) around Y’ axis (transformed Y 0 axis after 0 xA (t) rotation); 

- rotation ©^(t) around Z” axis (transformed Z,, axis after 0 xA (t) and 0 yA (t) rotations); 

- translation R A projected along the three axes X A , Y A , Z A . 

The rotation matrix [RoJ and the translation vector {R A } from system 0 (inertial) to 
system A (body A) are the following: 
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c(d yA )c(e M ) c(e M )j(e Xi4 )j( 0 vA ) + c(0 tA >5(0^ ) -c( 6 xA )c(0^ ) + s(d xA )s(e zA ) 

[ R(M J ~ — ^ )c(@zA^ ~ s(d xA )s(d y A)s(dzA ) ^'^:A ) "*" ^'^xA )^'(^v4 ) 

*(0m) -c(Oy A )s(O xA ) c(e xA )c(d yA ) 

where the superscript identifies the relevant reference frame. The translation vector is the 
R a vector, expressed in the A frame. This vector can be computed as follows: 

Ra = [RoaK-Xa -yA - Z a} T 

Where the vector {-x A , -y A , -z A } is the opposite of the vector {R A } expressed in the 0 
frame. 

We have: 

{-c(0 xA ) c(6. a )x a - [ 0 ( 6 ^ )s( 0 xA )s( 0 yA ) + c(0 xA )s(0 zA )]y A + [c(^ )c(^ )s(0 vA ) - s(0 xA )s(0 :A )]z A 
{/?/}-• c{e, A )sW :A ).x A - [c(0 tA )c(0^ A ) - s(0 xA )s(0 xA )s(0 :A )\y A - {<-(0^ )s(0 xA ) + f (0 xA )s(0 yA )s(0 zA )lz A 

-s(0 yA )x A + c(0 xA )s(0 xA )y A -c(0 xA k (0 vA )z A 

For the opposite transformation, which leads from the body A coordinate system to 
the inertial frame, we use the transpose of the rotation matrix [R A0 ] = [R 0 a] T an d the 
vector is simply the R A vector expressed in the inertial coordinates {x A , y A , z A }. 

We have 

c(e xA )c(d yA ) -c(d yA )s(d xA ) s(d yA ) 

[Rao]= cW^sie^ + dd^sie^sid^) c(0 xA )c(O^)-s(0 xA )s(O yA ).s(0 :A ) -c(6 vA )s(6 :A ) 

[-c(0 vA )C{0 :A )s{ e yA )+s(6 xA )s(d :A ) c(0 :A )s( 6 xA )s(d yA ) + c(d xA )s(d zA ) c(d yA )c(0, A ) 

yA 

ZA 

In summary, the transformations are as follows: 

Inertial to body A: 

[RoaMi- 0 } + {Ra A } 

Body to inertial: 

[RaoJ. {1^} + {R a °} 

The same procedure is followed for the other bodies. 
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Degrees of Freedom and equations of motion 

The transformation matrices have been written for each body for deriving the 
equation of motion for the different bodies. The body coordinate systems are centered in 
the body’s geometrical center, and fixed with the body itself. The position of the CM is 
described using three coordinates (x CMA , y CMA , Zc MA for body A). Each system has the x 
axis along the longitudinal axis of symmetry, corresponding also to the spin axis, the y 
and z axes are radial axes rotating with the body. 

Each body has all the six degrees of freedom, i.e., the three translations of the body 
frame origin, and the three rotations according to the Euler’s angles described above. 

To derive the equations of motion the translation accelerations of the CM and elastic 
forces have been written in the inertial system. The rotational equations are then 
projected onto the body axes. 

Translational accelerations 

In the following we describe the procedure for deriving the equations of motion by 
taking the two-body sensor as an example. The procedure to derive the equations for the 
three-body sensor is exactly the same with the added complexity of a higher number of 
equations. 

The CM vector in each body’s reference frame is: 

Body A. {r A } = {x CMA Zcma} 

Body B: {r B } = {x*. 

MB YcMB ^CMb } 1 

The translational acceleration expressed in the inertial reference is obtained using the 
formula: 

} = [^Tio] ■ { r 4 } + 2 [^4 o] ' } + [ R Ao ] ■ {*A } + {^4 } 

Where 

{a A } is the acceleration of body A in the inertial frame; 

[R a0 ] is the rotation matrix from A to inertial frame; 

{r A } is the coordinates vector of body A CM in the A frame; 

{R a 0 } is the translation vector expressed in inertial coordinates 

( ) indicates the time derivative and the dot (.) matrix multiplication. 

The same equations are used for the other bodies. 
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Elastic Forces 


To evaluate the elastic forces we define first the points where the springs are attached 
to each body. The position of these points are expressed in body frame. The distance 
between connected points is then computed from the actual position during the motion, 
projected in the A frame, and multiplied by the stiffness vector {k xA , k yA , k zA }. By 
following this procedure, it is possible to assign different stiffness to each degree of 
freedom. Finally the elastic forces are projected onto the inertial coordinate system, and 
added to the equation of motion. The points connected by the springs are as follows (all 
distances are in meters): 


{p A1 } = {0,0,-75xl0-y 
(Pa2> = (0,0,75x10‘ 3 } t 
{Pb,} = {0,0,-65x10- 3 } t 
{ P b2 } = {0,0,65 x10- 3 } t 

The initial lengths of the two springs constraining the motion arc as follows: 
{1.} = (Pa. - P Bi } = <0,0, -10x10 3 } T 


{1 2 } = {Pa:-Pb2} = {0,0,10x10- 3 } t 

The figure below shows the position of the points where the forces are applied 



Figure 59 Locations of attachment points 
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{F ia } = [R ao ] • ([K] . ({p A1 } - [T ()A ].{T B0 .p BI > - {!,})) 
{F 2a }= [R ao ] • ([K] . ({p A2 } - [T 0A ].{T B0 .p B2 } - {1 2 })) 
{F A }={F IA } + {F 2A } 


In the above notation the operator “T 0A .x” indicates the combination of the 
multiplication by the rotation matrix [Ro A ] applied to the vector {x} plus the translation of 
the vector {R A A } which yields: 

{T oa .x} = [Ro A ].{x} + {R a a } 


The expression T^XT^-Pb,) projects the coordinates of the point p BI (expressed in the 
body B frame) onto the body A coordinate system. Subsequently, the force is projected 
onto the inertial frame using the [R A0 ] matrix. A similar procedure is used to evaluate the 
forces acting on any other body pairs. 

The equations of translational motion for the two bodies A and B are: 


m A {a A }-{F A } = 0 


m B {a B } - {F b } = 0 


where m A and m B are the masses of body A and B, a A and are the accelerations, and F A 
and F b are the elastic forces: 


{F a } = {F ia } + {F 2A } 


{F b } = {F 1B } + {F 2B }. 


Rotational accelerations: 

The angular velocity of each body is computed by using the rotation matrices that 
transform the coordinate system from the inertial to the body frame and conversely for 
the opposite transformation. The rotational velocity matrices of bodies A and B are 
derived by using the Cartan’s formula as follows: 

[ w a] = [^o/i]*[^o] 

[ <w b] = [^Ob]'[^Bo] 

where [R« A ] and [R A0 ] are the rotation matrices from coordinate system 0 to A and vice- 
versa. [Ro B ] and [R^J are the correspondent matrices for body B. 

It should be noted that the former expression lead to the skew symmetric matrix of 
the angular velocity from which the components of the angular velocity vector {u> A } can 
be readily extracted. 


102 



The rotational inertia forces are obtained using the formulas: 

Va ] ■ {«U } + [^A ] ' Va ] ■ {<°A } 

[ l B ] ‘ { W B } + \. m B ] ' [h ] ‘ Wb } 

Where {co A } and [oo A ] represent the angular velocity vector and matrix, respectively, 
for body A (and similarly for body B), [1 A ] and [I B ] are the inertia matrices for body A 
and B which, assuming principal axes, have the form 



I A.xx 

0 

0 ‘ 


Ib xx 

0 

0 ' 


0 

l Ayy 

0 

; M- 

0 

hyy 

0 


0 

0 

*Azz. 


0 

0 

hzz. 


Elastic torques 

The elastic torques acting on each body are computed by using the definition of the 
attachment points and the expression of the elastic forces previously defined. The 
expressions of the torques in the body reference frame are as follows: 

{T a } = {p A j } x ([Ro A ] • {F 1A }) + {p A ^} x ([R 0A ] - {F 2A }) 

{T b } = {p B i } X ([RobJ -{Fib}) + {Pb?} X ([Rob] • {F 2 b}) 

where x indicates the external product of vectors. 

Finally, the expressions of the rotational equations of motion are: 

t^]' W} + ["/t]‘[ / /l]'{ W A}-{ 7 A} = 0 

[/ B ]- {(0 fi } + [cug] - [fg]' ~{Tb) = 0 

The equations of motion for the three bodies sensor are evaluated at the same way as 
for the two bodies. The same quantities are derived and the same structure for the 
equations is obtained. Figure 60 shows the positions of the constrained points for the 
three bodies sensor. 

The points p CA1 , p CA2 , p CB1 and p CB2 are the points of the external case C connected 
with the points p A1 , p A2 , p B1 , p B2 , respectively. In this case four springs are used to connect 
the two bodies to the case (two springs per each body). 
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Figure 60 Locations of attachment points 

The coordinates of the points are the following (each point is expressed in its own 
body system): 

(PaJ = {-40xl0 \ -60x10 \ 0} T 
{pBi } = {-250x10 \ -60x10 \0} T 
{PcAi> = {-50xl0'\ -60xl0'\ 0} T 
{PcBi> = {-270x10 \ -60xl0'\0} T 

{^Al } = (PaI ~ PcAl } 

0b|} = (PBI — PcBl } 

(PaJ = {40x1 0' 3 , -60xl0'\0} T 
{p BI } = {250x1 O' 3 , -60xl0‘\ 0} T 
{Pcai } = {50x1 0' 3 , -60xl0'\0} T 
{Pc B1 } = {270x 10\ -60x1 0' 3 , 0} T 

{{aj} = {Pa2 — Pca:) 

{Ibj) = {Pb: — Pcb’} 
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Two Bodies 


The integration of the equations of motion has been performed using different values for 
the CM position and initial angular velocity of body A. 

A frequency analysis has been performed for each result, to evaluate the frequency 
content of the differential rotation around the z axis, which leads to the detection of a 
differential signal by the capacitors. 

The parameters in the equations are the following: 

m A = mass of body A 

m B = mass of body B 

I Ax = moment of inertia of body A about the x axis 

I A = moment of inertia of body A about the y axis 

I Az = moment of inertia of body A about the z axis 

I Bx = moment of inertia of body B about the x axis 

l By = moment of inertia of body B about the y axis 

1 B/ = moment of inertia of body B about the z axis 

k x = lateral stiffness of spring connecting bodies A and B (in A coordinates) 

k y = lateral stiffness of spring connecting bodies A and B (in A coordinates) 

k z = lateral stiffness of spring connecting bodies A and B (in A coordinates) 

k Bz = torsional stiffness of spring connecting bodies A and B 

x C ma = x position of CM of body A (in body A coordinates) 

y CMA = y position of CM of body A (in body A coordinates) 

z CMA = z position of CM of body A (in body A coordinates) 

x CMB = x position of CM of body B (in body B coordinates) 

Ycmb = y position of CM of body B (in body B coordinates) 

Zcmb = 2 position of CM of body B (in body B coordinates) 

The variables for the equations (DoF) are: 

x A (t) = x coordinate of the geometrical center of body A (in inertial coordinates) 

y A (t) = y coordinate of the geometrical center of body A (in inertial coordinates) 

z A (t) = z coordinate of the geometrical center of body A (in inertial coordinates) 

x B (t) = x coordinate of the geometrical center of body B (in inertial coordinates) 

y B (t) = y coordinate of the geometrical center of body B (in inertial coordinates) 

z B (t) = z coordinate of the geometrical center of body B (in inertial coordinates) 

0 xA (t) = x rotation of body A 
0 yA (t) = y rotation of body A 
0 zA (t) = z rotation of body A 
0 xB (t) = x rotation of body A 
0 yB (t) = y rotation of body A 
0 zB (t) = z rotation of body A 

Numerical values have been assigned to the parameters and then a study of the influence 
of the position of the CM with respect to the geometrical center of each body (CM 
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parameters), also combined with an initial angular velocity of body B. The values used 
for the fixed parameters are in the following. 

We analyze in the following representative cases with prograde precession (i.e., ICx > 
ICy). The reason for favoring prograde precession over retrograde precession is that 
prograde precession brings about significantly smaller accelerations due to imperfect 
release (see later on) than retrograde precession. 

Two-Body Dynamics Simulation 

Prograde precession 
IAx < IAy 

mA = 30 kg, mB = 1 kg, 

IAx = 95/100 kgm 2 , IAy = 146/100 kgm 2 , IAz = 146/100 kgm 2 , 

IBx = 17/100, IBy = 17/100, IBz = 17/100, 
kx = 45000 N/m, ky -> 35000 N/m, kz -> 35000 N/m, 
kx2 = 45000 N/m, ky2 = 35000 N/m, kz2 = 35000 N/m, 
k0x = 1000 Nm/rad, k0y = 1000 Nm/rad, k0z = 54 Nm/rad; 

spin frequency = 0.3 Hz 

initial velocity 0 Z of body A = 10 2 rad/s 

initial nutation angle = 0 


xCMA = 0, yCMA = 0, zCMA = 0, 
xCMB = 10 6 m, yCMB = 0, zCMB = 10 6 m 


/ 2. 99716 i 
7.38688 
8.08957 
42. 8045 
42. 8045 
1 48. 5357 j 

Figure 61 Natural frequencies 
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Figure 67 FFT of 0 zB - 0, A vs. frequency (Hz) 



Figure 68 Zoom of 0 zB - 0 zA FFT vs. frequency (Hz) 
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Three-Body Dynamics Simulations 

The parameters in the equations are as follows: 

m A = mass of body A 

m B = mass of body B 

m c = mass of body C 

I Ax = moment of inertia of body A about the x axis 

I Ay = moment of inertia of body A about the y axis 

l Az = moment of inertia of body A about the z axis 

I B x = moment of inertia of body B about the x axis 

I B> = moment of inertia of body B about the y axis 

I Bz = moment of inertia of body B about the z axis 

l Cx = moment of inertia of body C about the x axis 

I Cy = moment of inertia of body C about the y axis 

I Cz = moment of inertia of body C about the z axis 

k Ax = lateral stiffness of spring connecting bodies A and C (in C coordinates) 

k Ay = lateral stiffness of spring connecting bodies A and C (in C coordinates) 

k Az = lateral stiffness of spring connecting bodies A and C (in C coordinates) 

k Bx = lateral stiffness of spring connecting bodies B and C (in C coordinates) 

k By = lateral stiffness of spring connecting bodies B and C (in C coordinates) 

k Bz = lateral stiffness of spring connecting bodies B and C (in C coordinates) 

k BxA = torsional stiffness of spring connecting bodies A and C 
k„ xB = torsional stiffness of spring connecting bodies B and C 
x cma = x position of CM of body A (in body A coordinates) 

Ycma = y position of CM of body A (in body A coordinates) 

Zc MA = z position of CM of body A (in body A coordinates) 

x cmb = x position of CM of body B (in body B coordinates) 

ycMB — y position of CM of body B (in body B coordinates) 

Zf-MB z position of CM of body B (in body B coordinates) 
x cmc x position of CM of body C (in body C coordinates) 
y CM c y position of CM of body C (in body C coordinates) 

Z( Mt z position of CM of body C (in body C coordinates) 

The variables for the equations (dof) are: 

x A (t) = x coordinate of the geometrical center of body A (in inertial coordinates) 
y A (t) = y coordinate of the geometrical center of body A (in inertial coordinates) 
z A (t) = z coordinate of the geometrical center of body A (in inertial coordinates) 
x B (t) = x coordinate of the geometrical center of body B (in inertial coordinates) 
y B (t) = y coordinate of the geometrical center of body B (in inertial coordinates) 
z B (t) = z coordinate of the geometrical center of body B (in inertial coordinates) 
x c (t) = x coordinate of the geometrical center of body C (in inertial coordinates) 
y c (t) = y coordinate of the geometrical center of body C (in inertial coordinates) 
Zc(t) = z coordinate of the geometrical center of body C (in inertial coordinates) 
0 xA (t) = x rotation of body A 
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0 yA (t) = y rotation of body A 
0 zA (t) = z rotation of body A 
0 xB (t) = x rotation of body A 
0 yB (t) = y rotation of body A 
O^ft) = z rotation of body A 
0 zC (t) = z rotation of body C 
0 xC (t) = x rotation of body C 
©^(t) = y rotation of body C 


Prograde precession 

ICx < ICy 
spin freq = 0.3 Hz 


values: 

mA = 1 kg, mB = 1 kg, mC = 30 kg. 

IAx = 17/100 kgm 2 , IAy = 17/100 kgm 2 , IAz = 17/100 kgm 2 , 

IBx = 17/100 kgm 2 , IBy = 17/100 kgm 2 , IBz = 17/100 kgm 2 , 

ICx = 95/100 kgm 2 , ICy = 146/100 kgm 2 , ICz = 146/100 kgm 2 , 
kAlx = 45000 N/m, kAly = 35000 N/m, kAlz = 35000 N/m, 
kA2x = 45000 N/m, kA2y = 35000 N/m, kA2z = 35000 N/m, 
kAqx = 61 .68 Nm/rad, kAqy = 1000 Nm/rad, kAqz = 1000 Nm/rad, 
kB 1 x = 45000 N/m, kB 1 y = 35000 N/m, kB 1 z = 35000 N/m, 
kB2x = 45000 N/m, kB2y = 35000 N/m, kB2z = 35000 N/m, 
kB0x = 6168/100 Nm/rad, kB0y = 1000 Nm/rad, kB0z = 1000 Nm/rad 

initial 0 Z of body C = 10 2 rad/s 
Initial nutation angle = 0 

xCMA = 10 6 m, yCMA = 0, zCMA = 10 6 m, 
xCMB = 0, yCMB = 0, zCMB = 0, 
xCMC = 0, yCMC = 0, zCMC = 0 
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M2. 1084 
43. 4894 
4. 39904 
27. 2175 
2. 99982 
42. 5542 
3. 48568 
44. 0754 
48. 344 
50. 0619 
4. 34997 
. 26. 7789 

Figure 69 Natural frequencies 

The simulations have been run for 100 s (which is higher than the actual value) for 
the sole reason of providing a good frequency resolution for the FFT plots. 



Figure 70 z A - z c (m) projected in body C frame vs. time (s) 
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Concluding Remarks 

The analysis carried out allowed to evaluate the influence of CM position and initial 
angular velocity errors in the dynamic of the 2-body and the 3-body sensors. The two 
systems have been modeled using 12 DoF, and 18 DoF respectively (six dof for each 
body). 

2- body results 

The position of the CM for body A (external case) did not affect the rotation about the z 
axis (measurement axis), while the initial non zero condition on the angular velocity of 
the same body caused a differential signal to be detected. 

Analyzing the frequency content of the signal, though, no frequency peak was noticed at 
the spin frequency (0.3 Hz), while the natural frequency for the rotation about the z axis 
(3 Hz) and the precession frequency (0.18 Hz) were present. 

These frequencies can easily be changed modifying the torsional stiffness about z axis 
and changing the inertia ratio (IAz — 1 Ax)/1Az, to which the precession frequency is 
directly related. 

Note that the spin frequency appears in the single 0, A or 0 zB starting from an initial 
angular velocity of body A about y or z axes, but it is rejected when the difference 
between the two signals is taken. 

Applying all the initial conditions together (both CM positions and angular velocities) led 
to analogous results, resulting in a rejection of the spin frequency in the differential 
rotation. 

3- body results 

Likewise, in the three-body detector the position of the CM of the external case (body C) 
along the spin axis does not influence any of the variables (in the absence of gravity 
gradients). The displacement of the CM along the axes y and z cause the excitation of the 
y and z variables only, while the rotations are not affected. Moreover the differential 
displacement is not excited. 

An initial condition on the angular velocity (both along the y or z directions) excites all 
the dof of the system, both individual and differential; the frequency pattern of the 
variables doesn’t show evidence of peaks at the spin frequency, neither in the individual, 
nor in the differential component of the displacements. 

After applying all the different conditions together, the frequency analysis again shows 
that the differential signal (along z direction) doesn’t show any peak at the spin 
frequency, even if these peaks are present in the individual variables (y and z) due to 
initial displacement of the CM of body C along y and/or z directions. 
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Detector Requirements Development 
Model for the Gravitational Perturbations Acting on Proof Masses 

Overview 

The test body (proof mass), falling inside the capsule, is subjected to non-negligible 
gravitational attraction by the capsule. In particular, the higher order gravitational 
potential plays a major role because of the elimination of the zero-order potential due to 
the equivalence principle. 

The model consists of a spinning test body inside a capsule: hollow cylinder covered with 
flat caps. The test body is released at the axis of the cylinder, and deviates from the axis 
during its fall. Therefore, our task is to compute the gravitational force and torque acting 
at the neighborhood of the fall. The fact that both the test mass and the capsule are closed 
finite bodies, increases the complexity of the problem. We took advantage of the size of 
the test mass; smaller than the cylinder radius, to derive an asymptotic analytical solution. 
We are mostly concerned with the force/torquc at the modulation frequency. In other 
words, our model should evaluate the forcc/torque in a rotating body frame, rather than 
the static field in the capsule frame of reference. For that purpose we built a semi- 
analytical model. The main advantage of our model is its flexibility. It can handle any 
configuration of test mass as well as any additional mass distribution in the test chamber. 

The purpose of this work is two folds. We need a working gravitational model for future 
simulation for the dynamics of the test body. We also need to choose the inertia 
properties of the test body that minimize the gravitational disturbance. 

In the following we presents the analytical approach, analytical results, numerical 
analyses, and the interpretation of results. Appendices to this section are available in the 
first Annual Report #2 xxxii for this grant. Appendix A explains the computer code. 
Appendix B and C present analytical solutions for the gravitational attraction due to the 
cylinder and the caps, and check the validity of our numerical model. Appendix D 
proposes an approximation of the gravitational field and its gradient by using radial basis 
functions. 

General Approach 

The purpose of this work is to model and analyze the gravitational attraction between 
the capsule and the test body. There are at least three ways to approach this problem. The 
first approach is to compute the force between each capsule mass element and a test body 
mass element, and to perform a double summation on these forces. If N is the number of 
capsule mass elements, and N B is the number of test body mass elements, then the 
cardinality of the computation is N ® N B (the complexity of the computation is the 
product of the complexity of a single element computation by the cardinality). The 
second approach is a double integration over the bodies. The drawback of the first 
approach is the heavy computational effort, especially if we need, in the future, an online 
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computation in a dynamical process. The weakness of the second approach is in the 
cumbersome analytical computation. Moreover, the integration will result in an 
asymptotic series of complicated functions. This will require intense computation. In 
addition, asymptotic series introduce a truncation error into the computation. We adopted 
a third approach for this model. 

Our computational approach is to consider the capsule as a discrete ensemble of 
lumped mass, where the resultant force and torque acting on the test body are the 
cumulative force and torque due to each capsule mass element while the test body is 
terated with a series expansion. In other words, the interaction is between a finite body 
and a point mass. This approach is a tradeoff between the previous approaches, its 
cardinality is N . The main advantage is the flexibility in modeling any capsule shape, 
according to future requirements. 

The purpose of the following computation is to come up with simple closed-form 
analytical expressions, for the force and torque acting on the spinning test mass due to the 
capsule gravitational attraction. 

Gravitational Model 

The gravitational potential for finite size bodies is: 
d M 

V - -G / J -dM 

m m r ° 

o B 

Where in our case is the test mass (proof mass), and 9l 0 represents the attracting 
bodies, like the capsule (cylinder and caps) and the Earth. 

For the purpose of simplicity we will proceed with a representation of the test body, 
and a single element of the attracting mass, M t . 

The gravitational potential at a representative element mass M, due to the test body 
dM B 

is: V ( M ) = -GM J , where r is the radius vector between an element mass 

M b r 

of the test body and dM 0 as shown in Figure 78. 
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Figure 78 Gravitation model for test mass and outer attracting bodies 

Assuming a sufficient discretization, the total force acting on the test body due to all 

N 

outer mass elements is: F B s J VV(M,) , where N is the number of outer mass elements. 

f-l 

The model utilizes two coordinate systems. The first is the capsule frame, denoted by: 
{ X,Y,Z } • X, Fare so far arbitrary, while Z coincides with the symmetry axis of the 
main cylinder. The test body frame, denoted by{jc,j,z} , is attached to the test mass. So 
far, the origin and the orientation of the coordinate system are arbitrary. The gravitational 
potential will be expressed in body frame. The inertia coefficients of the body are 
constants in this frame. 


The analytical procedure is to express r as r = R- p 


1 

the following power series: — 



, and to expand the potential by 


Here P n is the Legendre polynomial of degree n. 
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The derivation consists of the following steps. Substituting the inverse radius 
approximation into the potential, carrying out the integration, while expressing vectors in 
body frame, as follows. R = (le x + me y + ne,)R and p=* xe x + ye y + ze z , where 

{e x e Y e.} are body unit directions, and /, m. n are direction cosines between R and the 
body axes. The resulting potential is a summation over all mass elements M,: 


„ % gm x m 2 gm x m 2t ,- _ _ n 

V = £{ - -\lx + my + nz\ 


y[(3/ 2 - l)Jxx + (3m 2 -l)Jyj +(3 n 2 - 1)J^ +6(/mJ Jl y + In J xz + mniy Z )] 

2 R 

~ [1(51 — 3) J XO + fTl (5 Ml — 3) Jyyry + W(5/J^ “ 3) 

+ 3m (5/ 2 - 1) Jjcxy + 3n(5/ 2 - 1) + 3l(5m 2 - 1) J^, v 

+ 3n(5m 2 -1) + 3/(5/? 2 - l)J X zz + 3m(5n 2 -1)J V2Z 

+ 6/mnJ^] + 0(-^-)> 

R 


where j xPy q. r are the inertia integrals defined as: 


Lp v q,r = Jx p y q z r dm 
' w m b 

For the second order, it is common to use the moment of inertia. The following 
relations relate them to the inertia integrals: 


J xx 2 ^vy + ^zz 

Jyy-\( I «+ I zz- I ») 

J xx ~ + — ) 

Jxy-Ixy 

J xz = I x: 

Jyz= Ax 
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We define <5 = {x,y,z} as the offset between the origin of the frame and the center of 
mass. It should be emphasized at this point that the origin does not necessarily coincide 
with the center of mass. Although the test bodies are typically axisymmetric, and ideally 
the geometrical center is the center of mass, mass and/or mounting imperfections can 
displace the center of mass from the geometrical center. 

Note that if the origin is at the center of mass, the body frame coincides with the 
principal axes, and 0(1 / R 4 ) is neglected, then the resulting potential is the so-called 
MacCullagh formula: 


V = - 


R 2 R 3 X K ’ 


where 1^ I = P'lxx + tn^lyy + /i^I^ is the projection of the second order 
inertia tensor i on R . 

The force acting on the body is the gradient of the potential (the negative gradient is 
the force acting on M,). 

The resulting total force components, in body frame, arc: 


F r = G 



■ l).v + 3/ ( m- v + 3lj tij z 


+ l(|r ; . - I VT - U + 21 „)lj + 3m,.(5/ ( 2 - l)I xy + 3« ( .(5/ ; 2 - l)I r; + 1 51^.1^] 

R i 

+ ^5 1 ^( 35 /, 4 - 30 if + 3 ) + Ulm} - 3 )/. m . J yyy + 1 ( 7 n? - 3 )/• n. 

+ T (7/ / 2 - 3 > + f (7/ / 2 " 3)/ / + |(35/, 2 m, 2 -Sir -5 mf + l)J XV v 

+ jOn,? ~ 1) // n i iyy Z + 1(35 l?n} - 5/ 2 - 5n 2 + l)S x:: + y(7/i, 2 - 1) 

+ 3<7/ ( - 2 - \)m-n- J X j 2 ] + H OT} 
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F > 


N M R m ; Mu F o i 

“ G . 1 + TT [ 3 l i m i 7 + (3m / 2 “ ] >- v + 3 m i n i ? ] 

i = l Rj /?, J 

+ ~4 [( f r / ~ ‘ xx ~ l zz + 21 yy ,m , + 3/, (5m, 2 - Dixy + 1 5/,-m,-n,- I xz + 3n,.(5m, 2 - l)I yz ] 

+ (7/ / 2 ~ 3) l i m i J xxx + ~-(35m ; 4 - 30m ( 2 + 3) J yyy + |(7w,- 2 - 3)m,n,- J zzz 

+ ^(35 /,- m,- -5m,- -5/,- + l)J x *y +y(7/,- -llm,-/!,-.!^- +y(7m, 2 -3)/,-m,-J xyy 
+ —(7m,- - 3)m,m, Jyy Z + —(In 3 - l)(,-m,-J jrz , + — (35m, - 2 w,- 2 -5m,- 2 -5w, 2 + 1) J yzz 
+ 3(7m,- 2 — 1 ) / ( - J X y Z ] + HOT} 


N vf / jtf „ r ^ i 

F. = G 2 Mj { — — + — 3/. n- x + 3wi /i . v + 3(/z- - 1)- 
i = l /?/ fl/L ' ' 1 ' 1 J 


+ 4 fxv Iy V + 2I„)/i,- + 15/,-m,-n,-l X y + 3/,-(5n,- — l)I t . +3m-(5n- — 1)I V . ] 

+ yj I [|( ( ?/ / 2 - 3)/,- 1 (7m, 2 - 3)m,- n,- J vv> . + ^(35w,- 4 - 30n,- 2 + 3)J =: 


+ ^(7/ 2 - Dmrn,J m . + |(35/, 2 n, 2 - 5/, 2 - 5 n, 2 + l)J„ ; + y (7m, 2 - 1)/,- «,- J. n . v 
+ |(35m, 2 „, 2 - 5m, 2 - 5«, 2 + 1) + y (7«, 2 - 1)/,- + y(7„, 2 - 3)m,.„,.J VK 
+ 3(7n, 2 - 1) Ijirtj J rv , ] + H OJ} 

where 


+ ^Ij cx)/,- +(Jxx + ^zz + (I xr + lyy /1,- 

Note that the force has the following order of magnitude: 

— 2 ^ 

F « ZG^2-{0( 1) + 0(-£-) + O&r) + Oi^r) + HOT} 

i Rr R i R t 2 /?, 

Where L g is a typical length of the test body, — « 1 , and — « 1 . 
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Gravity Torques 


The torque, acting on the test body, can be computed with respect to the center of 
mass or with respect to the geometrical center. Both coincide when the body has perfect 
symmetries. In reality there are mass imperfections that cause the center of mass to 
deviate from the geometrical center. Since we have no a-priori knowledge of the 
imperfections, our reference point for computing the torque is the geometrical center. 



Figure 79 A Model for the Torque 

There are few approaches for computing the gravity torque. Our approach is to 
integrate the torque due to mass element in each body. The torque with respect to an 
arbitrary point is, therefore: 

T 0 = Tc.m + dxF B 
where. 


T CM -GJ fdx^dM B dM Q * 

B r 


dxR 

G inij f — — dM B 

i=\ di B r 
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F B =G / / ^dM B dM 0 = c|m, J 

,= I SR s r 3 

There are three approaches to compute the integrals. The choice of the approach 
depends on the size of each body. If both and 9t 0 are small compared with R, than 
we can expand r, and obtain an analytical approximation. If only 9ig is small, as in our 
case, we can approximate 9^ but have to integrate numerically in 9t 0 . If both are large 
relative to r, then the only way is to perform a double numerical integration. In our case, 
not only 3t 0 is large, it is also surrounded by so there is no topological way to 
expand r in 9t 0 . 

N j 

T =3 G 2 AT. { [m n (I -I ) + / n I -l m 1 + (/? - rn )/ ] 

x / = | 7 ^3 / / zz yy i i xy i i xz /' i yz 

i 

1 2 2 2 

+ — — J[n (5m - \)J -m(Sn -\)J + n (5/ - \)J +10 ImnJ 

2 R. i i yyy i i zzz i i xxy i i i xyy 

- m (5/2 - \)J -10 ImnJ + m (lo7 - 5m* + \)J 

i i xxz i 7 7 XZZ i I 7 yyz 

-n (10/7 -5n 2 +\)J +10 l(n 2 -m 2 )J ] + H.O.T} 

i i i yzz i i /' xyz 

N 1 9 9 

T * 3G J MJ [/ w (/ - / ) - m n I + l m I +(/ - m )I ] 

J i = 1 f ' / i xv zz / / xj ? / / yz i i xz 

i 

+ [-« (5/ 2 - 1)7 + / (5« 2 - 1)7 - 10/ m n J - n ( 5m 2 - 1)7 

2R i i xxx i i zzz i i i xr v / / xyy 

- / (10 n 2 - 5/ 2 + 1)7 + a? (1 0/ 2 - 5 w 2 +1)7 + / (5m 2 - 1)7 + 10/ m n 7 

/ f * xxz i i i xzz i i yyz i i i yzz 

+ 10 m(l 2 -n 2 )J ] + H .O.T} 

i i i xyz 
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T 

z 


N \ 2 2 

= 3G 2 M { [/ m (I -I )- (m -/ )/ +m n I -l n I ] 

■ _ j i r 3 i i yy xx i i xy i i xz i i yz 


+ -J— [m (5 1 2 - 1 )J - / (5m 2 - 1 )J +1 (10m 2 - 5/ 2 + 1 )J 
2 R 4 i i xxx i i yyy i i i xxy 

i 

- m (10/ 2 - 5m 2 + 1 )J + 10/ m n J + m (5 « 2 - 1)J - 10/ m n J 

i / i xyy i i i xxz i i xzz i i i yyz 

-l{5n 2 -\)J +10« (m 2 -l 2 )J ] + H .O.T} 

i i yzz i i i xyz 


The order of magnitude of the torque is: 

T « + O(^) + HOT} 

i R, R, R, 

If the origin is at the center of mass and the body frame coincides with the principal 
axes, then the resulting torque is reduced the well-known expression (based on 
MacCullagh formula): 

N M 

7>3G(I zz -I yy )2-f/n,n,- 

7= 1 /C|* 

N m 

T y - 3G (I xx -I zz ) l-fhnt 

i=\Rj' 

N M. 

T z as 3G (I yy - I xx ) h m, 

i=[Rj 


Proposal for frequency decomposition 

We would like to propose a more general approach for a closed-form solution of the 
force/torque in terms of their frequencies. 

Let us expand the force in potential orders as well as a Fourier series in the spin 
angle, assuming a pure spin around the axial axis of the test body. 

F = G$-^ 2 [a { „ d) cosn(p + b { „ d) sin ncp]}dM 0 

91 „R d = o «=i 


125 


where d is the degree of the inertia. Let J x i v j k = / x'y'z k dM B than d = i + j + k. 

Assuming that we have carried out the integration (analytically or by numerical means), 
F can be expressed as: 


F = G 2 {Ax) d) + ^ [A { n a> cosncp + Bl d) sin n(p]} 

d = 0 n = 1 


d + 1 


\(d). 




(d) 


where A <d) = / ^dM 0 , and B { f = f ^ dM a 


9ln R 


9tn R 


These coefficients may be viewed as a generalization of the inertia coefficients of 9t 0 . 
If 9f 0 is topologically connected (i.e., is outside9( 0 ), and if 9t 0 is sufficiently 
small, then one may asymptotically expand these integrals. This will result in the inertia 
coefficients of9i 0 . In our case, these coefficients can be computed numerically, or 
sometimes analytically, when is simple enough. 

For a more general rotation, one needs to express the direction cosines in terms of the 
general transformation matrix (rather than the single-axis rotation). 


Order of Magnitude and Similarity Analysis 

The dominant sources of the gravitational attraction are the capsule, the Earth, and to 
a lesser degree the Moon. Let us perform a qualitative comparison of each source. First, 
we should note that the force acting on the test mass is not a homogeneous function of the 
mass and the radius (as a force acting on a point mass). This is because the parameters 
R and M g contribute as: 


Fee 



or 


f«f (0) + f (1) + f< 2 > + f <3) . 


where F a)) represents the 2-Body term, and so on (the order enumeration is according to 
the corresponding Inertia order). Therefore, it is impossible to express the ratio of forces 

from different sources as F(Ml,R])/ F(M2M2)cef((M\/ Ml) m ,(Rl/Rl) n ). 

On the other hand, each order is homogeneous, that is 
F (k \Ml,Rl)/ F (k \M2,R2 ) = (Ml/ M2)(R2l R\) (k) 

Therefore, in order to compare different sources of attraction, we must consider each 
order separately. 
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The following qualitative discussion concerns the force difference between the two 
test bodies. Let us denote /() as a generic function, then a generic expression for the 
force acting on test mass B } is: 

f b ‘ - + APy;>o(pr) - + /< 

We present two models for the force difference. In the first model the two test bodies 
are centered, and the inertia components of the bodies are different. In the second model, 
the bodies have equal inertia, but are not centered. These two models are idealization of 
the two sources of difference, while the reality is probably a combination of these. 

The force difference for the first model is: 

AF = jP e 2 _ f b i = / ( AM;)o|-^-j + + /(A/ 2 ;)o|-^-j + 

The force difference for the second model is represented by a differential. For central 
deviation e , the force differential is: 

dF*F B ^ -F B ‘ a VF-e=/(M;)^j + /(p;)^j + /(/ 2 ;)^j + /(/ 3 ;)^j 

The ratio between force terms of order k from two different sources 1 and 2, is: 

F a) (M t ,R t ) _ M, //? 2 V 
F {k \M 2 ,R 2 ) 

The following table compares the force order of magnitude due to different sources. 
Note that we eliminated order 0 (2-Body), since the test body is in free fall. 


Table 1 Ratio of force terms for few orders for the dominant gravitational sources 


Inertia order 

Capsule / Earth 

Capsule / Moon 

Earth / Moon 

1 

4*10' 3 

8* 10 4 

2*10 7 

2 

3*10 4 

3 * 1 0 13 

1*10 9 

3 

2*10" 

1*10” 

6*10’° 

4 

1*10 18 

4*10 30 

4*10 12 
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It is obvious that the effect of the Moon is negligible. The Earth is dominant only 
with respect to the gravity gradient. Again, as the order under consideration becomes 
higher, the effect of the capsule increases. 

To gain more insight into the nature of gravitational attraction let us examine the 
force. It can be rewritten in the following form: 




Where c represents the direction cosines. The notations c° and c E show the nature of 
the degree of the direction cosines (odd or even). For example, a square of a particular 
direction cosine has an even degree, while a product of the three direction cosines has an 
odd degree. Now, suppose that the test mass is at the center of the cylinder. Since for 
each direction cosine to a mass element there is an opposite direction cosine, all the odd 
terms are cancelled out. Therefore, the only non-zero contribution is due to the even 
terms. 


Another important issue is the role of the test body mass and size. For this purpose we 
perform a similarity analysis, as shown in the following. 

We formulated the gravitational potential as an asymptotic series in MR n . Each term 
in the series has an inertial tensor of degree n - 1. For example, the first term is of order 
n = 1 and its inertia tensor is of order 0 (the test body mass). The third term is of degree 
n = 3 and its inertia tensor is of order 2. This non-uniformity introduces a problem, as 
mentioned earlier, when one computes the forces for a particular test mass, and wants to 
deduct the forces for a different scale of test mass. Mathematically speaking, if F(Mg^) 

and F(M bj) are the forces due to two different test masses, then the ratio between the 

forces is not a homogeneous function, that is, F(M ai ) / F(M B2 ) * g((M Bi I M B 2 ) m ) 
where g is a function and m is the degree of the homogeneity. Note however that each 

term of the series is homogeneous by itself. For example, the first term is homogeneous 
of degree one, that is, the ratio between forces equals the ratio between the masses. 

We will distinguish between two cases. In both, the mass distribution of the different 
test masses can be scaled. The first case is when the difference in masses is due to a 
different density. Since the forces are homogeneous of degree one in density, the ratio 
between the forces is equal to the ratio between the masses. In the second case the density 
is the same, and the different mass is due to different sizes. In this case we will apply a 
different similarity to each term. Let L and M be scaling factors for the length and the 

mass, respectively. If the density of different test bodies is the same, then L a M ' 3 . Let 

m - n - 1 be the order of the moment of inertia. Its similarity dimension is L m M , or 

Af I+m/3 , or L m+3 . The acceleration similarity is M m/3 or L m . The following table 
summarizes the similarity dimension for each term in the potential. 
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Table 2 Similarity relations for the gravitational potential terms 


Potential Order 
n : \IR n 

Inertia Order 
m = r 7-1 

Inertia 

Similarity 

Acceleration 

Similarity 

1 

0 

M (L 3 ) 

1 (1) 

2 

1 

M 473 (L 4 ) 

M ,/3 ( L ) 

3 

2 

M s/3 (L 5 ) 

M 213 (L 2 ) 

4 

3 

M 2 (L 6 ) 

M <Z 3 ) 

5 

4 

M 1 ' 3 (L 7 ) 

M 4/3 (I 4 ) 


To summarize, given the force (or the acceleration) due to a particular test body, we 
can deduct the force (acceleration) due to another geometrically scaled test mass. The 
acceleration ratio (of two different test masses) versus similarity dimensions (mass and 
length) is illustrated in the following plot. 




Figure 80 Scaling of Acceleration of Test Mass 

The maximum potential order considered in this analysis is the hexadecapole (n = 5) 
which involves the fourth-order inertia integrals. We computed numerically the forces 
associated with the hexadecapole for a body with equal second-order inertia integrals (in 
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order to minimize the quadrupole term). The results show that, at the worst expected 
location inside the capsule, the maximum acceleration associated with the hexadecapole 
term for a 1-kg proof mass are less than 10 16 g. In other words, for the accuracy goal of 
this experiment there is no need for belted cylinders (which would reduce the 
hexadecapole component even further). Proof masses with equal second-order inertia 
integrals (or alternatively moments of inertia) and construction accuracy of order ten 
microns are sufficient to make the contributions of all the higher-order gravity terms 
negligible. 


Effect of Inertia imperfection on a Spinning Test Mass 

Our goal is to minimize the gravitational forces acting on the test mass. The force 
equations reveal that the dominant term under our control is the term corresponding to the 
second order inertia. The CM offset effect may be bigger to start with but it can be 
reduced by centering. The direct way to minimize the forces is to require equal second 
moments of inertia. The residual forces then are due to mass imperfection, and to higher 
order inertia integrals. The purpose of the following discussion is to explore the effect of 
each term for a spinning test body. 

First we will evaluate the mixed-inertia terms due to the imperfection. 

We assume that the imperfection is due to disturbance in mass and in length. The 
nominal test body is a perfect cylinder, with radius R„ and length L B . In order to minimize 
the force we require that I xx - Wy = ^zz ■ This constraint results in a given aspect ratio of 
the cylinder as follows: L B = ^3R B . Given the mass density, the mass or the sizes of the 

cylinder are now functions of a single free parameter. For example, given the mass, M B , 

1 

( 3A/ 

— - 

jtp B 


The similarity dimension of the k-order inertia is: [/] = ML k . Thus, the perturbed k- 
ordered inertia is: [<5/J = L k 6M + kML k ~ l 6L , or 


<5/ 


dM 

+ k 

6L 

— 

= 





M . 


L 


Next we consider the effect of the spin. The test body is spinning about its x-axis with 
a frequency that will be regarded as a fundamental frequency, or PI (period one). We are 
primarily concerned with PI which is the period of the measured signal. For that purpose 
we will analyze/^ by substituting the direction cosines, shown below, into the force 
expression (note that y is the sensitive axis of the accelerometer). 
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/ = / (/; PI); m = m(t;P \ ); n = n(t;P\) 
i i i i i i 

It is important to note that the forces are expressed in body frame. The direction 
cosines are PI -periodic. Moreover, if the body fall such that the body x-axis is parallel to 
the capsule X-axis, only m and n are periodic, as will be seen from the following 
discussion. 

The periodicity in time is introduced through the direction cosines that represent the 
orientation of the body frame relative to the capsule frame. The transformation between 
the two is: 


x' 


1 

0 

0 ' 

X ' 

y 

= 

0 

cd 

sd 

Y 

z 


0 

-sd 

cd 

Z 


The resulting direction cosines are: 

/ = — = — * 1 ( 6 ) 

R R 

v Ycd + Zsd 

m = — 

R R 

z - Y sd + Z cd 

n ~~R R 

/ is a cyclic function only if the body deviates from the X-axis. If the test body is 
perfect, then only the point mass term and the first term of 0(1/ R A ) contribute to Fy ■ 
Otherwise, introduce higher harmonics. If the power of the direction cosines is even 
(N e ), than the additional harmonics are: P 0 ,P„---,P W£ . For odd power (N a ) the 
additional harmonics are: P t ,P } ,--,P No . The dominant time dependency of a perfect 

body is introduced through the first order attraction on a point mass. It results in a 
fundamental frequency, because the gravitational attraction is static in the capsule frame, 
while the measurement is in the rotating body frame. 

The expected frequencies are the following. PI from the point mass termO(l/P~ ). 
PI and P2 from the offset term 0(\/R ' ). Since F. contributes P2, the term of 0(\/R i ) 

contributes PI and P2. However, because of the almost axial-symmetricity of the test 
body, F. is almost a constant, so the dominant frequency of 0(l//? 4 ) should be PI . The 

term of 0{i/R' ) contributes PI , P2, P3 and P4, where P2 should be dominant. 

The analysis above is for the x-axis parallel to the X-axis, that is, for / = constant . 
Otherwise 1 = 1(8) an d higher frequencies are involved. The explicit orientation 
dependency of period one (PI) becomes: 
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F Body (6;P)= G l M 

y 1 ,=i 


J_f ( 2 ) +— f (3) + -Lf< 4 > 

5 2 * F 3 >’ F 4 >> 


+ J-F< 5 > 
f 5 y 


where 


F^ =M (Y. cos0+Z.sin0) 

y B i i 

F (3) =3M d ^ (f.cos0 + Z.sin0)J 

y B ‘ i i 

fl 4) = || [AX H. + 1(3 B + C)(T 3 + Y.z. 2 )] - BY. Jcos0 
+ j| [AX. 2 Z. + 1(35 + C)(Y. 2 Z. + Z 3 )] - BZ. Jsin0 
fJ 5) = 1(7 A'. 2 -3)(^cos0 + Z sine)^/^. 


and X f — X / R, Y : — Y R , Z, — Z! R are the direction cosines in the capsule frame. 

An interesting case is the effect of the attraction of a point mass located on the Z- 
axis. It may represent the Earth attraction on horizontally falling body, or the caps on the 
top or the button of the capsule. 

Let the attracting massM, be at Z1 = +5, in the plane X = Y = 0 . Using the 

relations: Z, = 5, sgn(Z, ), Z, =sgn(Z,), Z 2 = 1 , and substituting the attracting mass 
coordinates results in the following simple equation: 


7 Bod\ 


(0-,Pj,/ = O)= GMj 


Mr. 1 1 

2 + z^ (7B + 5C) 
Rf 8 R; 


Isgn(Zj)sin0 


Again, the relations above are the terms corresponding to PI. The other frequencies 
may be observed from analytical expansions, or from a numerical frequency analysis, as 
demonstrated in the next section. 
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Numerical Investigation of an Imperfect Spinning Test Body 

The scenario under investigation is shown in Figure 8 1 . The capsule consists of a 
hollow cylinder, covered by two flat caps and a total mass of 500 kg. The test body mass 
is 1 kg with sizes: R B = 0.0412m and L B =0.0713m. The geometrical and mass 
inhomogeneity errors are: 6L/L = 10 4 and 6M/M = 10 5 . Consequently, the offset 
components are: x = 1 .7 x 10 -6 m , y = 2 x 10 -6 m , z = 1 .7 x 10 -6 /n. The components of the 

inertial matrix are: I xx = 8.5 x 10 -4 kg-m 2 , lyy = 8.501 x 10 -4 kg-m 2 , I zz = 8.5015 x 10~ 4 
kg-m 2 . Note that the nominal values of those are the same. The mass perturbations 
prevent the elimination of the following inequalities: I xx + Iyy ~2/ Z z *0, 

Ixx + Itz~^yy * 0 , Iyy + hz~ 21 xx * 0 • 

I xy = 3 x 10 -8 kg-m 2 , I xz = 3 x 10 -8 kg-m 2 , Iyz = 4 x 10 -8 kg-m 2 . The components 
of the third-order inertia tensor are: I xxx = 5 x 10 -10 , /yyy = 8 x 10 _l °, I zzz = 8 x 10 -10 , 
/ xxy = 6 X 10 I XXZ = 6 X 10 I XYY = 7 X 10 ^ ilyYZ = 8 X 10 , I Y77 =6x10 , 

1 yzz = 8 x 1 0“ 1 0 , I XY z = 7 x 10- 10 (all in kg-m 1 ). 



1.2 [m] 


Figure 81 Coordinate systems for the test body and the capsule (dimensions account for 
the thickness of the capsule walls) 
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The gravitational force and the torque acting on the test mass were computed at 9 
points inside the capsule. Three points along the axial axis, three points along an axis 
deviated by 0.1 m from the axial axis, and three points along an axis deviated by 0.2 m 
from the axial axis. The ninth point is the farthest away, thus representing the worst case. 



-1 - 0.8 - 0.6 - 0.4 - 0.2 0 0.2 0.4 0.6 0.8 1 


Figure 82 Tested Points inside the Capsule 

We are mostly interested in F y ( F y is in phase with F z ) and in 7^ . Because of the 
free fall, we exclude the 2-Body terms. The following plots show the angular history of 
the force/torque for a single revolution with respect to the x-axis, and the corresponding 
spectra. The latter are normalized, such that the highest amplitude is equal to one. The 
results are given in the following 18 plots. The analysis of the frequency spectrum 
requires a closed look at the equations for the force and torque, considering the particular 
location of the test mass relative to the capsule. Some of the results are non intuitive. The 
general approach to estimate the results is to specify each term according to the degree of 
the direction cosines (odd or even), and consider possible cancellations due to anti- 
symmetricity with respect to the radius vector between the test mass and each capsule 
element. The results are shown in the following figures. 
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spin angle (rad) frequency ( HZ) 


Figure 83 Ordered Force and harmonics at Point 1 due to Capsule Attraction 



spin angle [rad] frequency [HZ] 


Figure 84 Ordered Torque and harmonics at Point 1 due to Capsule Attraction 
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Figure 85 Ordered Force and harmonics at Point 2 due to Capsule Attraction 



spin angle [rad] frequency ]HZ] 

Figure 86 Ordered Torque and harmonics at Point 2 due to Capsule Attraction 
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Figure 87 Ordered Force and harmonics at Point 3 due to Capsule Attraction 



spin angle [rad] frequency [HZ] 


Figure 88 Ordered Torque and harmonics at Point 3 due to Capsule Attraction 
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xIO 14 



spin angle (rad) frequency (HZ) 

Figure 89 Ordered Force and harmonics at Point 4 due to Capsule Attraction 



spin angle (rad) frequency (HZ] 


Figure 90 Ordered Torque and harmonics at Point 4 due to Capsule Attraction 
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spin angle (rad) frequency (HZ) 


Figure 91 Ordered Force and harmonics at Point 5 due to Capsule Attraction 



spin angle {rad) frequency [HZ] 


Figure 92 Ordered Torque and harmonics at Point 5 due to Capsule Attraction 


139 













spin angle (rad) frequency (HZ) 


Figure 93 Ordered Force and harmonics at Point 6 due to Capsule Attraction 



spin angle [rad] frequency [HZ] 


Figure 94 Ordered Torque and harmonics at Point 6 due to Capsule Attraction 
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spin angle (rad) frequency (HZ) 

Figure 95 Ordered Force and harmonics at Point 7 due to Capsule Attraction 



spin angle [rad] frequency [HZ] 

Figure 96 Ordered Torque and harmonics at Point 7 due to Capsule Attraction 
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Figure 97 Ordered Force and harmonics at Point 8 due to Capsule Attraction 
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Figure 98 Ordered Torque and harmonics at Point 8 due to Capsule Attraction 
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Figure 99 Ordered Force and harmonics at Point 9 due to Capsule Attraction 
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Figure 100 Ordered Torque and harmonics at Point 9 due to Capsule Attraction 
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The frequency spectrum for the force and the torque are shown in the following 
tables. The dominant frequencies are in bold face. 


Table 3 Frequency spectrum along the capsule axis 



First Order 

Second Order 

Third Order 

F 

y 

P0,P2 

P0,P1, P2,P3 

P2, P4 

B 

P2 

P2 

P2, P3, P4 





Table 4 Frequency spectrum to the side of the capsule axis 



First Order 

Second Order 

Third Order 

D 

P0,P2 

PO, P1,P2,P3 

P1,P2,P3,P4 


P2 

P2 

P1,P2,P3 


Table 5 Frequency spectrum close to the cap 



First Order 

Second Order 

Third Order 

F r 

PI, P2 

P1.P2.P3 

P1,P2,P3,P4 

T 

X 

PI, P2 

P2 

P1,P2,P3 


The orders in the above tables are the following. Order 1- offset or inertia order 1, 
0(1/ /?’) for the force. Order 2- inertia order 2, degenerates to MacCullagh formula for a 
perfect body (or with respect to the principal axes), 0(1 /R 4 ) for the force. Order3- 
inertia order 3. 0(1 /R 5 ) for the force. 

It is important to emphasize that the relative orders depend on the mass imperfection. 
For example, we should not rush to a conclusion regarding the contribution of the offset. 
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The dominant contribution for a perfectly symmetrical body is from the second order. In 
this case the torque is in agreement with the MacCullagh formula. The spectrum 
associated with the second order can be obtained by noticing that the direction cosines, m 
and n, have period one: m -* P\ , n — * PI . Therefore, 

Tx <=(Izz- Iyy) m n -* P2 , Ty « {Ixx - Izz) n -* P\ , Tz «= (Iyy - Ixx) m-* PI. 

Note that Iyy =Izz, thus both the offset term and MacCullagh term are amplified by 
the offset. Tx is null for a perfectly symmetrical body. Ty and Tz show various 
frequencies. However, Ty and Tz acting on a perfectly symmetrical body contribute PI 
according to MacCullagh formula. Also note that order 4 may play a role as well, because 
we assumed equal principal inertia of order 2. This assumption nulls the perfect body 
contribution of the second order but it does not cancel completely the term of order 4. 

The frequency PI is dominant only at the extreme location (point 9 inside the 
capsule), due to the offset. We investigated the vicinity of this point, and found that the 
amplitude of PI decreases as we move farther from the cap. 

Concluding Remarks 

We presented an analytical and a numerical analysis of the gravitational perturbations 
acting on the test mass due to the capsule. The outcome of this work is a closed-form 
formulation for the force/torque as well as a flexible and interactive computer code. 

Our analysis predicts that the gravitational perturbations, acting on the test mass due 
to the capsule attraction, are safely within the limit required by the experiment if 
construction criteria, outlined below, are followed.. Moreover, along most of the test 
mass trajectory, the perturbations are far smaller than the limit, and their frequencies are 
different from the modulation frequency. 

The conclusions from our investigations concerning the construction of the test 
masses are: 

- The test masses should be smaller than about 10 cm and lighter than about 2 kg. 

The second-order principal moments of inertia should be all equal within 
construction tolerances, i.e., 6L/L < 10' 4 . 

There is no need for belted cylinders (as used in the STEP satellite experiment) 
for the accuracy target of our experiment. 

The mass distribution of the test masses should be as uniform as possible within 
6M/M < 10- 4 . 

The investigation carried out is essential for the definition of the tolerable sizes, 
masses, inertia characteristics, and construction accuracy of the sensing masses. 
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Release/Leveling Mechanism Development 
Introductory remarks 

The leveling and release mechanism must be able to release the detector with initial 
conditions that do not impair the accuracy of the science measurement. Once the effects 
of the higher-order mass moments on the proof masses are made negligible through the 
right shapes and sizes, the most important external perturbation is the Earth’s gravity 
gradient. Other perturbations are of course present but they can be made negligible 
through suitable thermal design, reducing the pressure inside the detector, and shielding 
the detector from magnetic disturbances. 

The noise components associated with the Earth’s gravity gradient manifest 
themselves at twice the spin frequency and (depending on the orientation of the spin axis) 
at the spin frequency. The latter are the damaging components. The gravity component 
of importance to us is the g xy where x is the spin axis and y the sensitive axis of the 
accelerometer. As shown previously, this component is proportional to the product 
sin((j»)6 x where <J> is the elevation angle of the body symmetry axis with respect to the 
horizontal plane (defined by the local gravity) and 6 X is the centering error along the spin 
axis between the CMs of the proof masses. The formulas derived previously for the 
Earth’s gravity gradient components will be utilized at the end of this section to set a 
limit on the tolerable tilt angle at release. 

Additional harmonic components come from the rotational dynamics of the 
instrument package. These harmonics are related to the inertia characteristics of the 
package and the rotational velocity errors at release. In summary, the leveling and 
release mechanisms must provide an orientation of the spin axis close to horizontal (to 
reduce the Earth’s gravity gradient component) and rotational velocity errors sufficiently 
low to avoid saturation of the accelerometer output. 

The detector has also its own elastic dynamics (as shown in a previous section) which 
is excited by the conditions at release. We can conservatively assume that the release 
will excite the elastic dynamics of the detector up to its end of scale. The experiment 
strategy is to damp the elastic oscillations through electrical dissipative forces (see later 
on for experimental results) exercised for a few seconds after release. Once the elastic 
(natural) oscillations are abated to a level well within the dynamic range of the 
instrument, the electrical dissipative forces are removed and the detector operates as a 
high-Q detector. After the natural oscillations are abated, the oscillations of the proof 
masses will be forced by the rotational dynamics of the detector during the fall. 

In order to understand the effects of initial errors at release upon the differential 
accelerometer output, we need to develop a simplified model of the accelerometer. This 
model must contain the key dynamical elements but must also have a sufficiently simple 
formulation that shows analytically the origin and frequency content of the proof masses 
acceleration. 
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Simplified dynamical model 

The differential accelerometer consists of proof masses that are about twenty times 
lighter than the mass of the instrument package. Moreover, the proof masses are 
expected to move with respect to the CM of the whole package by only microns during 
the early transient phase and by nanometers during the measurement phase of the fall. 
Consequently, the rotational dynamics of the package is unaffected (and this will be 
confirmed in follow-up analyses) by the tiny motions of the proof masses. On the 
contrary, the rotational dynamics of the package drives the motion of the proof masses. 
In summary and with a good degree of approximation, the rotational dynamics of the 
package can be assumed to be steady. The solution of the attitude motion of the rigid 
body is available in close form for a free spinning body. Consequently, we need to write 
the equation of a proof mass that is mechanically constrained to the rotating (and 
wobbling) instrument package in order to understand the origin and frequency content of 
the acceleration output of a single proof mass. The analysis can be extended to two proof 
masses. 

The general expression for the acceleration of a test mass at point P with respect to O 
in a rotating system F can be written in matrix form as: 


w 

X 

y 

■ + [to][a>] 

X 

y 

► + 2 [a)]- 

y 

► + < 



z 








= {a P -a 0 } 


(60) 


where [co] is the angular rate matrix of the reference system and {x y z} T the position of 
the point P with respect to the center of the reference system O (placed at the instrument 
package CM). 

For a rigid body the equations are projected onto the axes attached to F itself (body 
axes) where the position vector is assumed fixed with respect to the center of F. 

Consequently, eqn. (60) transforms into: 

B 

= { a P~ a o} B (61) 


Where the angular rate matrix projected onto the body axes can be computed from the 
rotation matrix that relates the body axes to an inertial system: 
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After substituting eqn. (62) into (61), we derive the well-known acceleration matrixes 
that provide the acceleration gradient at the point P with respect to O as observed by the 
rotating observer (or equivalently a proof mass of the detector): 
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Euler acceleration matrix 


Consequently, the (apparent) acceleration gradient tensor (in body axes) between a 
point P and the package CM at O is: 

[A] = [C] + [E] (63) 

In order to compute the expression of to x , (D y , o 2 we need to consider the dynamics of 
the instrument package that houses the detector. In other words, we must solve the Euler 
equations for the torque-free case of a rigid body. We will consider an inertially 
axisymmetric body and regard x as the longitudinal axis of inertial symmetry. The non- 
symmetric body also has a known solution but it is more complicated as it involves 
elliptic integrals. 

After setting v = component of the angular rate along the axis of symmetry, the well- 
known solution of the Euler equations is obtained as follows: 


ft) Jr =V 

i <o v = -ft), cos(Qr) ( 64 ) 

ft). = w, sin(Qr) 


in which the origin of t is at the time when a) y attains its maximum value and where: 


Q = 
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body precession rate 
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a) t = -JcUy 2 + w z 2 = -JcOyo 2 + o) z0 2 transverse angular velocity 

The quantities a> y0 and co z0 are the initial components of the angular velocity 
orthogonal to the symmetry axis. In our case, they can be interpreted as the rotational 
velocity errors at release of the instrument package. 

After substituting eqn. (64) into the matrices [C] and [E], we obtain the component of 
the overall acceleration gradient matrix: 


a xy = a yx = -ft), cos(Qt)(Q - v) 
a x; = = -ft), sin(Qt)(Q - v) 

a yy = ^-ft), 2 [l - cos(2Qr)] + v 2 

1 'y 

a yz = a zy = --ft), sin(2Qr) 
a„ = •^■ft), 2 [l + cos(2Qf)] + v 2 

The acceleration vector a measured at P and projected onto the body axes, is simply: 

W = MW (65) 

where {6} = [6 X , 6 y , 6J T is the position vector from O to P (we have changed the notation 
to the 6s to highlight the fact that we are dealing with very small distances). The 
component of the acceleration along the sensitive axis (ie, z-axis in this formulation) of 
the accelerometer is as follows: 


a. = -t«),(v -Q)cos(Qr)5 A . + —ft), sin(2Qr)6 v + 


v 2 + ^-ft), 2 (l + cos(2Qr)) 




( 66 ) 


This acceleration is the dominant acceleration experienced by the proof mass once the 
natural oscillations have been abated. In fact the amplitude of the residual natural 
oscillations can be made orders of magnitude smaller than the magnitude of the 6 vector 
through initial damping. Equation (66) highlights several important points as follows: (1) 
the rotational velocity errors at release (encapsulated into oa,) combine with the 
displacement error 6 to produce an acceleration output along the sensitive axis; and (2) 
this acceleration component is modulated by the body precession rate and not by the spin 
frequency. This fact was observed previously in the results of the general elastic model 
of the detector. The conclusion is very important because it implies that the errors at 
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release do not impact directly a possible EP signal violation (that appears at the spin 
frequency). The overall advantage of the gyroscopic body (with non-spherical inertia 
ellipsoid) is that the precession rate is different from the spin rate. Moreover, the 
precession rate can be chosen at will by selecting the ratio of the moment of inertia of the 
instrument package in order to maximize the ability of extracting the signal from the 
dynamics-related noise. Since both co p and 2 appear in the expression of the 
acceleration, it is advisable to make w p non commensurate with v (the spin rate). The 
selection of the moments of inertia ratio (and hence to p ) and the ability to extract a signal 
with a strength at the threshold sensitivity of the detector will be carried out in next year 
analysis. 


Derivation of requirements 


The first requirement derives directly from the analysis carried out previously and it is 
related to the strength of the components of the Earth’s gravity gradient at the spin 
frequency. The resulting acceleration component [from eqn. (43)] is 


<5a c < 3-^-cos(0)sin(0)sin(vr)<5j +--^cos 2 (0)[sin(2vr)<5 v 




+ -^r — cos 2 (0)-l 
R> 2 


(67) 


for an instrument package spinning about x, a sensitive axis along z, and <)> the elevation 
angle of the symmetry axis with respect to the horizontal plane defined by the local 
(gravity) vertical. Equation (67) highlights the three components produced by the Earth’s 
field at the spin frequency, at twice the spin frequency, and a dc term. The strength of the 
gravity gradient component at the spin frequency must be substantially less (say 10 15 g) 
than the signal at threshold. A product <|>6 X of 0.1 deg-micron will safely meet the 
previous condition. Consequently, we can either be more relaxed on the leveling/release 
mechanism (i.e., (|> < 1 deg and 6 X < 0.1 micron) or on the centering between the two 
proof masses along the spin axis (i.e., <j> < 0.1 deg and 6 X < 1 micron). Both options will 
be kept open for the time being because they involve several technical considerations 
related to the mechanization of the leveling/release mechanism and the calibration of the 
differential accelerometer. A choice between the two options will be searched for after a 
dynamic analysis of the release mechanism and laboratory tests on the differential 
accelerometer prototype to assess the difficulties involved with the accurate centering of 
the proof masses along the spin axis. 

The requirement on the rotational velocity errors at release can be readily obtained from 
the following considerations. First, the rotational velocity error must be sufficiently 
small not to saturate the instrument output. In this case we do not have to worry about 
the instrument sensitivity because with appropriate inertia characteristics there will be no 
components at the spin frequency related to the precession dynamics of the instrument 
package. In this case the stronger component is associated with the first term of eqn. (66). 
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A rotational velocity of 1 deg/s will only produce a signal of order 10' 9 g which is well 
within the dynamic range for the expected values of spin and precession frequencies. 

However, we must also worry about the amplitude of the precession (i.e., the nutation 
angle) because through coning of the body axis the accelerometer will sense the Earth’s 
gravity gradient. The frequency sensed through coning is not exactly the spin frequency 
because the coning adds a frequency modulation at the precession frequency. 
Nevertheless, we will assume conservatively that the amplitude of the nutation angle 
must be less than previously indicated for the gravity gradient related noise. For an 
axisymmetric body the amplitude of the body-axis coning (i.e., the nutation angle) is as 
follows: 


tan(0)«^ 


l a V 


( 68 ) 


where w t is the transverse angular velocity at release, I a and I, are the moments of inertia 
about the transverse and the symmetry axis (x-axis in our case), respectively. If we 
assume an upper bound for the nutation angle of 1 deg, a spin frequency of 0.5 Hz and a 
(worst-case) inertia ratio I/I,, « 3, we obtain u>, < 1 deg/s. Likewise, if we had assumed a 
maximum 0 < 0.1 deg then we would have obtained o), < 0.1 deg/s. In summary, we will 
establish the following two sets of requirements 


a) 6 X < 1 jim (centering error between proof masses CMs along the spin axis) 

4» < 0. 1 deg (horizontality of detector) 

to, < 0.1 deg/s (angular rate error at release) 

b) 6 X < 0. 1 //m (centering error between proof masses CMs along the spin axis) 

<)> < 1 deg (horizontality of detector) 

w t < 0.1 deg/s (angular rate error at release) 

We will start by base lining the release and leveling mechanism according to option a 
while keeping option b open if the tight centering of the proof masses of the detector 
turns out to be feasible. 
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Thermal Analysis/Issues 


Following the procedure of the STEP project, we have estimated the pressure 
requirement inside the detector for keeping the radiometer effect and the damping due to 
residual gas to tolerable levels. 

Radiometer effect 

The tolerable pressure inside the detector is related to the maximum allowable 
thermal gradient through the radiometer effect, i.e., the acceleration produced by gas 
molecules emitted from regions with different temperatures’" 1 * 1 ”. This effect is likely to 
produce an acceleration error modulated at the signal frequency because, as the detector 
rotates, it could see (static) regions with slightly different temperatures. The acceleration 
produced by the radiometer effect is a = p (AT/AL)/(2pT) where p is the pressure, T the 
temperature and p the density of the sensing mass. For p = 2800 kg/m 3 (Aluminum) and 
an acceleration error of lfr 15 g, we could have for example: p = 10 9 Pa and a 
temperature gradient of 0.3 K/m, or a higher pressure and proportionally lower 
temperature gradient. As a result, we can place a requirement on the product of the two 
as follows p*(AT/AL) < 3xl0 10 Pa-K/m. In other words, the requirement on the pressure 
can be traded for the requirement on the temperature gradient. Early estimates of 
temperature gradient across the detector diameter are very small and, consequently, a 
higher detector pressure can be tolerated in a detector that is well thermally shielded. 

A thermal gradient across the detector also affects the resonant frequencies of the 
sensing masses through variations of: the Young's modulus of the material, the geometry 
of the torsional springs, and the moments of inertia of the sensing masses. Consequently, 
a thermal gradient changes differentially the resonance frequencies of the sensing masses 
and ultimately affects the common mode rejection factor. 

The common-mode rejection factor is related to the temperature variation AT across 
the detector as follows CMRF ■* la + apJAT where a is the thermal expansion coefficient 
and cxe the thermal coefficient of the Young's modulus. For Aluminum at low 
temperatures a « cxe, cxe *= -3.5xl0 4 K’ 1 *** ,v while CMRF is required to be s 10 -4 . 
Consequently, the tolerable temperature gradient for suspensions of the proof masses that 
are separated by many centimeters is of order 1 K/m. This value is less stringent than the 
tolerable value of the thermal gradient dictated by the radiometer effect and, 
consequently, is superseded by the previous effect. 


Thermal configuration of instrument package 

There are 3 elements that operate at different internal temperatures in the instrument 
package, namely, the detector, the preamplifier and the electronic box. The detector 
should be at the lowest possible temperature with an upper limit of 10 K. The Gallium 
Arsenide FET preamplifier can operate from low temperature (> 10 K) up to room 
temperature with slightly reduced performance at high temperature. The electronic box 
must be kept close to room temperature. The three boxes also dissipate different amounts 
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of power. The power dissipated by the detector W, is in the nW range which is 
negligible. The power dissipated by the preamplifier and the electronic box are typically 
W 2 «= 100 mW and W 3 ® 500 mW which are both non negligible. A simplified thermal 
scheme of the three elements is shown in Figure 101 . 



Figure 101 Simplified thermal scheme 

Thermal issues and design 

Thermal issues are not solved by simply meeting the thermal gradient requirement 
and lowering the pressure. The experiment has two cold parts, that is, the proof masses 
plus casing (i.e., the sensor) and the preamplifier and a component at almost room 
temperature (inside its case), that is, the electronics box which includes battery and 
various housekeeping functions. 

The sensor has negligible thermal dissipation while the preamplifier and electronic 
boxes have non-negligible thermal dissipations. We need to evaluate the following: (1) 
cooling down time for the sensor and preamplifier before launching the balloon; (2) the 
tolerable temperature drifts of the detector, the preamplifier and the electronics during the 
measurement phase after release; and (3) the tolerable temperature drifts of the same 
units before release during the check out phase. 

Experimental data was provided by our partners at IFSI on the power consumption of 
the units involved and their sensitivity to temperature variations. Not surprisingly, the 
sensor and the preamplifier have the highest temperature sensitivity (expressed as change 
of acceleration output per degree K), of order 10 s g/K and 10' 9 g/K, respectively, 
followed by the electronics which is a few orders of magnitude less sensitive. If the two 
sensitivities quoted above are taken at face value, the temperature variation over the 
measurement time should be less than 0.1 pK! In reality, the sensor and preamplifier will 
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be grounded to the cryostat temperature before release and then they will rely on their 
own thermal inertia during the 25-30 s measurement phase. The temperature drift is 
slow, with a time constant that is typically two orders of magnitude bigger than the 
measurement time. Since variations of less than pK are unrealistic to obtain, the question 
can be put as follows: how well an acceleration signal at the expected sensitivity can be 
extracted from a slow varying background (dependent on the temperature drifts) which 
increases by several orders of magnitude over the measurement time? 


Answering this question is important to define more realistic requirements for the 
thermal design, in general, and for solving point 2 outlined above, in particular. In order 
to answer this question we have investigated techniques of signal extraction from a noisy 
signal which also drifts by a very large amount due to temperature variations. 

This exercise is not meant to be a comprehensive work on the techniques that will be 
adopted for extracting the signal from noise but rather it is limited to the issue described 
previously. In summary, we have taken a sample 0.5-Hz signal with a strength of 2x10 15 
g (at 95% confidence level), added to a random noise with V30xl0' 15 g rms and to a 
double-exponential (other functions have also been tested) drifting acceleration that 
ramps up to an acceleration intensity many orders of magnitude bigger than the signal. 
The double exponential (or similar function) represents the different heating rates of the 
sensor and the preamplifier with two different time constants TC, and TC 2 where we have 
assumed TC, > TC : to account for the sensor’s bigger thermal inertia We have practiced 
a basic technique to extract the signal, without assuming any knowledge of the functional 
expression of the thermal drift, for increasing steepness of the thermal ramp (that is the 
rates of change of the two temperatures). 

Figure 102 shows a conservative situation in which the thermal-induced acceleration 
ramps up to 2x10 8 g, that is, its maximum value is 7-orders of magnitude bigger than the 
expected signal. The signal was extracted with reasonable ease by detrending the signal 
+ noise with a 6-order polynomial and computing the FFT of the detrended signal. This 
extraction process would not succeed for much stronger thermal drifts, although no 
additional effort was devoted to improving the process because the results obtained were 
already satisfactory. After considering these results, we can establish thermal drift 
requirements for the relevant units as follows: 


Sensor temperature drift: < 0.0067 K/s (i.e., ~0.2 K in 30 s) 
Preamplifier temperature drift: < 0.017 K/s (i.e., -0.5 K in 30 s) 


(69) 


For the temperature of the electronic box we can safely assume that it is sufficient to 
keep the temperature variation below a few degrees over a 30 s maximum free-fall time. 
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Figure 102 Example of extraction of signal from random noise plus a strongly-dominant 
slow-drifting acceleration noise (Temp-1 = sensor and Temp-2 = preamplifier). 

Preliminary Thermal Results 

A preliminary thermal analysis was carried out by using a finite-difference thermal 
analysis code to address the two more critical thermal issues, that is: (a) the cooling down 
of the detector and (b) the temperature rise of the preamplifier during the measurement 
phase. 
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Results of the cooling down are shown in Figure 103 under different assumptions of 
radiation and conduction through a cold strap with a conductance of 0.02 W/K 
connecting the case of the detector to the cryostat. The most relevant results are those 
shown in Fig. 103(a) while Figs. 103(b) and 103(c) represent cases that were run to 
assess the relative contributions of radiation and conduction, respectively. As shown by 
/ Fig. 103(a) the cooling time needed, with the cold strap, to reduce the temperature below 

10 K is shorter than a week. 


r 
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Figure 103 Temperatures of detector during cool down under various assumptions 

The second critical issue that we addressed is the temperature rise of the preamplifier 
during free fall (i.e., the measurement phase). The preamplifier was singled out as more 
critical, from the thermal point of view, than the detector because it has high power and 
low mass while the detector has negligible power dissipation and higher mass. We 
assumed that the preamplifier is at its desired temperature value (ie, 10-20 K) before the 
instrument package release and that a thermal shield is interposed between the 
preamplifier box, which is attached to a thermal radiator, and the detector. 
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Figure 104 Temperature of preamplifier and thermal shield during free-fall phase 


Figure 104 shows the results for different thicknesses of the thermal radiator. 
Clearly, in order to keep the temperature variation below 0.5 K during a 30-s time, a 1- 
cm (thick) X 20-cm (dia.) radiator is required. Additional thermal analyses will be 
performed to compute temperature variations of the three connected boxes and, 
consequently , define in more detail the thermal design beyond this preliminary analysis 
of critical issues. 
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Mechanical Report 


Release sequence 

The release of the package takes place in 2 stages. Initially the capsule is released from 
the balloon. Then, very soon after, the sensor package is released, in two steps, within the 
vacuum chamber. 

Capsule Release 

The combination of wind shear during ascent, the extremely low air resistance at float, 
and the limited time at float make it necessary to assume that the balloon will be rocking 
at the time of package release. Since it is required that capsule is vertical and non-rotating 
before release, the release strategy must account for this. The present concept, developed 
thus far, is to actively separate the orientation dynamics of the capsule from the behavior 
of the balloon. That is, we have placed a 3-axis gimbals between the capsule release 
mechanism and the portion of the gondola that stays with the balloon (capsule 
leveling/release mechanism). 

As the balloon rockets back and forth, and rotates, the gimbals will be driven to ensure 
that the capsule angular orientation is unchanged. The sole effect being that the capsule 
will be translating back and forth, and moving up and down slightly. During this process 
the linear motion will be tracked from the ground. In order to ensure that we release the 
capsule when it is as close to unaccelerated as possible, we will release it at one of the 
extremes of the balloon pendular swing. 

Since the sensor release does not have any orientation adjustment capability, the capsule 
release mechanism must be designed to meet all the alignment and stability requirements 
that cover the release of both the capsule and the sensor. 

Sensor Release 

The sensor release concept involves a two stage support (sensor release mechanism). The 
first stage is a simple spit that holds the sensor package through release of the capsule 
from the balloon. Once the capsule is in free-fall the sensor release mechanism pulls back 
far enough to disengage the spit, and draw it to a position that the sensor will not hit it on 
the way down. At this point the sensor is held between 6 springs, 3 on each side. These 
springs will be sized so that they can not impart more force than is allowed by the 
rotation stability requirements. Any transient sensor motion induced in the spit-support 
release will be allowed to damp out in the 6 spring. Once the instrument package is 
spinning at a stable rate, and oscillation from the initial stage of the release have damped, 
the mechanism second stage will pull back further, and the instrument will drop. The 
next stage of the work on this design will be to create a detailed simulation of the entire 
system, in order to examine its behavior in detail. 
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Gondola Structure 


Interface Flange^ 

(To Gondola Structure) 


Drive Motor 
(About Y Axis) 



Drive Motor 
(About Y Axis) 

Drive Motor 
(About X Axis) 


Release Mechanism 
( Gondo 1 a /Pay 1 oad ) 

Interface Flange 
’(To Payload Structure) 


Payload Structure 


Einstein Elevator 
3-Axes Gimbals Mechanism 


Figure 105 Schematic of leveling mechanism 



Einstein Elevator 
Instrument Package, Spin-Up 
& Release Mechanism 


Figure 106 Instrument package spin-up and release mechanism 
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Shown w/ Release Spindle 
In Engaged Position 


Shown w/ Release Spindle 
In Disengaged Position 



Einstein Elevator 
Instrument Package. Spin-Up & 
Release Mechanism 


Figure 107 Front view of release mechanism 

Stabilizing Spring (3) 



Release Spindle Engaged To 
Instrument Package 


Figure 108 Detail of release mechanism (with spindle engaged) 
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Stabilizing Spring (3) 
In Free Length Mode 



Figure 109 Detail of release mechanism with spindle disengaged and springs still 
engaged 


Mechanical/Thermal Design 

One of the key advances this year was in the area of the thermal design and increased 
detail in the thermal model. In order to assist in this area the main supporting structure for 
the sensor package had to be sized. Though most of the instrument design remains at the 
concept level, it became important to examine the design forces that controlled this 
component's dimensions. The sizing involved balancing the need to minimize the heat 
transfer, and therefore the support’s cross sectional area, while at the same time minimize 
the support deflection, or therefore the support’s structure moment. Studies of viable 
cryogenic materials suggested that a Kevlar composite would be the best material to use. 
Since we intent to use a composite material, and by the nature of the design, one with a 
thin wall, a further constraint on the design was a limitation on the allowable 
compression load. Coupling these requirements we were able to size the main shaft in 
such a way that is can support the instrument without allowing out-of-spec instrument 
deflection, while dropping thermal conductivity well below the allowable level. 

Below are the present support parameters: 

o Length = 0.12m 
o Diameter = 0.05m 
o Wall thickness = 0.0012m 
o Kt/m = 1 .9x 10 5 W*m/K, at 10°K 
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Cool down Techniques 

Though this will be covered in more detail in the thermal section, the preferred method 
for lowering the instrument to cryogenic temperature has been re-examined and changed. 
Initially we had examined using radiation alone to reduce the instrument operation 
temperature to below 10 K. It was quickly clear that this approach was inadequate. Next 
we examined having a cold strap in place prior to launching the balloon. We would then 
rely on just radiation from the cold Dewar walls to maintain the temperature. Though this 
is effective, it introduces issues related to manipulating the cold strap. 

To avoid this we have since base lined a third approach, one involving cooling the 
instrument down by flowing cooled gaseous helium (Ghe) through the vacuum chamber 
until the instrument is cool. We would then pumping out the remaining GHe from the 
inside of the Dewar. This give us a fast way to cool the instrument without having to 
manipulate a cold strap. 

Dewar Operation Techniques 

One of the larger instrument dynamics stability issues that needs to be addressed is the 
effects of the liquid cryogens during free fall: 

o the fact that they may slush during the fall 

o the fact that they will continue to boil during the fall 

o the fact that they needs to be vented, which will cause a propulsion effect. 

In order to avoid these effects we are beginning to examine novel Dewar operation 
techniques. The standard approach for maintaining temperatures at this level is a Dewar 
with a vacuum shroud and 2 cryogen shrouds, LN2 and LHe. But the LN2 is heavy. It can 
slush, it tends to boil similarly to water, and when it is vented it produces some 
propulsion, all of which will cause some dynamic reaction in the Dewar that could 
disrupt the experiment. The LHe, on the other hand is quite light and is far less likely to 
cause a dynamic disturbance. 

For these reasons we have decided to examine operating the Dewar, at least through the 
experiment phase, without LN2. The first approach that we examined was allowing the 
LN2 to completely vent prior to dropping the experiment package from the balloon. This 
is an acceptable approach, though it does produce some limitation on system timing. We 
would rather be able to begin the experiment when the system is settled, rather than 
waiting longer for the LN2 to boil off. Instead we are going to examine the possibility of 
allowing the LN2 boil off on the ground, and then flowing the boiled off LHe into LN2 
shroud. 

Dewar Layout 

We have moved the Dewar opening from the top of the chamber to the bottom. The 
initial Dewar design had the opening for inserting the instrument package into the top of 
the Dewar. This meant that all the feedthroughs into the system, as well as all the service 
feedthroughs for the Dewar itself were broken anytime the Dewar was opened. In 
addition, the load path from the balloon to most of the instrument mass goes through this 
connection. By moving the opening to the bottom, all the feedthroughs remain intact 
when loading the system, and only the instrument crash protection is supported through 
the Dewar cover (bottom plate). 
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Thermal Analysis 


The simplified model discussed previously has been further developed, particularly in the 
area of isolation of the sensor from the (relatively) warmer preamp and electronics 
module. In fact, the electronics has to be kept above a minimum temperature (for this 
model, we have assumed -65C) to prevent damage to the components, and possibly 
needs to be warmed further to a “cold-start” temperature at the initiation of experiment 
calibration. 

The design consists of a sensor supported by a Kevlar-reinforced shaft (discussed 
previously) that is suspended by the rotation fittings across the dewar diameter. We have 
assumed the electronics module and preamp are balanced disks supported by this same 
shaft, both to one side of the centrally located sensor module. Between these various 
elements we have places radiation shields to minimize views of the warmer elements to 
the sensor. The model representation is shown in Figure 111. The electronics are shown 
in red as a result of their maintained temperature; the preamp is just visible between the 
two largest disks. 



'=08. 1 | 


Temperature [K], Tine = Q. sec 


Figure 111 Thermal model representation 
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We have assumed that during the entire process (including cool down of the experiment) 
that the electronics is maintained above -65C, as illustrated by the red disk at the left end 
of the figure. 

There are two fundamental conductive paths between the various elements in the model: 
the supporting shaft and the wires connecting the electronic elements. The hollow shaft 
has been sized to support the masses with sufficient stiffness to allow the experiment to 
be balanced. The wires (currently assumed to be 4 pairs between each module) have 
been sized to provide less than 1-ohm resistance per circuit. Pure-metal wire like copper 
or silver have sharply peaked thermal conductivity in this temperature regime (10-25K) 
as well as rapidly changing electrical conductivity, so we have assumed for this analysis 
that the wires are Constantan, which has very little change in electrical characteristics 
from room temperature down to 4K, but its thermal conductivity drops by a factor of 100, 
providing effective thermal isolation at the experiment temperature. The conductors in 
the model have been sized to accurately represent the shaft and wires, and all the 
important materials are modeled with temperature-varying conductivity and specific heat. 

It has been shown in earlier studies that the experiment cannot be cooled effectively by 
radiative exchange with the dewar alone - it requires the addition of a thermal strap to 
achieve the desired starting temperature in a reasonable period of time. However, a 
mechanical strap (or straps) adds several complications to the design. The largest mass is 
the sensor itself, yet this is the least desirable location to affix a strap. We are currently 
considering an alternative approach: to flow gaseous helium over the experiment to cool 
it. This has the advantage that after starting conditions are achieved; the helium can be 
pumped out, allowing us to take advantage of the effective isolation of the mechanical 
hardware. 

We have run a study of the warm-up of the preamp and experiment with nominal power 
inputs and the previously mentioned initial condition of the electronics module at 208K 
(-65C). Using an assumed cooling of the preamp based on area ratio to the sensor and 
heat flow through the wires, the equilibrium temperature of the preamp prior to startup is 
approximately 20K. Using these initial conditions, and wires and support tube modeled 
as described above, the temperature rise of the preamp and the temperatures and 
estimated gradients in the sensor are shown in Figures 1 12 and 113. 
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Figure 112 Temperature rise of preamplifer during calibration 
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Figure 113 Temperature rise at various locations of sensor 
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This study has yielded several significant conclusions for the design of the experiment: 

• Properly sized alloy wires provide sufficient electrical conductance without 
significant impact on the thermal performance 

• The most significant thermal conductance is the shaft supporting the sensor and 
other modules. 

• Radiative heating of the sensor can be well controlled with concentric disks 
between modules. These disks, if of significant thickness, also provide thermal 
mass to slow the heating of the sensor. 

• Temperature gradients across the sensor package at the end of the preamp 
warmup phase appear to be quite small (see Fig. 1 14). 

It should be noted that while the temperature rates of change fall within the previously 
described limits (current estimate of maximum temperature in sensor is 0.0004 K/sec), 
improvements can be contemplated which would control these rates better. For example, 
the model currently has the shaft connected directly to the center of the aluminum endcap 
of the sensor module, essentially a worst-case solution. A set of Kevlar straps between 
the shaft and the sensor could provide high stiffness with lower conductance. However, 
such improvements, as well as better fidelity modeling of the inner structure of the 
sensor, require more sensor design detail. 





l Eflpera-t.ii.rg L£L Tine = 1 £00 sec 

Figure 114 Sensor temperature distribution at end of 20 min. 
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Instrument Electronics overview 


The overall communication and control system for the experiment is shown below. 
Minimum electronics beyond that already in the detector are added to provide a data 
collection and status link to electronics mounted outside the vacuum chamber. The link 
electronics on the detector is a battery-powered RS-232 to IR transceiver comprised of a 
1 10° field of view IR emitter-detector pair driven by an IR encoder-decoder and linked to 
an identical pair on the side of the chamber. The total power required for this portion of 
the interface should be less than lOOmW. Presently the electronics will need to run at -40 
°C, but we will investigate if we can go lower. 

The power and telemetry control box will be a comprised of a Book-sized PC with an 
Ethernet interface to a procured telemetry radio, such as the Freewave FGR1115RE 
900MHz spread spectrum radio. A small IR to RS-232 module provides the interface to 
the IR transceiver inside the capsule. The externally mounted system will run from a 
single +12VDC 85 W power source. The unit is expected to occupy a volume of 40cm x 
30cm x 25cm. 

The gondola will contain a standard SIP module for ground communication and 
control. We will need to provide an interface between the SIP telemetry interface and a 
radio that links to the experiment. This interface will essentially be the same design as the 
power and telemetry control box in the instrument. 

Other means of communication within the capsule were considered, including a low 
power 2.4 GHz radio link, and this option could be a viable alternative if the IR link is 
not feasible. The major advantage to this approach would that the link does not depend on 
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System Requirements Development 


Tabic 6 summarized key requirements of major subsystems derived thus far. 


Table 6 Development of Requirements vs. Design Drivers 


Design Driver 

Scienee- 

originated Reqmt. 

Engineering Requirements 

Tolerable differ, 
acceleration 

Detector/ 

Package 

Cryostat 

Drop System 

Free-fall (time > 20 s) 

N/A 

Transient 
damping time 
<5 s: 

N/A 

Low-speed 
ballistic coeff. 
Po > 5000 kg/m 2 

Amplifier noise (white) 

< 10’ 14 g/VHz 

Preamp. 
T n < 60 mK 

N/A 

N/A 

Brownian noise (white) 



LHe cryostat 

N/A 

Temperature gradients 
(Radiometer effect] 

< 10 15 g 

p( AT/Ax) < 
3x10'"' Pa-K/m 

TBD 

N/A 

Acceleration noise inside 
capsule in free fall 

< 10 12 g 

CMRF < 10 4 

Pressure inside 
capsule 
p t < 10 6 mBar 

Structural and 
attitude 
freqs. » co 

Earth’s gravity gradients 

< 10 12 g (2w) 

< 10 15 g (co) 

Centering of 
proof masses (6 X ) 
along spin axis 
within 1 pim 

N/A 

Verticality before 
release 

6 X 4> < 0.1 pm-deg 

Cryostat's gravity 
gradients 
(distributed mass) 


Centering of 
proof masses 
w ithin 10 urn 

Cryostat internal 
diameter ^ 1 m 

N/A 

Gravity gradients of lump 
masses on board capsule 

— 

Centering of 
proof masses 
w ithin 10 urn j 

N/A 

Mass-distance | 
exclusion zones 
(see pp. 53-54) 

Magnetic disturbances 

< 10 16 g 

Use Niobium 
alloy blanket 
around detector. 
Degauss proof 
masses 

Temperature of 
package 
T < Tc 
(T c = critical 
temperature) 

Limit magnetic 
moments outside 
sensor package 
M m < TBD A-m 2 
and r > TBD m 

Higher-order mass 
moments 

< 10 "'g 

Proof masses 
with almost equal 
moments of 
inertia (6M/M - 
6L/L < 10 -*) 

Mass of cryostat 
<500 kg. 
Intern, dia. ^ 1 m 

N/A 

Centrifugal gradients due 
to skewed spin axis 

< 10 15 g 

Centering of 
proof masses 
6 X < 1 //m 

N/A 

Verticality w ithin 
0.1°; rate errors 
at release 
<;0.1 °/s 


10 = signal angular frequency; to,, = detector resonant angular frequency 
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Experimental Activity (IFSI/CNR) 1 

Our partners have been working very hard in carrying out experimental activity on 
critical aspects of the detector development. The team at the Institute of Space Physics 
(IFSI/CNR) has already built a differential accelerometer prototype and carried out 
significant laboratory measurements on the prototype. 

The prototype was designed with the goal of exploring key aspects of a differential 
accelerometer and not for carrying out a preliminary test of the Equivalence Principle in 
the laboratory. As such, the prototype has two sensing masses of the same material with 
their centers of masses close together but not perfectly coincident. The prototype has 
several features in common with the instrument that we expect to develop for the flight 
experiment as follow: same capacitive pick-up system; same elastic suspension of the 
sensing masses; same measurement chain to extract the differential signal from the 
accelerometer. 

Figure 116 and Figure 117 show schematics of the differential accelerometer 
prototype. Each sensing mass (blue and orange elements in the figure) is shaped to 
compenetratc (with leeway) into the opposing one so as to bring the two CMs close 
together (perfect coincidence is not a requisite for this prototype). 



Figure 116 Exploded view of differential accelerometer prototype 


1 Section contributed by V. Iafolla S Nozzoli and M. Fiorcnza of IFSI/CNR funded through Italian Space 
Agency (ASI). 
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Figure 117 Cross section of (assembled) differential accelerometer prototype 



NJ 


Figure 118 Simplified electrical diagram of signal pick-ups 

Each sensing mass can rotate (through clastic restraints) about the pivot axis on each 
side of the accelerometer. The two pivot axes are parallel to one another. The sensitive 
axis of the differential accelerometer is perpendicular to the pivot axes and to the 
longitudinal axis of the instrument (spin axis). Each sensing mass has two fixed 
capacitor plates for signal pickup. 
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Figure 1 19 Pictures of differential accelerometer prototype 
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The output signals of the sensing masses are independently amplified and sent to the 
data acquisition system for conditioning, filtering and comparing. This prototype 
differential accelerometer builds on the heritage of (single) high-sensitivity 
accelerometers that were built at IFSI through the years. These accelerometers were 
tested extensively in the laboratory and in the field"". Key characteristics of the 
differential accelerometer prototype are summarized in Table 7. 


Table 7 Key characteristics of differential accelerometer prototype 


Item 


Sensing mass 

0.22 kg 

Quality factor 

2900 

Resonance frequency 

18.12 Hz 

Preamplifier noise temp. (K) 

0.76 

Preamplifier type 

AD743 (commercial) 



1 time (day) 


Figure 120 Luni-solar tide signal filtered with 49-hour filter (1 mas = 2.78xl0' 7 g) and 
measured with single accelerometer in August 1998. 
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The single accelerometer, in its various implementations, has demonstrated a 
remarkable sensitivity and repeatability. The accelerometer was able to resolve a 
modulated gravity signal (at 0.3 Hz) with a strength of less than 10' 9 g in the presence of 
seismic noise. Another single accelerometer is used to measure the Luni-Solar tides 
(working as a tilt meter) in the Gran Sasso laboratory (Italy). Excellent tidal signal over 
periods of months have been measured by the IFSI team. A month-long sample of a 
Luni-Solar tide in August 1998 is shown in Figure 120 as an example of the performance 
of the single accelerometer in the field. 

Damping of transient oscillations 

The differential accelerometer prototype was developed and built to test, among 
others, key aspects of differential acceleration measurements as follows: (a) abatement of 
the natural dynamics excited by the instrument release into the capsule; and (b) rejection 
of the common-mode signals. Point (a) above is critical for the success of our 
experimental scheme. For an instrument that is required to resolve acceleration of 10 14 g, 
the release is an abrupt event that pushes the detector well beyond its saturation point. 
Moreover, in a (flight) instrument with a Quality factor (Q) of order 10 5 and a resonance 
frequency of a few Hz, the transient oscillations would take a very long time to decay to 
within the instrument’s dynamic range. The strategy that we plan to follow in the 
experiment is to reduce the Q factor during the first few seconds of fall from 10 5 down to 
a few units. The reduction of the Q factor is accomplished by inserting a resistance in the 
feedback control loops of the accelerometer (see Fig. 121). The technique has been 
tested successfully on the prototype accelerometer. 



Figure 121 Electrical diagram of one test mass pickup system with resistance added to v 

the feedback loop 

Figure 122 shows the oscillation amplitude (from the oscilloscope output) of the \ 

accelerometer after an excitation. The accelerometer has a resonant frequency of 18.5 Hz 
and an (undamped) Q, = 2900. The introduction of a 50 Mohm resistive load in the 
feedback loop, reduces the Q from its undamped value to Q,’ = 441 and, consequently. 
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the oscillation amplitude decays much more rapidly. Subsequently, the resistance is 
removed to demonstrate that this operation does not re-excite the oscillation amplitude. 



Figure 122 Oscillation amplitude decay vs. time with and without resistance inserted in 
the feedback loop 


It is worth reminding that the decay of the proof mass oscillations is driven by the 
total quality factor of the electro-mechanical system which is as follows: 


1 1 1 . . I , conRC 

— = + — with — = A ^ 

Qt Qm Qe Qc 1 + (iOqRC) 


(70) 


where k is the electro-mechanical coupling factor 


A = 


CE 1 2 
mw o 


that is the ratio of electrical energy to mechanical energy of the oscillator. The electro- 
mechanical coupling factor of the instrument prototype is a low A = 0.01 . The flight 
instrument will have a value much closer to unity thanks to a lower value of u) (l and 
higher value of the capacitance C. Consequently, an electrical resistance of the value 
adopted for this test coupled into the feed-back loop of the flight instrument will imply a 
reduction of the value of Q, to a few units. In summary, this damping technique is able to 
provide the desired damping performance. 
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Common-mode rejection factor 

One important characteristic of a differential accelerometer is its ability to reject 
perturbations that are not differential, i.e., common-mode disturbances. This ability is 
quantified by the common-mode rejection factor (CMRF). 



Figure 123 Experimental setup to perturb differential accelerometer with periodic 
common-mode disturbances. 
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The IFSI laboratory has built a rotating system whose rotation axis can be accurately 
tilted off the vertical (see Fig. 123). The accelerometer is mounted inside the rotating 
device with the sensitive axis orthogonal to the spin axis. The rotating outfit is tilted by 
about 10 3 rad away from the vertical and spun at a slow rotation with a frequency of 0.15 
Hz. As the differential accelerometer slowly rotates, it measures a small component 
(proportional to the tilt) of the Earth’s gravity field at the rotation frequency. This 
perturbation is a common-mode disturbance which affects equally the two proof masses. 



i 


» 



Figure 124 Accelerometer outputs: (a) single acceleration from proof mass 1 and 2 and 
(b) differential acceleration. 


The differential output generated by the calibration signal is then treated by the 
software through spectral analysis of the individual signals to adjust the proportionality 
factors of the two sensing masses and the phase differences of their responses. 
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Figure 125 Spectra of individual and differential acceleration outputs: (a) after 
amplitude calibration only and (b) after amplitude and phase calibration. 

Figure 125 shows that after calibrating for amplitude and phase a 10 4 attenuation is 
readily obtained for the differential signal. This level of attenuation is effective not only 
at the perturbation frequency of 0.15 Hz but also over a larger frequency band. An 
attenuation of 10 4 or equivalently a common-mode rejection factor of 10‘ 4 meets the 
present requirement on the CMRF for the proposed tests of the Equivalence Principle. 
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Concluding Remarks 

The laboratory activity consisted in the construction of a laboratory prototype of a 
differential accelerometer. The laboratory prototype has been used to conduct key tests 
on the differential instrument. The team at IFSI demonstrated the ability of the detector 
to damp quickly transient oscillations by utilizing a resistive load in the feedback loops 
and then removing that load to reestablish a high quality factor of the detector. A rotating 
device with tilt control was also built at IFSI. This device was utilized to impart (through 
the Earth’s gravity) common-mode perturbations to the differential accelerometer. These 
calibration disturbances have been used to trim the acceleration outputs of the individual 
proof masses in order to obtain a common-mode rejection factor better than 10' 4 in a 
sufficiently large frequency band centered at the spin frequency. 
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